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SUMMARY 


Two-phase flow instability (Ledinegg Excursion) and critical heat flux in heated 
microchannels are of great concern in the design and operation of numerous practical 
systems. In order to prevent Ledinegg flow excursion and eventual burnout of the heated 
channel, limits are typically imposed on the channel power and coolant mass flux to 
avoid the possibility of operation within the negatively-sloped region of the coolant 
channel demand curve. 

The primary purpose of this investigation was to experimentally determine the 
effect of design and operational parameters on the onset of flow instability (OFI) and 
critical heat flux (CHF) in thin, horizontal, uniformly-heated annuli. The work was 
motivated by the need for such data for the design and safety analysis of the Accelerator 
Production of Tritium (APT) project under development at Los Alamos National 
Laboratory. To this end, thirteen different experimental setups were employed to 
determine the OFI and CHF behavior of annular channels with hydraulic diameters on the 
order of two millimeters (annular gap widths of approximately 1.0 mm). The OFI data 
were then compiled and used to create two correlations; one comparing OFI heat flux to 
saturation heat flux and the other comparing OFI mass flux to saturation mass flux. The 
CHF data acquired in this investigation were compared with several empirical CHF 
correlations available in the literature and used to demonstrate that OFI, rather than CHF, 
is the limiting phenomenon in microchannel heat transfer. 

All OFI experiments were conducted in the Georgia Tech MicroChannel Test 





Facility (GTMTF) with test section inlet temperatures ranging from 30 to 70°C, exit 
pressures from 0.345 to 1.034 MPa, surface heat fluxes from 0.45 to 2.5 MW/m 2 , inner- 
to-outer surface heat flux ratios from 0:1 to 1:0, and fully saturated non-condensable gas 
concentrations. The corresponding range of coolant mass flux was 85 to 2900 kg/m 2 s. 
A total of 66 OFI experiments were conducted in this investigation; together with the OFI 
data (129 experiments) of Blasick (1999), these data were used to develop the following 
correlations. 

q" 0FI = 0.942V m ,± 9.80% (S.l) 

G on = 1.064 * G sat ± 11.55% (S.2) 

The above correlations mean that for a specified mass flux, the heat flux necessary to 
cause OFI equals 94.2% of the heat flux required to cause bulk exit saturation of the 
working fluid. Alternatively, for a specified heat flux, the mass flux corresponding to 
OFI equals 1.064 times the mass flux required to cause bulk exit saturation of the 
working fluid. The test section heated length and inner-to-outer heat flux ratio were 
found to have insignificant effects on OFI. 

All CHF experiments were also conducted in the Georgia Tech MicroChannel 
Test Facility (GTMTF) with test section inlet temperatures ranging from 30 to 65°C, exit 
pressures from 0.345 to 1.034 MPa, coolant mass fluxes from 100 to 400 kg/m 2 s, an 
inner-to-outer surface heat flux ratio of 1:1, and fully saturated non-condensable gas 

xviii 






concentrations. The corresponding range of surface heat flux was 0.2309 to 1.0686 
MW/m 2 . A total of 35 CHF experiments were conducted in this investigation and 
compared with several CHF correlations available in the literature. In addition, the CHF 
data were compared with the OFI correlation developed during this investigation; the 
following relation was obtained. 

fcHF = 3.047 * q" OF[ ±l2.63% (S.3) 

Given a specified mass flux, the heat flux necessary to cause CHF was found to be 3.047 
times the heat flux required to cause OFI. This effectively demonstrates that OFI, rather 
than CHF, is the limiting phenomenon in microchannel heat transfer. The data obtained 
in this investigation provided the basis for quantifying the thermal margins for the 
proposed Accelerator Production of Tritium (APT) target bundle geometry under both 
normal operation and accident conditions. 
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CHAPTER I 


INTRODUCTION 


Forced convection using single-phase liquid water is a common cooling 
mechanism used in several high power density engineering applications. Typically, 
design specifications require coolant channels to remain completely filled with sub¬ 
cooled liquid water to ensure adequate cooling of heated surfaces. Despite these planned 
steady state parameters, accident or transient conditions can cause boiling to take place. 
For this reason, it is important in the design and operation of heated channels to 
understand and predict the onset of boiling and two-phase flow instabilities that can lead 
to channel dry-out and potential overheating of the system. The present work describes 
experiments that examine the onset of Ledinegg flow instability and critical heat flux 
caused by significant voiding in an annular channel. 


1.1 Onset of Flow Instability (OFI) 

Figure 1.1 shows a typical demand curve for a uniformly heated channel. In the 
high mass flux regime, sub-cooled single-phase conditions exist throughout the channel. 
Consequently, all heat transfer is accomplished through single-phase forced convection. 
As the mass flux through the channel is decreased, sub-cooled boiling followed by 


1 



Subcooled 



Figure 1.1: Schematic Diagram of the Demand Curve for a Heated Channel 

saturated boiling, and finally single-phase convection to superheated vapor takes place 
(Kennedy, et al., 1997). 

In Figure 1.1, the terms ONB, OSV, OFI, and S refer to the onset of nucleate 
boiling, onset of significant void, onset of flow instability, and saturation conditions, 
respectively. Prior to the OSV point, the coolant flow is primarily single-phase, sub¬ 
cooled liquid with some bubbles sticking to the walls of the heated channel. These 
bubbles are unable to withstand condensation and hence collapse when released from the 
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wall and mix with the bulk fluid flow. At the onset of significant void, the void fraction 
is approximately five percent and bubbles released from the channel walls survive 
condensation and become entrained in the bulk fluid flow; the flow becomes two-phase. 
Flow boiling ensues following the OSV point and a fully developed nucleate boiling 
regime dominates. For high mass fluxes, the onset of significant void and the initiation 
of bulk boiling occur at approximately the same point. As the mass flux through the 
heated channel decreases to the point where pressure drop across the channel is at a 
minimum on the demand curve (referred to as the onset of flow instability point), a flow 
excursion can occur resulting in a significant reduction in flow through the channel 
(Kennedy, et al., 1997). 

Referring to the demand curve shown in Figure 1.2, case A represents the 
situation that would exist if a positive displacement pump were used to circulate coolant 
through a heated channel (Lahey and Moody, 1979). A positive displacement pump 
operating at a given speed would force coolant through the boiling channel at a fixed rate 
regardless of the overall pressure drop. For this reason, the channel could be operated at 
position 1 by maintaining a constant pump speed. Case B represents the situation where 
numerous channels operate in parallel with the boiling channel. For this case, operating 
the system steadily and continuously at position 1 is impossible. A small negative 
perturbation in the system would immediately shift conditions to position 2 while a small 
positive perturbation would shift the flow state to postion 3. A flow excursion from 
position 1 to position 2 would result in a significant reduction in flow rate within the 
heated channel and may lead to dryout or physical burnout of the channel. Case C 
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represents the situation that would result if a centrifugal pump were used to circulate the 
coolant through the heated channel. Once again, steady and continuous operation at point 
1 is impossible since a small perturbation could result in a flow excursion to either 
position 2' or position 3'. Finally, case D represents a situation where the circulation 
pump has a steep characteristic so that position 1 could be maintained and Ledinegg 
instability would not occur (Kennedy, et al., 1997). 


A 



Mass Flow Ra*te, m 

Figure 1.2: Excursive Instability (from Lahey and Moody, 1979) 
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In order to prevent Ledinegg flow excursions in practical systems, limits are 
imposed on the channel power input to exclude the possibility of operation within the 
negatively-sloped segment of the demand curve at all anticipated flow rates. Typically, 
these power limits are conservatively estimated by using OSV models and/or 
correlations, where OSV is considered a precursor to the onset of flow instability (OFI) 
which, in turn, may be viewed as a precursor to the critical heat flux (CHF) (Kennedy, et 
al., 1997). In this research, the location of the OFI point was compared with commonly 
used OFI and OSV correlations for large diameter channels. In addition, CHF 
experiments were performed in order to demonstrate that OFI, rather than CHF, is the 
limiting phenomenon in microchannel heat transfer. 


1.2 Critical Heat Flux 

Internal flow critical heat flux (CHF) is one of the most important design 
limitations in many industrial systems. At CHF, a severe reduction in heat transfer 
between the heated surface and coolant occurs, frequently referred to as forced 
convection burnout (Roach, 1998). In situations where the surface heat flux remains 
constant, CHF is marked by an abrupt increase in the wall temperature of the heated 
surface. Forced flow CHF is a complicated phenomenon in which thermal and 
hydrodynamic processes are tightly coupled, but the basic aspects can be outlined as 
shown in Figure 1.3 (Stoddard, et al., 2000). 


5 



Figure 1.3 shows the variation of coolant flow pattern alongside a plot of the 
typical behavior of the wall and liquid tempeatatures. In this simplified example, single¬ 
phase liquid enters a heated, vertical tube at the bottom and proceeds to boil as it travels 
through the tube, exiting as superheated vapor. Nucleate boiling begins when the tube 
wall temperature reaches the saturation temperature of the coolant. As long as liquid 
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Figure 1.3: Typical Temperature and Flow Regime Variation in a Heated Tube 
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coolant remains in contact to some degree with the tube wall, the wall temperature 
remains essentially constant, until the dryout point is reached. Beyond this point, only 
vapor is in contact with the tube wall, and its temperature increases dramatically 
(Stoddard, et al., 2000). 

As mentioned previously, the complicated coupling of thermal and hydrodynamic 
effects makes it difficult to generalize or predict the exact progression of the flow pattern 
as the channel geometry changes. Furthermore, the exact mechanism that leads to CHF 
can vary depending on the exact nature of the flow, thus not all of the flow regimes 
shown in Figure 1.3 necessarily occur. For this reason, a large variety of data for 
different channel geometries and coolant conditions have been generated by past 
researchers, and several different mechanistic models have been developed for the 
prediction of CHF, with separate models addressing the various observed mechanisms. 
The CHF data collected in this investigation will be compared with existing CHF models 
available in the literature; a new CHF correlation will not be developed. As mentioned 
previously, the fundamental goal of this research is to demonstrate that OFI, rather than 
CHF, is the limiting phenomenon in microchannel heat transfer. 


1.3 Accelerator Production of Tritium (APT! Project 

The Accelerator Production of Tritium (APT) Project is one of the alternatives 
under consideration by the Department of Energy (DOE) for the production of the tritium 
necessary to replenish the United States nuclear weapons stockpile. Tritium decays at a 
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rate of 5.5 percent per year and therefore must be produced regularly. All previous 
tritium production occurred in several fission reactors located at the DOE Savannah River 
Site in South Carolina, which have since been shut down (Blasick, 1999). 

Tritium would be produced in the APT facility by irradiating a target assembly 
with a high-energy proton beam to create neutrons; He-3 gas would then be used as the 
tritium producing feedstock. The proposed APT target design employs D 2 0 (heavy 
water) cooled concentric tungsten annuli as the primary neutron source. This tungsten 
neutron source (TNS) consists of thirteen ladders, each separated by a 0.3 m void region 
as shown in Figure 1.4. Each ladder has several short rungs consisting of six concentric, 
Inconel-718 clad tungsten cylinders of varying thickness since the proton beam produces 
the highest heat fluxes on the front ladders. Hence, the tungsten cylinder wall thickness 
increases from 2.030 millimeters in the front ladder to 7.703 millimeters in the back. 
This causes an attendant increase in the rung tube diameter from 43.4 millimeters to 77.7 
millimeters. The number of rungs in each ladder assembly ranges from 28 in the front 
two ladders to 16 in the back ladders, and the rungs are offset from ladder to ladder to 
eliminate proton streaming. Each rung has a total length of 25 centimeters while the 
beam spot covers the central 19 centimeters (Pasamehmetoglu, et al., 1998). A cross- 
section of a rung from ladders one and two, the peak power ladders of the APT target, is 
shown in Figure 1.5. 
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Figure 1.4: Schematic of APT Target Assembly (from Pasamehmetoglu, et al., 1998) 


Table 1.1: Summary of Geometric Data for Ladders and Rungs in APT Target 


Dimension 

Ladders 

1 -2 

Ladders 

3-4 

Ladders 

5-6 

Ladders 

7-8 

Ladders 

9-13 

Number of Rungs 

28 

26 

22 

19 

16 

Rung Outer Diameter (cm) 

4.343 

4.679 

5.563 

6.477 

7.696 

Rung Pitch (cm) 

7.478 





Tungsten Volume Fraction 

58.5% 

■ 

Whk'JMk 



Number of Clad Tungsten Tubes per Rung 

6 

6 

6 

6 

6 

Clad Thickness (cm) 

0.0127 

0.0127 

0.0127 

0.0127 

0.0127 

Center Hole Diameter (cm) 

0.2032 

0.2032 


0.2032 

wmm 

Tungsten Ring Thickness (cm) 


lillillB 



mmm 

Total Ring Thickness (cm) 





KH 

Flow Gap (cm) 

0.1016 

0.1016 

0.1016 

0.1016 

0.1016 

Rung Tube Thickness (cm) 

0.0890 

0.0890 

0.0890 

0.0890 

0.0890 
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Figure 1.5: Scale Cross-Sectional View of a TNS Rung for Ladder Levels One or Two 


The resultant annular flow regions have a gap width of 1.016 millimeters. The 
inside and outside diameters are 6.6 and 8.6 millimeters, respectively, for the smallest 
proposed annular region in the peak power rung. The dimensions of the innermost 
annular channel are identical in all rungs of each ladder. An illustration of the innermost 
annular channel, the focus of this investigation, is shown in Figure 1.6. In addition, Table 
1.1 provides a summary of the geometric data for all of the ladders and rungs in the APT 
target as provided by Pasamehmetoglu, et al. (1998). 
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Figure 1.6: Simplified View of the Innermost Annular Coolant Channel 


Coolant water would be provided to the horizontal assemblies through larger 
diameter vertical inlet and exit plenums as shown in Figure 1.4. For this reason, each 
rung of each ladder would be connected in parallel with the other rungs of its ladder. In 
addition, all ladders would be connected in parallel via the inlet and exit coolant 
manifolds shown in Figure 1.4. 
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1.4 Motivation and Objectives 


As a result of the potentially severe consequences of Ledinegg flow excursion and 
critical heat flux, several experimental and theoretical studies have been performed by the 
nuclear industry to quantify the effect of various design and operational parameters on 
these limits. The existing experimental data and correlations are well established for 
channels with hydraulic diameters larger than four millimeters. The current APT design 
has interior coolant channels with a hydraulic diameter of roughly two millimeters. Little 
data exists in the literature for either OFI or CHF in such small diameter channels. To 
this end, Blasick (1999), conducted an investigation using light water (H 2 O) to determine 
the effect of various parameters on OFI for a simulated APT innermost annular coolant 
channel. The work by Blasick (1999) covers the low heat flux range corresponding to 
operating conditions in off-peak APT ladder assemblies. This investigation will 
complement the work of Blasick (1999) by determining the effect of various parameters 
on OFI at high heat flux conditions corresponding to the peak power ladder assemblies in 
a simulated APT innermost annular coolant channel using light water. Additionally, 
experiments will be conducted at low mass flux conditions to determine the effect of 
various parameters on CHF. 

The primary purpose of this investigation was to experimentally determine the 
effect of system parameters on OFI and CHF in heated microchannels. To this end, 
thirteen different experimental setups were employed to determine the OFI and CHF 
behavior of annular channels with hydraulic diameters on the order of two millimeters 
(annular gap of approximately 1.0 mm). The OFI data were then compiled and used to 
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create two correlations; one comparing OFI heat flux to saturation heat flux and the other 
comparing OFI mass flux to saturation mass flux. The CHF data acquired in this 
investigation were compared with several empirical CHF correlations available in the 
literature and used to demonstrate that OFI, rather than CHF, is the limiting phenomenon 
in microchannel heat transfer. No new correlations were developed for CHF in 
microchannels. In conclusion, the data obtained in this investigation were used to 
quantify the thermal margins for the proposed APT target bundle geometry under both 
normal operation and accident conditions. 

The remainder of this thesis is organized as follows. Chapter II provides a brief 
overview of the literature available on OFI and CHF in small diameter channels. Chapter 

III gives a detailed description of the instrumentation and equipment used in all 

/ 

experiments performed during this investigation. Chapter IV outlines the experimental 
procedures used to setup and perform all OFI and CHF experiments. All OFI and CHF 
experimental results are reported and compared with available models in Chapter V. 
Finally, Chapter VI summarizes the conclusions of this investigation and provides 
recommendations for future microchannel heat transfer research. 
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CHAPTER II 


LITERATURE REVIEW 


II. 1 Introduction 

Two phase flow instability can lead to serious safety problems in the operation of 
heated channels supporting forced and natural circulatory flows. For this reason, it is 
important in the design of heated channels to understand and predict the onset of boiling 
and two phase flow instabilities that can lead to channel dry-out and potential overheating 
of the system. Reviews of flow instability and critical heat flux have been published by 
McBeth, et al. (1964), Bowring (1972), Boure, et al. (1973), Ishii (1976), Bergles (1977), 
Yadigaroglu (1981), Groenenveld, et al. (1986), Bankoff, et al. (1991), and Caira, et al. 
(1993). The following is a brief overview of previous work on OFI and CHF. In 
addition, note that the majority of the OFI section of this chapter was derived from the 
literature review performed by Kennedy, et al. (1997) while the bulk of the CHF section 
of this chapter was adapted from Roach (1998). 

II.2 Onset of Flow Instability 

Two phase flow instabilities are generally categorized as either dynamic or static 
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(Boure, et al., 1973). Dynamic instabilities (e.g., flow regime oscillations) involve 
transient inertial dynamic and feedback effects and lead to oscillatory flows, which can 
be analyzed using the transient flow conservation equations. Modeling system behavior 
during dynamic oscillations, however, requires numerical solution of nonlinear 
conservation equations when perturbations are imposed on the flow (Gurgenci, et al., 
1983). Dynamic instabilities occur over a broad range of conditions and become serious 
safety problems when their amplitudes become large. 

Static instabilities result from pressure dropBflow rate characteristics of heated 
channels and can be analyzed using steady state flow conservation equations. The typical 
demand curve (steady state pressure drop versus flow rate) of a boiling channel has an 
inverted S-shape (Figure 1.1). The different two phase flow regimes for heated channels 
are shown in Figure 1.3. As described in Chapter I, when the channel is part of a forced 
circulation loop, the segment of the demand curve with negative slope can be unstable 
depending on pump characteristics. If the mass flow rate through the channel drops 
below the local minimum point on the curve, then a flow excursion can occur resulting in 
a significant reduction in the flow rate through the channel. Such flow excursions can 
lead to physical burnout of the heated channel (Roach, 1998). 

The onset of flow instability (OFI) is known to occur at a slightly lower flow rate 
than the flow rate corresponding to the onset of significant void. The onset of significant 
void (OSV) occurs when a further reduction in flow rate results in a sharp increase in 
coolant void fraction in the channel. Consequently, the pressure drop due to flow 
acceleration resulting from increased voidage renders the gradient of the demand curve 
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negative. OSV can therefore be considered a conservative criterion for avoiding classic 
Ledinegg flow excursions (Yadigaroglu, 1981). Furthermore, OSV signifies the 
initiation of bulk boiling in the subcooled flow boiling regime. If significant void 
fraction develops in a nuclear reactor, then reactor power may be affected (Kennedy, et 
al., 1997). 

In the past, OSV and void fraction have been extensively studied in channels 
experiencing subcooled boiling (Whittle and Forgan, 1967; Evangelisti and Lupoli, 1969; 
Chan, 1984; Rogers, et ah, 1987; Johnston, 1989; Ahmad, 1970; Dix, 1971; Levy, 1967; 
Staub, 1968; Unal, 1975). Most of the published studies, however, deal with OSV and 
OFI in relatively large channels. In this investigation, the onset of flow instability in a 
horizontal, annular microchannel subject to a uniform and constant wall heat flux is 
experimentally investigated. The efforts are motivated by the need to assess the thermal 
limits for the concentric, annular target coolant channels in the proposed Accelerator 
Production of Tritium (APT) system. The hydraulic diameter of these coolant channels is 
significantly smaller than the parameter ranges covered by the existing models and data 
(Kennedy, et ah, 1997). 

OSV, OFI, and other related phenomena may occur in research reactors where 
flow channels typically have small hydraulic diameters. However, systematic studies 
dealing with OSV and OFI in microchannels are scarce. Inasaka, et ah (1989) measured 
subcooled boiling pressure drop in channels with diameter of 1.0 mm and larger. They 
noted that the onset of nucleate boiling (ONB) point was well predicted by the correlation 
of Bergles and Rhosenow (1964), while the correlations of Saha and Zuber (1974) and 
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Levy (1967) both failed to predict OSV data well in small channels. Their measured two 
phase friction factor multipliers for capillaries, however, were considerably smaller than 
the values predicted by the Lockhart and Martinelli (1949) correlation. Inasaka, et al. 
(1989) attribute the latter discrepancy to the possible overprediction of quality in 
capillaries by Ahmad’s correlation (Ahmad, 1970). Hainoun, et al. (1996) recently 
reported a two fluid model capable of predicting OFI phenomenon in small channels. 
Their model applies the aforementioned correlation of Saha and Zuber (1974) for the 
prediction of OSV. 


II.2.1 Flow Instability Initiation 

During a hypothetical loss of coolant accident, the coolant flow rate through the 
heated channels decreases. Under normal operating conditions, the coolant remains 
subcooled and thus single phase. However, if a substantial drop in flow rate occurs, the 
onset of nucleate boiling (ONB) may ensue. Bankoff, et al. (1991) presented a 
comprehensive review of ONB and the mechanics of bubble formation and departure. At 
the ONB point, bubbles begin to form on the walls of the heated channel. Due to the 
subcooled status of the bulk coolant, these bubbles either collapse or slide along the 
channel walls; they do not separate and survive in the bulk flow. The presence of the 
bubbles on the walls causes the differential pressure across the test section to increase. In 
turn, the friction factor rises due to the increased differential pressure across the channel. 
When vapor generation becomes great enough, the bubbles detach from the walls and 
survive in the subcooled, bulk coolant flow. This is known as the onset of significant 
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void (OSV). Vapor bubbles in the bulk flow act to lower the average density of the fluid 
within the channel. The increased apparent friction factor coupled with the increased 
acceleration term causes the differential pressure across the channel to rise. In parallel 
channels with an imposed external pressure difference, the flow can easily become 
unstable and possibly dry out due to a flow excursion (Bankoff, et al., 1991; Lahey and 
Moody, 1979). 

In practice, limits are imposed on channel power input over all anticipated flow 
rates to avoid operation within the negative sloped region of the demand curve. This 
measure prevents Ledinegg flow excursions during both normal and accident conditions. 
At low flow rates, the power limits dictated by the flow excursion are considerably lower 
than those corresponding to critical heat flux (CHF). Therefore, OFI is the limiting 
criterion for establishing the maximum power level for safe operation of the system. For 
high mass flux, however, CHF may become more limiting than OFI. Generally, the 
design and safety analysis of practical systems require that both criteria be satisfied with 
an adequate safety margin corresponding to a designer specified confidence level 
(Kennedy, et al., 1997). 


II.2.2 Onset of Significant Void 

Many models and correlations exist which predict OSV. Some of these models 
are based on theoretical arguments while others are based on empirical analysis. In 
general, two phase flow phenomena are difficult to model theoretically which explains 
the development and use of empirical correlations. Examples of commonly accepted 
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formulations include those reported by Bowring (1962), Saha and Zuber (1974), and 
Unal(1975). Some correlations such as Levy (1967), Staub (1968), and Rogers, et al. 
(1987) incorporate a theoretical model with empirical adjustments. 


II.3 Review of OSV Models 

In this section, various theoretical and empirical OSV models published in the 
literature are reviewed; Bankoff, et al. (1991), reported a comprehensive examination of 
the subject. The more widely reported correlations and their ranges of applicability are 
summarized in Tables 2.1 through 2.10, while Tables 2.11 through 2.13 summarize 
recent OFI correlations produced at Georgia Tech by Kennedy, et al. (1997), Roach 
(1998), and Blasick (1999) in conjunction with the APT project. 


II.3.1 Theoretical OSV Models 

Theoretical OSV models are based on vapor bubble mechanics in vertical, 
upward, subcooled flow. The mechanics involved in analyzing bubble behavior focus on 
the heat transfer and fluid flow at the bubble-liquid interface. Bubble departure is 
determined by performing a force balance. Three aspects of the bubble-liquid interface 
combine to determine when the bubble will depart the channel wall: buoyancy of the 
bubble within the liquid, shear force exerted on the bubble by the flowing liquid, and the 
surface tension which holds the bubble to the wall (Kennedy, et al., 1997). 
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Of the three OSV correlations described by Bankoff, et al. (1991), the one 
proposed by Levy (1967) is the simplest. Summarized in Table 2.1, the Levy correlation 
is valid for water under the following conditions: mass flux range of 130 to 1420 kg/m 2 s, 
operating pressure between 0.41 and 13.6 MPa, and heat flux from 0.24 to 1.91 MW/m 2 . 

The model of Staub (1968) is described in Table 2.2 and is applicable for two 
channel geometries: circular with hydraulic diameter between 0.006 and 0.20 meter and 
rectangular with hydraulic diameter between 0.004 and 0.066 meter. The correlation is 
valid for both water and R-22, but no range is specified for the applicable heat fluxes. 

The correlation of Rogers, et al. (1987) is more specific and therefore more 
restricted than those of Levy and Staub. As tabulated in Table 2.3, this correlation is only 
valid for annular tube geometry with three specific hydraulic diameters: 0.009, 0.012, and 
0.017 meter. Also, this correlation only holds for tube lengths of 0.48 and 0.60 meter. 
The correlation is valid for water under the following conditions: mass flux range of 70 to 
450 kg/m 2 s, operating pressure between 1.50 and 1.55 MPa, and heat flux from 0.3 to 1.2 
MW/m 2 . 

Even though the analytical models of OSV are similar, it is interesting to note that 
different bubble shapes are assumed. Levy assumed a spherical shaped bubble, Staub 
assumed a hemispherical bubble that took on the contour of the channel, and Rogers, et 
al. assume that the bubble was a Atruncated sphere, with a contact angle at the surface 
equal to the equilibrium contact angle (Bankoff, 1991). 
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II.3.2 Empirical OSV Models 

Often it is easier to use empirical correlations as opposed to theoretical models for 
OSV prediction due to the complex mechanisms that exist in two phase flow. Several 
empirical models are summarized in Tables 2.4 through 2.13. Note that for each of these 
correlations, the liquid subcooling at the OSV point can be expressed as a function of 
pressure, hydraulic diameter, heated length, heat flux, mass flux, inlet temperature, and 
fluid properties at designated temperature and pressure. For example, a pressure term is 
used in the Bowring (1962) correlation while other correlations only consider the effect 
of pressure on fluid properties. It is also interesting to note that Bowring (1962) and 
Thom, et al. (1966) are the only dimensional correlations; Saha and Zuber (1974), Unal 
(1975), Ahmad (1970), and Sekoguchi, et al. (1974) are all dimensionless. In fact, the 
Saha and Zuber (1974) correlation includes both the Stanton and Nusselt numbers. This 
signifies that under low mass flux conditions, flow behavior is thermally controlled, 
while under high mass flux conditions, flow behavior is hydrodynamically controlled. 
The Unal (1975) correlation suggests a similar split into two different regimes that 
depend on flow velocity. 

The theoretical models described above have been compared with the empirical 
models. Both the theoretical Levy (1967) correlation and the empirical Bowring (1962) 
correlation predict that as liquid flow rate through the channel increases, the degree of 
subcooling at the OSV point decreases. Similarly, the theoretical Staub (1968) 
correlation was compared with the empirical models of Costa (1967), Griffith, et al. 
(1968), Christensen (1961), and Kroeger and Zuber (1968). The experimental data as 
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compared to the predictions of the Staub (1968) correlation were 15% higher for circular 
channels and 15% lower for rectangular channels. 

All models reviewed did not deal with microchannels. The size range was 
generally limited to channels with hydraulic diameter greater than 4 millimeters. Peng 
and Wang (1993) have provided evidence indicating significant differences between the 
characteristics of boiling microchannels and larger channels. The aim of this 
investigation was to examine the effect of microchannel geometry and operating 
conditions on OFI. Comparisons were made between experimental OFI data and 
predictions of the more widely used OSV models. 
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Table 2.1: Levy’s Theoretical OSV Model (1967) 


Author: 

Levy (1967) 

Model Type: 

Theoretical, Dimensionless 

Tube Shape: 

No Restriction 

Tube Diameter: 

No Restriction 

Tube Length: 

No Restriction 

Tube Orientation: 

Vertical 

Flow Direction: 

Upwards 

Medium: 

Water 

Mass Flux: 

130 to 1420 (kg/m 2 s) 

Pressure: 

0.41 to 13.6 (MPa) 

Heat Flux: 

0.24 to 1.91 (MW/m 2 ) 


Model 


Dimensional Distance to Tip of Bubble: 

y; = 0 . 015 ^^- 
A/ 

Heat Transfer Coefficient: 

h = 0.023—Re^ 8 Pr“ 
D h ’ ' 


Wall Shear and Friction Factor: 

f £ 10 6 

/ = 0.0055 1+ 20,000— +- 

l D k Re f 



— = 10" 4 
D h 
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Table 2.1 (continued) 


Heat Flux Parameter: 



Conditions for Onset of Bubble Departure: 

0 < Yb + < 5 T s -T fd =^-Q?r f Y; 

5 < Yb + <30 T s -T f d =^-5e|p r/ + ln 1 + Pr r 

Y B + > 30 T s - T f d = y - 5 < 2 [ + ln[l + Pr r ] 


l 5 

+ 0.5 In 
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Table 2.2: Staub’s Theoretical OSV Model (1968) 


Author: 
Model Type: 


Staub (1968) 

Theoretical, Dimensionless 


Tube Shape: 

Tube Diameter: 
Tube Length: 
Tube Orientation: 
Flow Direction: 


Circular, Rectangular 

0.006 to 0.020 (m), 0.004 to 0.066 (m) 

Unknown 

Vertical 

Upwards 


Medium: 
Flow Rate: 
Pressure: 
Heat Flux: 


Water, R-22 

0.5 to 8.0 (m/s), 0.1 to 0.9 (m/s) 
0.17 to 6.8 (MPa) 

Unknown 


Model 


Bubble Diameter: 


fp,V}D, Dl(p f -p,) 

32 gqf(p) \2qf(fl) 

f{j3) = 0.02 for V f < 0.45 (m/s) 
f{j3 ) = 0.03 for V f > 0.45 (m/s) 


British Units 


Dimensionless Distance to Tip of Bubble: 


y;=o.5d b 




Wall Shear and Friction Factor: 


P f Vj \ 2 D 
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Table 2.2 (continued) 




Table 2.3: Rogers, et al. Theoretical OSV Model (1987) 


Author: 

Rogers, et al. (1987) 

Model Type: 

Theoretical, Dimensionless 

Tube Shape: 

Annulus 

Tube Diameter: 

D h = 0.009, 0.0012, 0.0017 (m) 

Tube Length: 

0.48, 0.60 (m) 

Tube Orientation: 

Vertical 

Flow Direction: 

Upwards 

Medium: 

Water 

Mass Flux: 

70 to 450 (kg/m 2 s) 

Pressure: 

1.50 to 1.55 (MPa) 

Heat Flux: 

0.3 to 1.2 (MW/m 2 ) 


Model 


Dimensionless Distance to Tip of Bubble: 


Where: 


A = 



(l + cos£ 0 ) 

_4 7V c , g 


c , = 2 + 3cos# 0 - cos 3 0 O 
c 2 = tt-0 o + cos 0 O sin 0 Q 
c 3 = sin 0 Q [cos(<9 0 -10) - cos fo+10)] 


u 


Q 

c d 


58 

= 0 o +5 
= 1.22 
24 
Re fl 


+ 0.14 

for 20 < Res < 40 
for 4 < Re B < 20 
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Table 2.3 (continued) 


Re 5 


PfUr&j) 

M f 


Calculation of u r : 

y; =0.57; (l + cos 6> 0 ) 


Y r + <5.0 
5.0 < Y r + < 30 
Y r + > 30 


u ;=y; 

u; = 5.0 In Y r + -3.08 
M r + 2.51n7 r + +5.5 


Wall Shear and Friction Factor: 

f - = 0.046 Re -020 

PfV} 


Calculation of Tb + : 
Y b + < 5.0 

5.0 < Y b + < 30 
Y b + > 30 


r;=Pr,r; 


t;= 5 


n =5i 


Pr, + lml + Pr 


i. 


i 


£^__i 




j 


Pr 7 + ln(l + 5 Pr 7 )+ 5.0 In 


y + 

30 


Conditions for OSV: 


Where: 



Nu = 0.023 Re™ Pr™ 

F r =1.11 + -0.036 for 0 O = 30° 
F r = 1.06 + -0.022 for 0 O = 80° 
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Table 2.4: Bowring’s Empirical OSV Correlation (1962) 


Author: 

Model Type: 

Tube Shape: 

Tube Diameter: 
Tube Length: 
Tube Orientation: 
Flow Direction: 

Medium: 

Mass Flux: 
Pressure: 

Heat Flux: 


Bowring (1962) 
Empirical, Dimensional 


Water 

Unknown 

1.1 to 13.8 (MPa) 

Unknown 


Unknown 

Unknown 

Unknown 

Unknown 

Unknown 


Model 


Where: 


r t 

T,-T„=ny 

¥ f 

7 = 10” 6 (l4.0 + / > ) 
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Table 2.5: Thom, et al. Empirical OSV Correlation (1966) 


Author: 

Thom, et al. (1966) 

Model Type: 

Empirical, Dimensional 

Tube Shape: 

Unknown 

Tube Diameter: 

Unknown 

Tube Length: 

Unknown 

Tube Orientation: 

Vertical 

Flow Direction: 

Upwards 

Medium: 

Water 

Mass Flux: 

Unknown 

Pressure: 

5.14 to 6.80 (MPa) 

Heat Flux: 

Unknown 


Model 


T. ~ T U= O' 02 *, -FT 
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Table 2.6: Ahmad’s Empirical OSV Correlation (1970) 


Author: 

Model Type: 

Tube Shape: 

Tube Diameter: 
Tube Length: 
Tube Orientation: 
Flow Direction: 

Medium: 

Mass Flux: 
Pressure: 

Heat Flux: 


Ahmad (1970) 

Empirical, Dimensionless 

Circular, Rectangular 
0.004 to 0.007 (m) 

Unknown 

Vertical 

Upwards 

Water 

Less than 8000 (kg/m 2 s) 

0.17 to 0.5 (MPa), 8.3 and 11.0 (MPa) 
Unknown 

Model 


h = 2.44-^-Re 0 , 5 Pr^ 

A 7 7 
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Table 2.7: Dix’s Empirical OSV Correlation (1971) 


Author: Dix (1971) 

Model Type: Empirical, Dimensionless 

Tube Shape: Annulus 

Tube Diameter: 0.018 (m) 

Tube Length: 1.07 (m) 

Tube Orientation: Vertical 

Flow Direction: Upwards 

Medium: R-114 

Mass Flux: 100 to 210 (kg/m 2 s) 

Pressure: 0.32 and 0.85 (MPa) 

Heat Flux: 0.004 to 0.03 (MW/m 2 ) 

Model 

T s -T f d =0.00135^Re^ 
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Table 2.8: Saha and Zuber Empirical OSV Correlation (1974) 


Author: 

Saha and Zuber (1974) 

Model Type: 

Empirical, Dimensionless 

Tube Shape: 

Circular, Rectangular, Annulus 

Tube Diameter: 

0.01 to 0.24 (m), 0.004 to 0.012 (m), 0.012 to 0.026 (m) 

Tube Length: 

Unknown 

Tube Orientation: 

No Restriction 

Flow Direction: 

No Restriction 

Medium: 

Water, R-22,R-114 

Mass Flux: 

Unknown 

Pressure: 

0.1 to 13.8 (MPa), 1.1 to 13.8 (MPa), 0.32, 0.85 (MPa) 

Heat Flux: 

Unknown 


Model 


Nu = 



for Pe < 70,000 


St = -t— -r = 0.0065 


for Pe > 70,000 


Where: 


Nu GD k C 


Pe = — = RePr = 
St 


h^pf 
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Table 2.9: Sekoguchi, et al. Empirical OSV Correlation (1974) 


Author: 

Sekoguchi, et al. (1974) 

Model Type: 

Empirical, Dimensionless 

Tube Shape: 

Annulus 

Tube Diameter: 

0.002 to 0.003 (m) 

Tube Length: 

0.7 to 4.5 (m) 

Tube Orientation: 

Vertical 

Flow Direction: 

Upwards 

Medium: 

Water 

Mass Flux: 

310 to 2100 (kg/m 2 s) 

Pressure: 

0.14 to 1.6 (MPa) 

Heat Flux: 

0.05 to 1.75 (MW/m 2 ) 


Model 


. 0.65 


T s ~T f , d = 13.5 c 


fg 


Pf 


K h fg°J 
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Table 2.10: Unal’s Empirical OSV Correlation (1975) 


Author: Unal(1975) 

Model Type: Empirical, Dimensionless 


Tube Shape: 

Tube Diameter: 
Tube Length: 
Tube Orientation: 
Flow Direction: 


Circular, Rectangular, Annulus 
0.004 to 0.020 (m) 

Unknown 
No Restriction 
No Restriction 


Medium: 
Mass Flux: 
Pressure: 
Heat Flux: 


Water, R-22 

132 to 2818 (kg/m 2 s), 180 to 1391 (kg/m 2 s) 
0.1 to 15.8 (MPa), 1.2 and 3.3 (MPa) 

0.15 to 1.92 (MW/m 2 ), 0.02 to 0.06 (MW/m 2 ) 

Model 


n Vs 

q" 

f 

= a 

(a = constant) 

For Water: 

a = 0.24 

(Vf > 0.45 m/s) 


a = 0.11 

(Vf < 0.45 m/s) 

For R-22: 

OO 

o 

II 

a 

(V f > 0.45 m/s) 


a = 0.11 

(V f < 0.45 m/s) 


Heat Transfer Coefficient: 


hD 


k f 

hD„ 


h = 0.023Re™Pr° 4 


= 0.020 


(D V ' 53 




Re™ Pr™ 


for circular and rectangular tubes 


for annulus 
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Table 2.11: Kennedy, et al. Empirical OFI Correlation (1997) 


Author: Kennedy, et al. (1997) 

Model Type: Empirical 


Tube Shape: 

Tube Diameter: 
Tube Length: 
Tube Orientation: 


Circular and Prototypical 
1.131 to 1.445 (mm) 

0.16 (m) 

Horizontal 


Medium: 
Mass Flux: 
Pressure: 
Heat Flux: 


Water 

800 to 4500 (kg/m 2 s) 
0.344 to 1.034 (MPa) 
1.0 to 4.0 (MW/m 2 ) 


Model 


(a =10.0±8.7°C 

= 0.90 f ",.,±6.9% 

G 0 „=1.11G„±6.3% 
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Table 2.12: Roach Empirical OFI Correlation (1998) 


Author: 

Model Type: 

Tube Shape: 

Tube Diameter: 
Tube Length: 
Tube Orientation: 


Roach(1998) 
Empirical 

Prototypical 
1.131 to 1.445 (mm) 
0.16 (m) 

Horizontal 


Medium: 
Mass Flux: 
Pressure: 
Heat Flux: 


Water 

125 to 950 (kg/m 2 s) 
0.344 to 1.034 (MPa) 
0.1 to 0.5 (MW/m 2 ) 


Model 


Fully-Degassed: 




G o „=0.89G„,±10% 


Air-Saturated: 

?V, = 1.07?"„, ±10% 
G on =0.91 G„, ±10% 
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Table 2.13: Blasick, Empirical OFI Correlation (1999) 


Author: 

Model Type: 

Tube Shape: 

Tube Diameter: 
Tube Length: 
Tube Orientation: 

Medium: 

Mass Flux: 
Pressure: 

Heat Flux: 


Blasick (1999) 
Empirical 

Annulus 

1.448 to 2.002 (mm) 
0.19 (m) 

Horizontal 

Water 

85 to 1428 (kg/m 2 s) 
0.344 to 1.034 (MPa) 
0.124 to 1.0 (MW/m 2 ) 


Model 


Q"ofi = 0.900#" ia ,±l 1.3% 
G 0 fi — l-086G ia , ±9.90% 
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II.4 Critical Heat Flux 


Forced flow critical heat flux (CHF) is an important design limitation for 
industrial systems. Recent experimental studies have provided understanding of the basic 
phenomenology of various modes of CHF. In vertical, upward flow channels, CHF may 
occur in three different patterns. These basic patterns are illustrated in Figure 2.1. First, 
at high void fraction, annular-dispersed flow is common and ultimately results in channel 
film dryout (Diagrams “C” & “D” of Figure 2.1). Second, when the flow is in the slug or 
plug regime, the dryout of the liquid film layer between the slug and the wall leads to 
localized CHF events (Diagram “B” of Figure 2.1). Third, in high mass flux systems, 
departure from nucleate boiling will occur at low quality. Coalescence of bubbles on the 
heated surface leads to local critical void fraction which, in turn, leads to the evaporation 
of the liquid layer between the bubble region and the heated surface (Diagram “A” of 
Figure 2.1). For high mass flux systems using highly subcooled coolant, CHF is 
generally attributed to departure from nucleate boiling (DNB) by the mechanism just 
described. It has been determined that the critical heat flux will actually increase as 
channel diameter decreases. The formation of bubbles at the channel walls will be 
hampered in channels with very small diameters (i.e., microchannels), therefore, CHF 
will occur at higher heat fluxes. In this study, the coolant mass fluxes will be extremely 
low. It is therefore expected that channel dryout will be the primary cause of CHF. For 
CHF driven by dryout, the effect of diameter will be inverted (Roach, 1998). 
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II.5 Review of CHF Models 


In this section, various empirical and tabular CHF models published in the 
literature are reviewed. The more widely reported correlations and their ranges of 
applicability are reported in Tables 2.14 through 2.17. 


II.5.1 Empirical CHF Correlations 

The thermal-hydraulics group of the APT design team used the Bowring (1972) 
correlation to predict CHF in the peak power coolant channels. This correlation is ideal 
for the APT target coolant channels in that its heated length, pressure, and mass flux 
ranges cover the design parameter ranges under investigation. The extremely broad 
pressure (0.2 to 19.0 MPa) and mass flux (136 to 18,600 kg/m 2 s) ranges make this one of 
the most useful CHF correlations in the literature (Todreas and Kazimi, 1990). The 
Bowring (1972) correlation along with its ranges of applicability is listed in Table 2.14. 

The Caira (1993) correlation was developed with data from uniformly heated 
circular, rectangular, and annular vertical upflow channels. The correlation is applicable 
for 0.1 to 8.4 MPa pressure, 2000 to 9000 kg/m 2 s coolant mass flux, 0.3 to 25.4 mm 
diameter, and 10 to 61 cm heated length. In predicting CHF, the Caira (1993) correlation 
incorporates three dimensionless parameters strictly dependent on hydraulic diameter and 
mass flux. This correlation along with its ranges of applicability is listed in Table 2.15. 

The McBeth (1964) correlation was developed from a compilation of burnout data 
derived from uniformly heated round tubes with vertical upflow. There are two relations 


41 




predicting CHF, one for low velocity flow and one for high velocity flow. The minimum 
value between these two relations is the predicted CHF value. The correlation is valid 
over a pressure range of 0.1 to 12.1 MPa, a diameter range of 1 to 37.5 mm, and a heated 
length range of 2.54 to 366 cm (VIPRE - 01 Reference Manual). This correlation along 
with its ranges of applicability is listed in Table 2.16 

The CISE-4 correlation is a modification of CISE-3 in which the critical flow 
quality (x c ) approaches unity as the mass flux decreases to zero. This correlation is based 
on the quality-boiling length concept and is typically restricted to boiling water reactor 
applications. The CISE correlation can be applied to rod bundles by including the ratio 
of heated to wetted perimeters as shown in Table 2.17 (Todreas and Kazimi, 1990). 


II.5.2 Tabular CHF Correlation 

The Groeneveld (1986) correlation uses a dimensionless table to predict CHF in 
uniformly heated 8 mm diameter round channels. This value can then be converted to a 
CHF value for channels 2 to 16 mm in diameter. The dimensionless tabular value of 
CHF is determined by entering in the exit quality and two dimensionless quantities that 
are based on geometry and fluid properties (e.g. density, surface tension, diameter, and 
mass flux). This correlation involves a tremendous amount of interpolation and 
extrapolation, and hence is subject to inherent inaccuracies. The dimensionless 
parameters associated with the Groeneveld (1986) CHF correlation are given in Table 


2.18. 





Table 2.14: Bowring’s Empirical CHF Correlation (1972) 


Author: 

Model Type: 

Tube Shape: 

Tube Diameter: 
Tube Length: 
Tube Orientation: 

Medium: 

Mass Flux: 
Pressure: 

Heat Flux: 


Bowring (1972) 
Empirical 

Circular 

0.002 to 0.045 (m) 

0.15 to 3.7 (m) 
Horizontal 

Water 

136 to 18,600 (kg/m 2 s) 
0.2 to 19.0 (MPa) 
Unknown 


Model 


A-DGh f x 

a" =- — 

q CHF AC 


2.317 


r h r DG^ 


A = 


‘fg 


\ 4 J 


1 + 0.0143F,gVd 


_ 0.077 F,DG 


1 + 0.347F. 


G 


Y 


ti = 2.0-0.5P* 


P R = 0.145P 


Fx 


.1356, 

_ P R ex p[20.89(l-Pj+ 0.917] 


F, = 


1.917 


1.309F, 


F, = 


P ff 17023 exp[l6.658(1 - P fi )+ 0.667] 


1.667 


P| 316 exp[2.444(l -P R )+ 0.309] 


F 4 = F 3 P’ 


.649 
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Table 2.15: Caira, et al. Empirical CHF Correlation (1993) 


Author: 

Caira, et al. (1993) 

Model Type: 

Empirical 

Tube Shape: 

Circular and Annular 

Tube Diameter: 

0.03 to 25.4 (mm) 

Tube Length: 

10 to 61 (cm) 

Tube Orientation: 

Vertical 

Flow Direction: 

Upwards 

Medium: 

Water 

Mass Flux: 

2000 to 9000 (kg/m 2 s) 

Pressure: 

0.1 to 8.4 (MPa) 

Heat Flux: 

Unknown 

Model 

h f -h)J>e 
qcHF ~ l+ru»- 

0 = y o D y 'G y2 

0 = y 4 D ys G y6 

Q = y 7 D y *G y » 


y 0 = 10829.55 

y 3 = 0.978 

= 0.188 

y, = -0.0547 

y 4 = 0.188 

<N 

II 

£ 

y 2 = 0.713 

ys = 0.486 

y9 = 0.36 


y 6 = 0.462 

o 

II 

O 

VO 

i —»• 
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Table 2.16: McBeth, et al. Empirical CHF Correlation (1964) 


Author: 

McBeth, et al. (1964) 

Model Type: 

Empirical 

Tube Shape: 

Circular 

Tube Diameter: 

0.04 to 1.475 (in) 

Tube Length: 

1.0 to 144.0 (in) 

Tube Orientation: 

Vertical 

Flow Direction: 

Upwards 

Medium: 

Water 

Mass Flux: 

0.0073 to 13.7 (Mlbm/ft 2 hr) 

Pressure: 

15 to 1750 (psia) 

Heat Flux: 

Unknown 


Model 


f Q ^°' 51 


q-CHF ( loW - vd) = 0.00633*. D- t °‘ — (l - x) 


A-0.25D. 


Q chf vel) — ■ 


10 


V 10 J 


J I h / s x 


c 


q" C HF = min b"c//F (low - vel}, q" (high-vel)] 


A'=y 0 D?G y ' 
C'= y 3 D?G ys 


y 0 = 106.5 
y, = 0.847 
y 2 = 0.677 


y 3 = 60.3 
y 4 = 1.4 
y 5 = 0.937 
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Table 2.17: CISE-4 CHF Correlation (1975) 


Author: 

Model Type: 

Tube Shape: 

Tube Diameter: 
Tube Length: 
Tube Orientation: 
Flow Direction: 

Medium: 

Mass Flux: 
Pressure: 

Heat Flux: 


Gaspari, et al. (1975) 
Empirical 

Circular 

0.0102 to 0.0198 (m) 
0.76 to 3.66 (m) 
Vertical 
Upwards 

Water 

1085 to 4069 (kg/m 2 s) 
4.96 to 6.89 (MPa) 
Unknown 


Model 


= 


A 

D 


e \ 


L cr +b . 


1 


a = 


1 + 0.0001481 


f pX 
1-— G 

v ^ 


l-*- 


a - 


( G XA 


1000 , 


G* = 3375 


V P.J 


( 


Z? = 0.199 


\ 0.4 


-1 GD 


1.4 


V 


if G < G* 


if G > G* 


46 



Table 2.18: Groeneveld, et al. Tabular CHF Correlation (1986) 


Author: 

Groeneveld, et al. (1986) 

Model Type: 

Tabular 

Tube Shape: 

Circular 

Tube Diameter: 

2.0 to 16.0 (mm) 

Tube Length: 

Unknown 

Tube Orientation: 

Vertical 

Flow Direction: 

Upward 

Medium: 

Water 

Mass Flux: 

0 to 7500 (kg/m 2 s) 

Pressure: 

0.10 to 20.0 (MPa) 

Heat Flux: 

Unknown 


Model 


Look-up Table to determine CHF for an 8 mm diameter circular tube. 




GD 


0.5 


0.5 _ 0.5 

C7 p f 


CHF d = 


CHF ; 


8 mm 
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CHAPTER III 


INSTRUMENTATION AND EQUIPMENT 


All experiments in this investigation were performed using the Georgia Tech 
MicroChannel Test Facility (GTMTF). The GTMTF is a general purpose test facility 
designed to examine single and two-phase flow phenomena in microchannels. The 
facility was previously used by Kennedy, et al. (1997) and Roach (1998) to study OFI 
and CHF in microchannels with circular cross sections and with cross sections simulating 
coolant channels in a micro-rod bundle. Blasick (1999) later modified the GTMTF to 
accommodate test sections with annular flow passages simulating the innermost coolant 
channels in the APT target. This chapter describes the additional modifications made to 
the GTMTF as well as the test sections and power supplies used to perform high heat flux 
OFI and low mass flux CHF experiments for thin, concentric annuli. 

III.l Experimental Test Facility Overview 
The GTMTF design criteria, as elucidated by Kennedy, et al. (1997), Roach 
(1998), and Blasick (1999), were as follows. 
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1) To provide a stable, controllable, and measurable flow of water through the test 
section over a broad range of mass fluxes. 

2) To provide for a flexible system configuration and the possibility of running 
experiments with and without a test section bypass line. 

3) To allow for control and measurement of the dissolved non-condensable gas 
content within the system. 

4) To allow for fine control and measurement of the test section inlet temperature. 

5) To allow for fine control and measurement of the test section exit pressure. 

6) To provide stable, controllable, and measurable electrical power for heating of the 
test section. 

The remaining sections of this chapter provide details on the GTMTF hardware, 
instrumentation, and data recording setups used in this investigation to accomplish these 
goals. Also, the exact design criteria and assembly specifications are detailed for each 
test section configuration. A total of thirteen experimental setups were used to perform 
all of the experiments reported in this investigation: twelve OFI setups and one CHF 
setup. The experimental setup number was incremented each time the test section or any 
component of the GTMTF was modified. Consequently, many experimental setups 
utilize several of the same components. In order to simplify documentation of 
experimental setups, each component used in an experiment was given a descriptive 
name and number designation that can be cross-referenced with a figure and matching 
description as outlined in Section III.2. 
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III.2 Experimental Setups 


Each experiment run in this investigation utilized a uniquely definable 
experimental setup. As described above, a total of thirteen different experimental setups 
were employed in this investigation. Each setup was comprised of the following six 
separate, independently alterable configurations. 

1) Flow Loop Configuration 

2) Inner Heater Configuration 

3) Outer Heater Configuration 

4) Mounting Configuration 

5) Insulation Configuration 

6) Non-Condensable Gas Status 

Each time one or more of these system configurations were altered, the setup number was 
sequentially increased. Table 3.1 outlines the thirteen setups used and their 
corresponding system configurations. 


Table 3.1: Summary of Experimental Setups 


Setup 

Flow Loop 
Configuration 

Inner Heater 
Configuration 

Outer Heater 
Configuration 

Mounting 

Configuration 

Insulation 

Configuration 

Non-Condensable 

Gas Status 

1 

2 

i 

1 

1 

None 

Air Saturated 

2 

2 

1 

2 

1 

1 

Air Saturated 

3 

2 

i 

3 

1 

None 

Air Saturated 

4 

1 

2 

4 

1 

None 

Air Saturated 

5 

1 

3 

4 

1 

1 

Air Saturated 

6 

3 

3 

4 

i 

1 


7 

3 

3 

4 

1 

None 

Air Saturated (LP)* 

8 

3 

4 

5 

1 

1 

Nitrogen Saturated 

9 

3 

5 

5 

1 

i 


10 

3 

6 

6 

2 


Nitrogen Saturated 

11 

3 

7 

7 

2 

■EEH 

Nitrogen Saturated 

12 

3 

8 

7 

2 

None 

Nitrogen Saturated 

13 

4 

9 

8 

2 

None 

Nitrogen Saturated 


* The non-condensable gas status for Setup 7 was saturated with air at atmospheric temperature and 
pressure. The system was pressurized using a bellows type accumulator, so no new gas was forced into the 
working fluid during system pressurization. 
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III.3 Flow Loop Configuration 


A total of four flow loop configurations were used during this investigation. The 
flow loop was modified, depending on the experiment, in order to achieve the necessaiy 
flow rate control for a specific set of flow parameters. Before describing the flow paths 
employed by the different flow loop configurations, the mechanisms used by all four flow 
loop configurations, regardless of flow path, must be detailed. Figure 3.1 shows a 
schematic diagram of the GTMTF components and instrumentation used in all four flow 
loop configurations and Table 3.2 describes the significant hardware components in the 
system. 


III.3.1 System Instrumentation 

Table 3.3 provides a list of measuring instruments used in the GTMTF along with 
their manufacturer, model number, range, and accuracy. 

III.3.1.1 Absolute Pressure Transducers 

Four Rosemount model 1144A pressure transducers (Components 5 of Figure 3.1) 
were utilized in the GTMTF to measure test section inlet and exit pressure and pump inlet 
and exit pressure. All four transducers were mounted below the respective points where 
the pressure was to be measured in order to avoid potential incorrect readings caused by 
entrained air bubbles lodged in the transducer. 
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Figure 3.1 : Schematic Diagram of the Georgia Tech MicroChannel Test Facility 



COOLANT FLOW 

















Table 3.2: Listing of GTMTF Hardware Components 



Description 

Manufacturer 

Model 

A 

Multi-stage centrifugal pump 

Teel 

2PC45 

B 

Variable speed positive displacement pump 
(magnetically-coupled gear) 

MicroPump 

185 


Variable speed positive displacement pump 
(progressing cavity) 


Seepex 


0015-24MD 


E 

Pump bypass valve 

Whitey 

SS-12NBS8-SH 

F 

System throttling valve 

Watts 

CF8M 

G 

In-line heater 

In house construction 

- 

H 

Low flow measurement line metering valve 

Nupro 

SS-4LA 

I 

High flow measurement line isolation valve 

Watts 

CF8M 

J 

High flow measurement line throttling valve 

Watts 

CF8M 

K 

High flow measurement line metering valve (coarse) 

Whitey 

SS-3NBF4 


High flow measurement line metering valve (fine) 


High flow line / Test section line isolation valve _ 

High flow line / Test section bypass line isolation valve 


Test section heat exchanger (tube-in-tube) _ 

Throttling valve 


Bleed valve _ 

System heat exchanger (tube-in-tube) _ 

Saturation tank pressure regulation line isolation valve 


Accumulator pressure regulation line isolation valve _ 

Three-way valve 


Oxygen sensor housing _ 

Degassing tank isolation valve _ 

Degassing tank 


Degassing tank gas separation column _ 

Degassing tank gas bleed valve 


Variable-speed positive displacement pump (peristaltic) 
System water fill valve 


System pressure relief valve 


Degassing tank heat exchanger _ 

Degassing tank heater (cartridge style) 


Saturation tank _ 

Saturation tank liquid separation column 


Saturation tank gas bleed valve (integrated with flowmeter) 


Saturation tank compressed gas supply cylinder _ 

Saturation tank gas fill valve 


Saturation tank blow-off valve _ 

Saturation tank bubbler select valve 


Accumulator (bellows type, Teflon construction) 




Dv - 


Parker Hannifin 
Teel 


_ Hoke 

Parker Hannifin 
Whitey 


_Teel_ 

Whitey 


In house construction 

_Teel_ 

In house construction 


In house construction 
Whitey 


Masterflex 
Jones Mfg. Co. 


Parker 


In house construction 
Watlow 


In house construction 
In house construction 


Omega 


1315G4B 


CF8M 


DHTC-SS-8 

A743-CF8M 


1513G4S 

DHTC-CU-6 

SS-3NBF4 


A743-CF8M 

SS-45XS8 


A743-CF8M 


SS-IRS8 


7018-20 
150 WSP 


2030 CBBI-1 


5-65 C 


2030 CBBI-1 


MIMIlliKIBI 


H3020T 


OO 

Accumulator compressed gas supply cylinder 

- 

- 

PP 

Accumulator bleed valve 

Teel 

A743-CF8M 

QQ 

Accumulator depressurization valve 

Teel 

A743-CF8M 

RR 

Saturation tank recirculation pump 
(magnetically coupled gear) 

MicroPump 

000-655 

SS 

Saturation loop inlet valve 

Whitey 

SS-3NBF4 

TT 

Saturation loop exit valve 

Whitey 

IRS4 
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Table 3.3: Listing of GTMTF Instrumentation 


Label 

Description 

Manufacturer 

Model 

Serial Number 

1 

Rotameter, low range 

Fischer & Porter 

10A6132N 

8712A0041A1 

2 

Rotameter, high range 

Brooks Instruments 

1305D08A1A1A 

8603HC013 

3 

Liquid Flow Transducer 

EG&G/FTI 

FTO-3/RI-51 

F096107517RI5189010 

4 

Liquid Flow Transducer 

EG&G/FTI 

FT-8/RI-51 

804950RI5189030042 

5 

Pressure Transducer 

Rosemount 

1144A-0200A22 

Multiple 

6 

Thermocouple, Type-E 

Omega 

EMQSS-125G-6 

None 

7 

Pressure Gauge, Bourdon tube type 

Various 

Multiple 

Multiple 

8 

Differential Pressure Transducer 

Rosemount 

1151DP4E22B1 

11327273 

9 

Wattmeter, AC 

Ohio Semitronics 

GW5-020B 

49070 

10 

Wattmeter, DC 

Ohio Semitronics 

PC8-001-06-X5 

8043804 

11 

Rotameter 

Omega 

5-65 C 

None 

12 

Rotameter 

Fischer & Porter 

10A6131M 

8911A0675A17 

13 

Rotameter 

Fischer & Porter 

10A6131M 

8911A0675A17 

14 

Wattmeter, AC 

Ohio Semitronics 

GW5-020B 

53247 

15 

Wattmeter, AC with 

Ohio Semitronics / 

GW5-020B / 

53247/ 


external AC current transducer 

Omega 

CTT-050005 

None 

16 

Wattmeter, DC 

Ohio Semitronics 

PC8-001-06-X5 

8043804 

17 

AC/DC Current Shunt 

Empro Shunts 

C-1500 50 

None 

18 

AC/DC Current Shunt 

Empro Shunts 

C-1500 50 

None 

19 

Differential Pressure Transducer 

Rosemount 

1151DP5E22 

1072957 

20 

Data Acquisition / Switch Unit 

Hewlett Packard 

HP 34970A 

GT 0091257 


Table 3.3 (continued) 


Label 

Manufacturer Stated 
Range 

Manufacturer Stated Accuracy 

Calibration Source 

Calibration Date 

1 

1.5-21 gph, water 


IHHH 


2 

0.2 - 2.5 gpm, water 

0.1 gpm 

- 

- 

3 


0.1% repeatability 

Blasick, GT 

03/26/98 

4 

EEBElSsEjflHHB 

0.1% repeatability 

Dowling, GT 

10/12/99 

5 

0-1375 kPa 

0.25% cal. range 

Dowling, GT 

09/17/97 

6 

-200° - 425°C 

0.1 °C 

Dowling, GT 

11/04/96 

7 


2 psi 

- 

- 

8 


0.2% calibrated range 

Blasick, GT 

04/09/98 

9 

0 - 1 mA, 0 - 4 kW 

0.2% reading, 0.04% F.S. 

Manufacturer 

12/17/96 

10 

± 5 V, 0 - 25 kW 

1.0% reading F.S. 

Manufacturer 

01/21/98 

11 

0-9 gph, water 




12 

iRfiisrnnBi 

0.2 scfm, air 

- 

- 

13 


0.2 scfm, air 

- 

- 

14 

0-1 mA,0-4kW 

0.2% reading, 0.04% F.S. 

Manufacturer 

01/21/98 

15 

^ mm 

0.2% reading, 0.04% F.S. 

± 1.5% at 60 Hz with max burden 

Manufacturer 

01/21/98 

16 


1.0% reading F.S. 

Manufacturer 

01/21/98 

17 

0- 1500 A 

0.5% reading, with 50 mV drop 

Manufacturer 

June 1999 

18 

0- 1500 A 

0.5% reading, with 50 mV drop 

Manufacturer 

June 1999 

19 

0 - 750 inch H,0 

0.2% calibrated range 

Narrow, GT 

09/16/97 

20 

1 OOmV - 300V 
-150°- 1000°C 

0.0040% of reading 
± 1.0°C 

Manufacturer 

- 
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III.3.1.2 Differential Pressure Transducer 


Two differential pressure transducers (Components 8 and 19) were mounted on 
the GTMTF for the purpose of measuring the differential pressure across the test section. 
Both transducers were Rosemount style 1151DP transducers; however, Component 8 had 
a differential pressure range of 0 - 100 inch H 2 O while Component 19 had a differential 
pressure range of 0 - 750 inch H 2 O. Due to the high flow rates and associated large 
pressure drops required by this investigation, only Component 19 was utilized. 

The differential pressure transducer was mounted above the test section which 
caused entrained gas bubbles to become trapped in the transducer. For this reason, two 
toggle valves (Components Q) were installed and used as bleed valves to rid the 
differential pressure transducer of any trapped air pockets during the System Filling 
Procedure described in Section IV.2. 

111.3.1.3 Turbine Flowmeters 

Two EG&G turbine flowmeters (Components 3 and 4) were used to measure the 
flowrate through either the test section or the test section bypass line, depending on the 
flow loop configuration. The flowmeters were installed upstream of the test section in a 
horizontal orientation to ensure accurate readings. 

111.3.1.4 Rotameters 

Two rotameters (Components 1 and 2) were mounted in the GTMTF to allow 
visual verification of flowrate through the primary flow loop. In addition, two identical 
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low capacity rotameters (Components 12 and 13) were used to monitor and control the 
flow into and out of the saturation tank during experimentation. The saturation loop 
rotameters were used to prevent the saturation tank from either overflowing or drying out, 
while the principal function of the two primary loop rotameters was to provide a means 
of visually checking the flow for air pockets during the System Filling Procedure 
described in Section IV.2. Each rotameter was mounted vertically to ensure accurate 
readings and none of the four rotameters provided information essential to the analysis of 
experimental data. 

III.3.1.5 Thermocouples 

All thermocouples (Components 6) used in the GTMTF were Omega type-E with 
stainless steel sheaths. The thermocouples were mounted in the GTMTF using 
compression fittings to ensure the tips were positioned near the center of the water 
stream. During CHF experimentation (Setup 13), three fine gauge, type-E thermocouples 
(not shown) were mounted radially on the outer heater, 120° apart, one centimeter from 
the end of the test section heated length. In addition, one stainless steel sheathed, 
ungrounded, type-T thermocouple was mounted inside the inner heater, one centimeter 
from the end of the test section heated length. These thermocouples were used to detect 
CHF; a detailed explanation of their experimental arrangements can be found in Sections 
III.4 and III.5. 
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III.3.1.6 Pressure Gauges 

Four pressure gauges (Components 7) were used to provide visual verification of 
pressure levels in the system; these gauges were not used to provide the data necessary 
for analysis procedures. A pressure gauge manufactured by McDaniel Controls, Inc. and 
rated for 0 - 200 psig, was used to indicate the test section exit pressure. When filling the 
system during the System Filling Procedure outlined in Section IV.2, this gauge was 
monitored to ensure that the system did not pressurize above 20 psig. The peristaltic 
pump (Component AA) used to fill the system was not rated for high pressure and, thus, 
forcing the system above 20 psig could have damaged the pump. A second pressure 
gauge, manufactured by Span Instrumentation and rated for 0 - 300 psig, was used to 
monitor the pressure in the saturation tank. The reading on this pressure gauge, under 
normal operating conditions, always exactly matched the reading on the pressure gauge 
located at the test section exit. This fact allowed the operator to perform a rapid visual 
check to determine if an error had occurred in the system requiring termination of an 
experiment. 

The final two pressure gauges in the system were placed at the exits of the two 
centrifugal pumps used in Flow Loop Configuration Three. These pressure gauges, 
manufactured by U.S. Gauge and rated for 0 - 300 psig, were used to ensure that the 
centrifugal pumps were not overpowering the flow loop. The pressure gauge located 
after the second centrifugal pump was fitted with a pressure release valve. If the pressure 
at the centrifugal pump exit exceeded 300 psig due to an incorrect valving pattern in the 
system, the pressure release valve would trigger and prevent damage to sensitive 
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components in the GTMTF. 


III.3.1.7 Data Acquisition System (DAS) 

Data from all electronic sensors (Components 3 - 6, 8 - 10, and 14-19 of Figure 
3.1) were recorded using the Keithly MetraByte model DAS-8 data acquisition card with 
Keithly MetraByte model EXP-16 channel multiplexer/amplification boards. The DAS-8 
digitized all analog signals using a 10 kHz system clock and provided for software 
averaging of the digitized values to remove potential AC noise pick-up. The DAS-8 
system was run using an in house Microsoft Visual Basic 3.0 program on a 120-MHz 
Pentium computer that utilized Microsoft Windows 3.1 as an operating environment. 
This data acquisition program enabled the operator to control the time interval (1 s, 5 s, 
15 s, 1 m, 5 m, or 15 m) at which the data channels were sampled. This investigation 
only utilized the one and five second sampling time intervals. The DAS recorded all data 
for a particular experiment to a prescribed text file on the Pentium computer for future 
analysis. 

III.3.1.8: Hewlett Packard Data Acquisition / Switch Unit 

The Hewlett Packard Data Acquisition / Switch Unit (not shown in Figure 3.1) 
was used during CHF experimentation to increment power to the test section and to 
monitor three type-E thermocouples that were mounted on the outer heater. An input 
signal was provided to the unit from a Compaq DeskPro 575 personal computer (HP 
Benchlink Data Logger Software in a Microsoft Windows 95 operating environment) and 
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the HP unit, in turn, provided a signal to the power supply, which controlled the power 
input to the test section. The HP unit also monitored the outer heater thermocouples and 
automatically terminated the experiment by cutting the power supply signal if the outer 
heater surface temperature exceeded 300°C. The HP unit could sample at rates up to 600 
readings per second on a single channel and scan at rates up to 250 channels per second 
making it ideal for monitoring the outer heater temperature during a CHF experiment. 
All temperature readings from the outer heater thermocouples were recorded along with a 
real time stamp by the HP Benchlink software package on the Compaq personal 
computer for future analysis. 


III.3.2 Flow Loop Configuration One 

Flow Loop Configuration One was used for high flow / high heat flux OFI 
experimentation. As shown in Figure 3.2, one multi-stage centrifugal pump was utilized 
to force water through the high flow measurement line and through the test section. The 
FT-8 model high capacity flow transducer (Component 4 of Figure 3.2) was used to 
measure the test section volumetric flowrate which was controlled by the coarse and fine 
high flow line metering valves (Components K and L). Note that both the test section 
bypass line and the pump bypass line were closed creating a closed primary flow loop 
with no branching flow patterns. 


59 








Figure 3.2: Schematic Diagram of Flow Loop Configuration One 
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III.3.3 Flow Loop Configuration Two 


Flow Loop Configuration Two was used for low flow / low heat flux OFI 
experimentation. While testing new outer heater designs, several experiments (Setups 1 - 
3) were run, where only the outer heater was energized; no power was supplied to the 
inner heater. As a result, the working fluid was heated far less than necessary to make 
Flow Loop Configuration One effective for OFI experimentation. The FT-8 turbine 
flowmeter used by Flow Loop Configuration One to measure test section flowrate has a 
minimum flowrate of 0.2 gpm (45.4 L/hr). This value was far greater than the OFI points 
experienced in Setups 1-3 making alteration of the flow loop configuration necessary. 

As shown in Figure 3.3, one multi-stage centrifugal pump was utilized by Flow 
Loop Configuration Two to force water through both the low flow measurement line and 
the test section bypass line. This branched flow pattern significantly reduced the flowrate 
through the test section and allowed a great deal of flow control. The FTO-3 model low 
capacity flow transducer (Component 3 of Figure 3.3) was used to measure the test 
section volumetric flowrate, which was controlled by the low flow line metering valve 
(Component H). Likewise, the FT-8 model high capacity flow transducer was used to 
measure the test section bypass volumetric flowrate, which could be controlled by the 
coarse and fine high flow line metering valves (Components K and L). Note that the 
pump bypass line was closed while the test section bypass line was open, thus creating a 
closed primary flow loop with a flow pattern that branched immediately before the test 
section and converged just prior to entering the system heat exchanger (Component R). 
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III.3.4 Flow Loop Configuration Three 


Flow Loop Configuration Three was used for high flow / high heat flux OFI 
experimentation and was only a slight modification of Flow Loop Configuration One. As 
the OFI experimental power levels increased, OFI began to occur at or near the maximum 
flowrates achievable by Flow Loop Configuration One making modification necessary. 
A second centrifugal pump was added to the system in series with the original centrifugal 
pump. 

As shown in Figure 3.4, two multi-stage centrifugal pumps were utilized to force 
water through both the high flow measurement line and the pump bypass line. The FT-8 
model high capacity flow transducer (Component 4 of Figure 3.4) was used to measure 
the test section volumetric flowrate, which was controlled by the coarse and fine high 
flow line metering valves (Components K and L). Likewise, the flowrate through the 
pump bypass line was controlled using the pump bypass valve (Component E), but this 
flowrate was never measured or recorded. The pump bypass valve was only adjusted 
prior to the initiation of an experiment in order to either increase or decrease the 
maximum flowrate supplied to the test section. Note that the test section bypass line was 
closed while the pump bypass line was open, thus creating a closed primary flow loop 
with a flow pattern that branched immediately after the second centrifugal pump and 
converged in the degassing tank (Component X). 
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Figure 3.4: Schematic Diagram of Flow Loop Configuration Three 














III.3.5 Flow Loop Configuration Four 


Flow Loop Configuration Four was used for all CHF experiments. As shown in 
Figure 3.5, one variable speed positive displacement pump was utilized to force water 
through the low flow measurement line and through the test section. The FTO-3 model 
low capacity flow transducer (Component 3 of Figure 3.5) was used to measure the test 
section volumetric flowrate, which was controlled by the tachometer on the positive 
displacement pump. The low flow line metering valve (Component H) was partially 
closed for all CHF experiments to provide the positive displacement pump a small 
amount of back pressure (100 kPa). This back pressure helped reduce two-phase pressure 
and flow fluctuations encountered during CHF experimentation. Note that both the test 
section bypass line and the pump bypass line were closed creating a closed primary flow 
loop with no branching flow patterns. 
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Figure 3.5: Schematic Diagram of Flow Loop Configuration Four 
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III.3.6 Control of Flow Parameters 

Test section inlet temperature and exit pressure were controlled in the same 
manner for each of the four flow loop configurations described above. As explained 
earlier, modifications to the flow loop were required in order to achieve the necessary 
flow rate control for a specific set of flow parameters. This typically involved the 
addition of a pump to the system or slight modification of the test section flow path or 
flow rate instrumentation as described in Sections III.3.2 - III.3.5. Regardless of the test 
section flow path, however, the components used to pressurize the system were identical 
for all experimental setups except Setup 7, which is described in Section III.3.6.3 below. 

III.3.6.1 Test Section Inlet Temperature Preparation 

Test section inlet temperature was set and adjusted using three heat exchangers 
(Components 0, R, and DD of Figure 3.1) and/or two degassing tank cartridge style 
heaters (Components EE). The most significant temperature control system in the 
GTMTF was the system heat exchanger (Component R). The system heat exchanger was 
a high capacity, coiled, copper / copper, counter-flow, tube-in-tube heat exchanger that 
was cooled by the house water line (17° - 20°C) and located immediately after the test 
section exit. Due to its location, the system heat exchanger had to remove a large 
percentage of the heat in the working fluid since two-phase water exited the test section 
and had to be cooled prior to passing through the pump(s) for recirculation. A throttling 
valve (Component P) controlled the flow of house water through the secondary side of 
the system heat exchanger. This flow of house water was typically maximized through 
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the system heat exchanger to facilitate maximum cooling of the working fluid since 
manufacturer design specifications dictated that the inlet temperature to the pump(s) 
could not exceed 70°C. 

The test section heat exchanger (Component O) was located immediately before 
the test section. The test section heat exchanger was a coiled, stainless steel / copper, 
counter-flow, tube-in-tube heat exchanger, which was connected to a Fisher Scientific 
Isotemp Refrigerated Circulator (Model 9500) constant temperature bath. The set point 
of the constant temperature bath was selected manually by the operator and maintained 
automatically by the unit (± 0.2°C) during normal operation. The flow provided by the 
constant temperature bath through the secondary side of the test section heat exchanger 
was maximized at all times since the constant temperature bath could be set and 
maintained at a desired level. 

A low capacity, copper heat exchanger was located inside the degassing tank. 
This degassing tank heat exchanger (Component DD) was attached to the house water 
line and used to remove heat from the working fluid during high heat flux / high mass 
flux experiments where the system heat exchanger alone was not adequate. The flow 
through the secondary line of the degassing tank heat exchanger was controlled using a 
throttling valve (not shown). 

For approximately ninety percent of all experiments performed in this 
investigation, the three heat exchangers working in concert were adequate to maintain the 
test section inlet temperature. However, during high inlet temperature, low heat flux OFI 
experiments and high inlet temperature, low mass flux CHF experiments, the system heat 
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exchanger cooled the working fluid too much. As a result, the test section heat 
exchanger, which had a much lower capacity, could not maintain the required test section 
inlet temperature. To counter this effect, two variac controlled (0-6 kW), cartridge style 
degassing tank heaters (Components EE) were used to heat the working fluid prior to its 
introduction to the pump(s). These two heaters offered much finer control of system 
temperature than simple throttling of the secondary flow through the system heat 
exchanger. Note that in order to prevent component damage, care had to be used when 
operating these heaters. As mentioned earlier, the inlet temperature to the pump(s) could 
not rise above 70°C and thus, close observation of the degassing tank temperature and 
control of the degassing tank heaters was required in order to assure operation within 
such limits. 

Finally, a thermostatically controlled in-line heater (Component G) was located 
just before the test section inlet. This in-line heater was built in house and could deliver 
lkW of power to the working fluid; however, due to the control offered by the three heat 
exchangers and two degassing tank heaters mentioned above, the in-line heater was never 
used in any of the experiments performed in this investigation. 

III.3.6.2 Working Fluid Gas Saturation and System Pressurization 

All experiments in which data were acquired in this investigation utilized the 
same components to pressurize the system and saturate the working fluid with non¬ 
condensable gas. A compressed gas cylinder (Component JJ of Figure 3.1) with 
regulator was used to set the system pressure via the saturation tank (Component FF) by 
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means of flexible Swagelok tubing (1.70 MPa max) and the saturation tank gas fill valve 
(Component KK). The regulator (National Gas Cylinder model N1601 for air and The 
Matheson Company model H63 for nitrogen) was used to set the pressure inside the 
saturation tank. In turn, the open saturation tank pressure regulation valve (Component 
S) effectively equalized the test section exit pressure with the saturation tank pressure. 
During normal operation, there was no flow through this line. Gas saturation of the 
working fluid was accomplished by allowing the compressed gas to bubble through the 
saturation tank and escape at a controlled rate through the rotameter (Component 11) 
located atop the saturation tank liquid separation column (Component GG). The 
compressed gas bubbled through the saturation tank was either air or nitrogen, depending 
on the experimental setup. Note, as described in Section III.3.6.3 below, the non¬ 
condensable gas supplied to the system was switched from air to nitrogen after Setup 5 to 
avoid corrosion problems. Circulating the working fluid from the degassing tank, 
through the saturation tank, and back to the degassing tank ensured saturation of all water 
in the system. This path was commonly referred to as the saturation loop. The saturation 
tank recirculation pump (Component RR) circulated water through the saturation loop 
and the flow into and out of the saturation tank was monitored using two identical 
rotameters (Components 12 and 13). Equalization of the inlet and exit flow rates ensured 
that the saturation tank would neither overflow nor dry out; either situation resulting in 
termination of the experiment being performed. The exact saturation loop startup 
procedures are detailed in Section IV. 3. 









III.3.6.3 Minimization of Oxidation / Clarification of Setup 7 

Following higher power OFI experiments with Setup 5, a great deal of oxidation 
was observed on the surface of the inner heater. In an effort to minimize this type of 
damage, two new setups were tested to reduce the amount of dissolved oxygen present in 
the working fluid: Setups 6 and 7. For Setup 6, the gas saturation tank used to fix the test 
section exit pressure was pressurized using nitrogen rather than compressed air. No 
efforts were made to remove the ambient dissolved oxygen present in the distilled water 
due to contact with the atmosphere prior to its introduction to the system. In Setup 7, the 
exit pressure was set using an accumulator rather than the saturation tank. Again, no 
effort was made to remove the ambient dissolved oxygen present in the distilled water. 
Use of the accumulator pressurization method, however, was abandoned because 
maintenance of the test section exit pressure was impossible at high test section flow 
rates. Instead, all remaining experiments were pressurized using the saturation tank and 
compressed nitrogen. This pressurization method reduced the dissolved oxygen present 
in the working fluid while maintaining its non-condensable gas saturation status. Setup 7 
was never used to acquire data and, therefore, no further details will be provided about 
this setup. 


III.4 Inner Heater Configuration 

The “test section” referred to throughout this thesis was actually comprised of two 
separate pieces: an inner heater and an outer heater. The inner heater was mounted 
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concentrically within the outer heater to form an annular channel and the two heaters 
were typically powered by separate power supplies. This section describes all physical 
aspects of the inner heater and the power supply used to generate the required surface 
heat flux. As outlined in Table 3.1, a total of nine different inner heater configurations 
were used in this investigation. The inner heater configuration used in each experimental 
setup includes the following information. 

1) Inner Heater Label 

2) Power Supply Configuration 

3) Heated Length 

Table 3.4 categorizes the specific components used in these three categories for all nine 
inner heater configurations. Each time one or more of these three categories was altered, 
the inner heater configuration number was sequentially incremented. 


Table 3.4: Summary of Inner Heater Configurations 


Inner Heater Configuration 

Inner Heater Label 

Power Supply Configuration 

Heated Length (mm) 

1 

D 

None (unheated) 

N/A 

2 

E 

None (unheated) 

N/A 

3 

E, F, G 

2 

190 

4 

H 

2 


5 

I, J, L 

2 

180 

6 

K 

None (unheated) 

N/A 

7 

M 

3 

180 

8 

N 

3 

190 

9 

_o_ 

4 

190 


Note that the Inner Heater Label category begins with D. Since this investigation 
was intended to complement the work of Blasick (1999), the integrity of the labeling 
scheme she adopted was preserved whenever possible and she ended her work with Inner 
Heater Label C. Also, each individual inner heater used in an experiment was given a 
label (etched onto the copper electrode at the inlet end). When an effective inner heater 
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design was found, several duplicates were typically made which explains the multiple 
Inner Heater Labels shown in Table 3.4 for Inner Heater Configurations 3 and 5. 

In order to achieve the high power densities required by this investigation, a 
robust inner heater was necessary. All of the inner heaters for this research were 
manufactured in-house from standard, commercially available welded seam tubing and 
all were heated using direct resistive (Ohmic) heating. The electrodes used for the 
heaters were made from solid, Vi” diameter, free-machining copper rods (alloy Cl4500). 
Free-machining copper is an electrical-grade alloy notable for its machining properties 
and an electrical conductivity essentially equal to that of pure copper. The copper rods 
were cut to length and, using a lathe, the ends were turned down to a diameter equal to 
the inner diameter of the tubing selected for the heated portion (Inconel-600 or Hastelloy 
C-276). The copper electrodes were then press-fit into either end of the heater tubing. 
Next, a small hole was drilled and tapped through the tubing and into, but not through, 
the electrode. A brass screw (3-48 UNC) was then tightened into this threaded hole and 
sheared off even with the outer wall of the tubing. These brass stakes ensured that a 
tight, corrosion free conduction path was always present from the copper electrodes to the 
heater tubing. In general, this pinning method worked well, so that no oxidation of the 
copper or of the tubing in the region of contact was ever observed when selected heaters 
were disassembled after use. This is a remarkable fact given the high electrical currents 
(*4000 Amps) carried by the inner heater. Thus, it was found that the exact lengths of 
the electrodes and the precise amount of overlap between the tube wall and electrodes 
were not critical parameters. Extreme care was taken, however, to manufacture the inner 


73 







heaters in a symmetric way, so that measurements made relative to the exposed electrode 
ends could be used to verify that the heater was mounted with its heated region in the 
desired location. For all inner heaters, the heated region was considered to begin and end 
at the tips of the press-fit copper electrodes and the exposed ends of the copper electrodes 
were used to connect power supply leads to the inner heater. 

The following subsections describe the design specifications unique to each of the 
twelve different inner heaters (D - O) used in this investigation. These subsections will 
include a discussion on the heated length of each inner heater. The Power Supply 
Configuration category will be discussed in Section III.6 since three of the four different 
power supply configurations were used to heat both the inner and outer heaters. 


III.4.1 Inner Heater D 

Inner Heater D (shown in Figure 3.6) was inherited from Blasick (1999)and used 
in Setup 1 to test the feasibility of a proposed outer heater design. It was merely used to 
maintain an annular cross-section flow channel through the test section and, thus, was 
never heated. The relatively short length of Inner Heater D made construction of a new 
inner heater imperative before power could be applied to the inner heater. Inner Heater D 
was too short and thus left an inadequate amount of electrode length protruding from the 
ends of the outer heater on which to attach power supply leads. 
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Figure 3.6: Diagram of Inner Heater D 

III.4.2 Inner Heaters E. F. and G 

Inner Heaters E, F, and G (shown in Figure 3.7) improved upon the design of 
Inner Heater D, not only by extending the overall length of the heater, but also by adding 
a rigid alumina rod to the construction. The alumina rod provided stability to the inner 
heater and prevented it from bowing. Fortunately, such cylindrical rods were readily 
available in the exact diameter needed to provide a tight fit with the nominal !4” (6.35 
mm) diameter heater tubing with a 0.028” (0.7112 mm) wall thickness. Figure 3.7 shows 
the design of these heaters with the actual diameter of the Inconel-600 tubing (0.252” / 
6.40 mm). 

Heat fluxes above 1.5 MW/m 2 tended to cause a significant amount of oxidation 
at discrete points on the surface of this inner heater design after repeated use. It was 
unclear whether this damage was caused by electrical arcing (to the grounded inner 
surface of the AC outer heaters) or by extreme heating caused by highly localized film 
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boiling. In spite of the visible damage, these heaters remained very straight through all 
experiments as a result of the alumina rod insert. Typically, failure of this heater design 
occurred when the localized surface damage was severe enough to damage or break the 
alumina rod. Catastrophic failure of the heater (i.e., bowing and melting) could then 
occur. 


* ALL DIMENSIONS IN MILLIMETERS 



Figure 3.7: Diagram of Inner Heaters E, F, and G 


III.4.3 Inner Heater H 

Inner Heater H was produced from 14” diameter Inconel-600 tubing that had a 
wall thickness of 0.033” (0.8382 mm). This larger wall thickness made insertion of the 
alumina support rod impossible, which weakened the structural stability of the heater. In 
addition, the heated length of Inner Heater H differed from that of Inner Heaters E, F, and 
G. Due to limitations imposed by the outer heater mounting mechanism, the outer heater 
heated length was limited to 180 mm at that point in the experimental process. It was 
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determined that deviation from the nominal heated length of 190 mm was less important 
than matching the inner and outer heater heated lengths. Inner Heater H bowed and 
failed rather quickly due to its lack of rigidity and was never duplicated. 


* ALL DIMENSIONS IN MILLIMETERS 



Figure 3.8: Diagram of Inner Heater H 


III.4.4 Inner Heaters I. J. K. and L 

Inner Heaters I, J, K, and L preserved the 180 mm heated length used in Inner 
Heater H, but utilized V” diameter Inconel-600 tubing with a 0.028” (0.7112 mm) wall 
thickness. This change in heater tubing allowed the addition of an allumina support rod 
as outlined in the description of Inner Heaters E, F, and G in Section III.4.2. Shown in 
Figure 3.9, this inner heater design remained extremely straight in all experiments for 
which it was used, however, it was plagued by corrosion problems. 







* ALL DIMENSIONS IN MILLIMETERS 



Figure 3.9: Diagram of Inner Heaters I, J, K, and L 

III.4.5 Inner Heater M 

As shown in Figure 3.10, Inner Heater M matched the design of Inner Heater H 
exactly, except that Hastelloy C-276 was used for the heated section rather than Inconel- 
600. This material change was adopted in an attempt to alleviate the localized corrosion 
encountered with the Inconel-600. Hastelloy C-276 has better corrosion properties than 
Inconel-600, but only tubing with a wall thickness of 0.033” (0.8382 mm) was available, 
preventing the use of an alumina support rod. Inner Heater M showed less corrosion than 
the inner heaters constructed of Inconel-600, but it bowed unacceptably at high heat 
fluxes thus proving the necessity of the alumina support rod. 
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* ALL DIMENSIONS IN MILLIMETERS 



Figure 3.10: Diagram of Inner Heater M 


III.4.6 Inner Heaters N and O 

Inner Heaters N and O were adapted from the design of Inner Heaters E, F, and G 
for use primarily in CHF experimentation. As shown in Figure 3.11, Inner Heaters N and 
O provided the original nominal heated length of 190 mm and they included an alumina 
support rod. Also, the inlet end electrode was constructed from a 14” diameter solid 
copper rod as with all previously mentioned inner heater designs. 

In order to detect the approach of CHF, Inner Heaters N and O were fitted with an 
internally mounted, type-T, ungrounded thermocouple. A single thermocouple was 
permanently mounted in a short, 14” diameter, nylon sleeve using standard two-part 
epoxy. The nylon sleeve allowed mounting of the thermocouple to the inner heater while 
the epoxy provided a permanent, watertight seal around the thermocouple. The heater’s 
exit end electrode was fabricated from a 14” diameter, brass tube whose end was turned 
down and press-fit, but not staked, into the Inconel-600 heater tube. A Swagelok union 
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fitting was then attached to the end of the brass tube to allow mounting of the 
thermocouple. Next, the tip of the thermocouple was slightly bent and inserted an 
appropriate distance into the inner heater. The bend in the tip of the thermocouple 
ensured contact with the inner heater wall and the thermocouple was inserted such that it 
would contact the inner heater wall one centimeter before the end of its heated length. 
Once in position, the thermocouple was secured by tightening the Swagelok union around 
the nylon sleeve. The presence of the thermocouple in this inner heater design forced the 
alumina support rod to be shortened by one centimeter during construction. Unlike 
previous inner heater designs, the alumina support rod was in contact with the inlet end 
electrode, but not the exit end electrode. This design modification did not significantly 
weaken the inner heater and no detectable bowing of Inner Heaters N or O was ever 
observed. 

Note that this thermocouple mounting style allowed the same thermocouple to be 
easily transferred from one inner heater to another, thus eliminating data inconsistency 
and waste of material. In addition, the brass tube, which acted as the exit end power 
electrode, allowed insertion of the thermocouple, but it increased the impedance of the 
electrical connection. The increased resistivity of this electrode and the associated heat 
produced outside the test section area was countered by active cooling of the power 
supply connection as detailed in Section III.6. Finally, note that both the nylon sleeve 
and the epoxy were non-conducting materials, thus ensuring electrical independence 
between the inner heater and the thermocouple. 
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Figure 3.11: Diagram of Inner Heaters N and O 


III.5 Outer Heater Configuration 

The outer heater configuration describes all physical aspects of the outer heater 
and the power supply used to generate the required surface heat flux. As outlined in 
Table 3.1, a total of eight different outer heater configurations were used in this 
investigation. The outer heater configuration used in each experimental setup includes 
the following information. 

1) Outer Heater Label 

2) Power Supply Configuration 

3) Heated Length 

4) Wire Wrapping Style 

Table 3.5 categorizes the specific components used in these four categories for all eight 
outer heater configurations. Each time one or more of these four categories was altered, 
the outer heater configuration number was sequentially incremented. 


81 








Table 3.5: Summary of Outer Heater Configurations 


Outer Heater 
Configuration 

Outer Heater Label 

Power Supply 
Configuration 

Heated Length 
(mm) 

Wire Wrapping Style 

1 

H 

1 

55 

1 

2 

H 

1 

53 

1 

3 

H 

1 

123 

2 

4 

H 

1 

180 

3 

5 

E 

1 


3 

6 

L 

2 

185 

N/A (tube heater) 

7 

I,J 

3 

185 

N/A (tube heater) 

8 

K 

4 

185 

N/A (tube heater) 


Note that the Outer Heater Label category begins with E. Since this investigation 
was intended to complement the work of Blasick (1999), the integrity of the labeling 
scheme she adopted was preserved whenever possible and she ended her work with Outer 
Heater Label D. Outer Heaters F and G were fabricated but never used in this 
investigation, therefore, the specifications of these two outer heaters will not be reported. 
Also, each individual outer heater used in an experiment was given a label (etched or 
stamped on the inlet end). In later experiments, when an effective outer heater design 
was found and tube heaters were expendable, duplicates were typically made which 
explains the multiple Outer Heater Labels shown in Table 3.5 for Outer Heater 
Configuration 7. 

The following subsections will describe the design specifications of the six 
different outer heaters (E, H - L) used in this investigation. These subsections will 
include a discussion on the heated length of each outer heater. In addition, a subsection 
will be dedicated to the discussion of the various Wire Wrapping Styles utilized on the 
AC powered outer heaters. The Power Supply Configuration category will be discussed 
in Section III.6 since three of the four different power supply configurations were used on 
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both the inner heater and the outer heater. 


III.5.1 Outer Heater E 

Outer Heater E (shown in Figure 3.12) was constructed entirely of free-machining 
copper (alloy C14500). First, a piece of 2.5” (63.5 mm) diameter copper bar stock was 
cut to length (11.5” / 292.1 mm) and sent to Precision Group Inc. / Deep Hole Drilling 
Division to have a circular channel drilled through it. The circular channel was cut to a 
tolerance of ±0.002” (0.0508 mm) and hone finished to a surface roughness of 0.000032” 
(0.0008128 mm). This channel served as the outer wall of the annular test section once 
the inner heater was in place. Next, using the inlet and exit of the circular channel as 
center, the outer heater was placed in a lathe and the ends were turned down such that 
they would fit tightly into a %” Swagelok fitting. Also, the 2.5” diameter unfinished 
outer surface was turned down a few millimeters to ensure that the remaining heater 
material had a smooth surface finish and was concentrically centered on the drilled 
circular channel. These actions allowed application of an electrical insulation layer and 
ensured even heating once the heater wire was wrapped around the outer heater. Next, 
the inlet and exit ends of the outer heater circular channel were widened in order to 
ensure that the nominal APT annular channel length of 250 mm was maintained. Note 
that the Wire Wrapping Style and the placement of the power supply electrodes 
determined the heated length on Outer Heater E. 

After Outer Heater E was constructed it was sent to Cincinnati Thermal Spray 
who coated the 8.5” (215.9 mm) larger diameter center section of the outer heater with an 
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electrical insulation layer. The insulating layer consisted of a thin nickel-aluminum 
bonding layer followed by a 0.007 — 0.010 inch (0.1778 - 0.2540 mm) aluminum oxide 
layer finished with a thin epoxy scaler. This insulating layer served to electrically isolate 
the copper structure of Outer Heater E from the heater wire used to heat it. Once the 
outer heater was constructed and coated, it was wrapped with wire as described in 
Section III.5.3 and placed in the GTMTF for experimentation. 



* ALL DIMENSIONS IN MILLIMETERS 

Figure 3.12: Diagram of Outer Heater E 

III.5.2 Outer Heater H 

Outer Heater H was constructed and heated in an identical manner to Outer 
Heater E, however, Outer Heater H provided a larger channel diameter which resulted in 
a larger annular gap width. Figure 3.13 shows a diagram of Outer Heater H. 
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* ALL DIMENSIONS IN MILLIMETERS 

Figure 3.13: Diagram of Outer Heater H 


III.5.3 Wire Wrapping Style 

Wire Wrapping Style refers to the type of wire and the method of wrapping used 
for Outer Heaters E and H. Table 3.6 outlines the three different wire wrapping styles 
used on Outer Heaters E and H. 


Table 3.6: Summary of Wire Wrapping Styles 


Wire Wrapping Style 

Type of Wire 

Turns per Inch 

Number of Electrodes 

Contacts per Electrode 

1 

NiChrome 80, 
26 AWG 

22 

2 

2 

2 

NiChrome 80, 
26 AWG 

22 

3 

3 

3 

NiChrome 80, 
26 AWG 

22 

4 

3 


All three wrapping styles utilized the same type of wire and turns per inch on the 
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outer heater. The bare NiChrome 80, 26 AWG wire had a resistance of 2.647 Ohms per 
foot (0.00868 Ohms per millimeter) and was wrapped at 22 turns per inch (0.8661 turns 
per millimeter). A lathe operating at low speed was utilized to wrap the outer heater each 
time in order to ensure precise duplication of the turns per inch. A small plastic sleeve 
was slipped over the inlet end of the outer heater, which was in turn mounted in the chuck 
of the lathe. The lathe’s tale stock was then fitted with a live center which, in turn, was 
tightened into the exit hole of the outer heater. Next, a separate plastic piece was 
mounted to the lathe’s traversing mechanism to act as a guide for the wire. The end of 
the wire was then fed through the plastic guide and secured to the plastic sleeve. At this 
point, the appropriate gears were engaged on the lathe to produce the slowest outer heater 
turning speed and a traversing feed rate of 22 turns per inch. An operator then manually 
placed tension on the wire by holding the spool as tightly as possible while a second 
operator energized the lathe. The outer heater was allowed to wrap until the end of the 
aluminum oxide coated region was reached. At that point, the lathe operator disengaged 
the lathe while the wire holder maintained constant tension on the wire. The lathe 
operator then secured the wire to the outer heater using the specially designed steel 
Heater Wire Clamp shown in Figure 3.14. Once the wire was secured to the outer heater, 
the wire holder released the tension and the lathe operator cut away the excess wire. A 
hand-held Fluke 25 Multimeter was then used to check for electrical isolation between 
the wrapped NiChrome 80 wire and the copper portion of the outer heater. At that point, 
the outer heater was removed from the lathe and mounted in the GTMTF. 
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10.414 



* ALL DIMENSIONS IN MILLIMETERS 



Figure 3.14: Diagram of Heater Wire Clamp 


Once mounted, the wire on the outer heater had to be coated with a thin layer of 
Samereisen Electrical Resistor Cement (No. 78 paste). Normally, unsecured metal wire 
tends to lengthen as it is heated. When wrapped around the cylindrical outer heater, 
however, the loops of heated, bare wire became loose and came into contact with their 
neighboring loops. This contact, in turn, caused a short in the heater wire circuit. The 
short circuit caused uneven heating over the outer heater and typically lead to melting of 
the wire and damage to the aluminum oxide coating on the outer heater. The resistor 
cement was used to hold the wire securely to the outer heater, which ensured good 
thermal contact with the outer heater and prevented a contact short between two wire 


87 






loops. First, the small areas needed for the electrode contacts were masked off using 
tape. Next, the resistor cement was applied uniformly in an extremely thin layer over the 
entire outer heater. Care was taken to ensure that the resistor cement penetrated all gaps 
between the wire loops and that the coating layer did not become too thick. An ideal 
resistor cement layer just barely covered the heater wire such that, while wet, the heater 
wire could still be seen through the resistor cement. Trials showed that a thick coating of 
resistor cement tended to crack and expose highly localized pieces of heater wire to 
separation from the outer heater and melting. Once the heater wire was completely 
covered with resistor cement, the masking tape was removed and the cement was allowed 
to cure for twenty-four hours. After curing, the electrodes were put into place and the 
outer heater was ready for experimentation. 

The Number of Electrodes category outlined in Table 3.6 refers to the number of 
copper electrodes used to heat the outer heater for a given set of experiments. As many 
as four, 14” (6.35 mm) thick, nickel-plated copper electrodes (shown in Figure 3.15) were 
used to attain the required heated lengths for the OFI experiments that utilized AC 
powered outer heaters. Current passed into the nickel-plated copper ring and on to the 
outer heater wire through the nickel-plated brass screws that held the electrode in place. 
Knurled nuts were used to lock the brass screws in place as shown in Figure 3.15. 
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0 108.00 (APPROX CLEARANCE) 



* ALL DIMENSIONS IN MILLIMETERS 

Figure 3.15: Diagram of Outer Heater Electrode Ring 

Due to limitations in Power Supply Configuration One, the heater wire could not 
be heated as a single, continuous resistor. Instead, multiple electrodes were used to 
connect different heated zones of equal resistance in parallel on the outer heater. The 
tests run using Outer Heater Configuration One utilized a single, 55 millimeter heated 
zone (58.5 Ohms) with two electrodes. Outer Heater Configuration Two also utilized a 
single heated zone with two electrodes, except the heated length was only 53 millimeters 
(57.1 Ohms). Outer Heater Configuration Three had two heated zones of 60.0 and 60.2 
Ohms that combined electrically in parallel using three electrodes to form a single heated 
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length of 123 millimeters. Finally, Outer Heater Configuration Four utilized three heated 
zones of 55.8, 55.1, and 56.5 Ohms that combined electrically in parallel using four 
electrodes to form a single heated length of 180 millimeters. Electrical resistance 
readings were taken prior to each set of experiments performed and used to calculate and 
set the voltage required to generate the necessary heat flux on the outer heater 
( V = JPower * R eq ). Note that, due to the parallel nature of the electrical circuit formed, 

the three heated zones were exposed to the same voltage drop, but different amounts of 
current flowed through each zone due to the slightly differing resistances. These 
resistance readings fluctuated from one set of experiments to another since the test 
section was damaged and re-wrapped multiple times, however, the percent difference 
between the three resistance readings, i.e. the average heat flux over the different zones, 
never exceeded five percent. 

The Number of Contacts per Electrode simply refers to the number of nickel- 
plated brass screws that were actually in contact with the heater wire for each electrode. 
Due to the nature of the power configuration, one brass screw typically carried all or most 
of the current. For this reason, ideally, only one screw had to contact the heater wire 
while the other two could rest on the resistor cement used to electrically insulate the 
heater wire. Experimentation proved, however, that small imperfections in either the tip 
of the brass screw or the surface of the heater wire lead to rapid failure of a single contact 
point since there were no alternate paths through which the electricity could pass. Two 
contact points also proved to be less robust than required, so after Wire Wrapping Style 
One, all future wrapping styles utilized all three brass screws as electrical contacts. 
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III.5.4 Outer Heater I and J 


After repeated failures of Outer Heaters E and H, a tubular-style outer heater was 
adopted for use with a DC power supply. Figure 3.16 shows a schematic of Outer 
Heaters I and J. A single 3/8” (9.525 mm) diameter Hastelloy C-276 welded seam tube 
was cut in four places using an EDM machine to allow for the inlet and exit of cooling 
water. In addition, two sets of lines, one marking a 190.0 mm heated length and the other 
a 180.0 mm heated length, were lightly scribed into the surface of the heater using a 
lathe. These scribe marks allowed for simple electrode placement and connection at the 
desired heated length. 

Hastelloy C-276 was found to be the ideal material for the outer heaters due to its 
high strength at elevated temperatures. The single limitation of Outer Heaters I and J was 
the available inner diameter of the tube. Only an inner diameter of 0.304” (7.7216 mm) 
was obtainable, which created an annular gap width of only 0.660 mm. In addition, 
mounting constraints caused by the size of the outer heater electrodes (described in 
Section III.6) limited the maximum heated length to only 185.0 mm. The inlet end 
electrode was mounted on the 190.0 mm scribe line while the exit end electrode was 
mounted on the 180.0 mm scribe line, thus creating a 185.0 mm overall outer heater 
heated length. 
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* ALL DIMENSIONS IN MILLIMETERS 


Figure 3.16: Diagram of Outer Heaters I and J 

III.5.5 Outer Heater K 

Outer Heater K was constructed of Hastelloy C-276 and used to perform all of the 
CHF experiments conducted in this investigation. As shown in Figure 3.17, the 
dimensions of Outer Heater K were identical to those of Outer Heaters I and J, however, 
three type-E thermocouples were mounted one centimeter from the exit end 180.0 mm 
scribe line and used to detect the heater wall temperature spike associated with CHF. 
These thermocouples were mounted to Outer Heater K using a 'A” diameter clamp with 
ceramic pieces as 120° spacers. In addition to properly spacing the thermocouples, the 
ceramic pieces also electrically isolated and thermally protected the thermocouple leads. 
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(3) TO HP UNIT 

* ALL DIMENSIONS IN MILLIMETERS 



Figure 3.17: Diagram of Outer Heater K 


III.5.6 Outer Heater L 

Outer Heater L was constructed of Inconel-600 and used to test the feasibility of 







Inconel as an outer heater material. As shown in Figure 3.18, Outer Heater L had the 
same dimensions as Outer Heaters I and J. Only one experiment was attempted using this 
outer heater design. An OFI experiment with an unheated inner heater and 2.0 MW/m 2 
on Outer Heater L caused the outer heater to bend and glow. The experiment was 
terminated and this outer heater design was never duplicated. 



Figure 3.18: Diagram of Outer Heater L 


III.6 Power Sui 


Configuration 


As mentioned earlier, all test sections in this investigation were comprised of an 
inner heater and an outer heater, which together formed a heated annular flow passage. 
The heat flux for each heater could be set and maintained independently. Tables 3.4 and 
3.5 outline the Power Supply Configuration utilized by each heater for all experimental 
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setups. The following subsections describe the four power supplies used in this 
investigation and detail the power measurement points. Each circled numerical reference 
in the following power supply figures indicates a power measurement device detailed in 
Table 3.3. 


III.6.1 Power Supply Configuration One 

Power Supply Configuration One was only used to heat Outer Heaters E and H 
during OFI experimentation. As shown in Figure 3.19, AC power was measured before 
the transformer and thus was subject to electrical losses of three percent before reaching 
the outer heater load. Note that experimental measurements taken at the test section 
using independent voltage and current readings showed that the transformer caused 
electrical losses of three to five percent. A three percent power loss fraction was used as 
a conservative estimate for all analysis procedures as described in Section IV.7. 

Power Supply Configuration One controlled the in-house 60 Hz, 240 Volt AC 
electrical power input (in percentage) using two coupled Technipower type W 50 Variacs 
with 100 % representing a 240 Volt ouput from the Variacs. The power output from the 
Variacs was then measured by an Ohio Semitronics AC wattmeter with an external 
Omega Engineering AC current transducer (Component 15 of Figure 3.1). Next, the 
voltage was stepped up (2:1) using a General Electric transformer (catalog number 
9T21B9109) and supplied to two stationary nodes located near the outer heater using 
insulated, 6 AWG, stranded copper cables. Power was then supplied directly to the outer 
heater via two 1/16” (1.5875 mm) diameter, insulated copper wires connected to each of 
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the Outer Heater Electrode Rings described in Section III.5.3. Note that the copper 
portion of the outer heater was rigidly connected to the GTMTF via stainless steel 
Swagelok fittings and hence was electrically grounded. The wound heater wire was 
electrically isolated from the copper portion of the heater by the aluminum oxide layer 
described earlier; this was the portion of the outer heater to which power was supplied. 


VARIAC 

(DOUBLE GANG UNIT) 


WATTMETER 

HOUSING 


VOLTAGE 

TRANSFORMER 



OUTER HEATER 
LOAD 


Figure 3.19: Electrical Diagram of Power Supply Configuration One 


III.6.2 Power Supply Configuration Two 

Power Supply Configuration Two was used independently by both the inner 
heater and the outer heater for different experimental setups. Referring to Tables 3.1 and 
3.4, Inner Heater Configurations 3-5 used Power Supply Configuration Two while the 
outer heater was heated by Power Supply Configuration One. Likewise, referring to 
Tables 3.1 and 3.5, Outer Heater Configuration 6 used Power Supply Configuration Two 
to test Outer Heater L while the inner heater was left unheated. 

As shown in Figure 3.20, Power Supply Configuration Two used a DC power 
supply produced by Rapid Power Technologies (serial number 1198224) to supply 
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electrical power to the connected heater. The power level produced by this DC power 
supply was controlled using a set of two rheostats. The first rheostat controlled the 
maximum current produced by the power supply (in percentage) with 100 % representing 
a potential maximum current of 2500 Amps. The second rheostat controlled the voltage 
drop produced by the power supply (in percentage) with 100 % representing a maximum 
electrical potential of 25 Volts. Note that these two rheostats are not shown in Figure 
3.20. Power Supply Configuration Two measured the power output using an Ohio 
Semitronics DC wattmeter (Component 16 of Figure 3.1). When used to power the outer 
heater, Power Supply Configuration Two supplied power directly to the heater via one 
fifteen foot long, flexible, 1.25” (31.75 mm) diameter, 1000 MCM, insulated, copper 
cable connected to each of the two Outer Heater Electrode Plates (shown in Figure 3.21) 
that were mounted to the heater. When used to power the inner heater, electrical power 
was sent to a stationary positive node via a 1000 MCM, insulated copper cable. From 
there, power was applied to the inner heater using two flexible, insulated, AWG 4/0 
copper welding cables that were attached to the inner heater via a one inch wide copper 
electrical connection block. After passing through the inner heater, the inner heater 
power loop was completed by the negative side, which was a mirror image of the 
positive. 
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DC PDWER 
SUPPLY 



Figure 3.20: Electrical Diagram of Power Supply Configuration Two 



NOTES: 

1. Extremely smooth (honed) finish required 
on bored hole. 

2. All dimensions in millimeters unless 
otherwise specified. 


Figure 3.21: Diagram of Outer Heater DC Electrode Plate 
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III.6.3 Power Supply Configuration Three 

Power Supply Configuration Three was used simultaneously by both the inner 
heater and the outer heater during the last five OFI experiments performed in this 
investigation. In order to increase the surface heat flux on the outer heater above 2.0 
MW/m 2 , a DC powered tube heater design had to be adopted. This tube heater (Outer 
Heaters I - L), then had to be connected electrically in parallel with the inner heater to the 
Rapid Power Technologies (RPT) DC power supply described in Section III.6.2. 

Since Power Supply Configuration Three used the RPT DC power supply to 
supply electrical power to both the inner and the outer heater simultaneously, a method of 
connecting the heaters electrically in parallel had to be devised such that the inner-to- 
outer surface heat flux ratio remained unity. Due to the larger heater wall thickness of 
the outer heater tube, the outer heater leg of the power supply loop had a lower electrical 
resistance. Note that the outer heater, due to its larger diameter and accompanying larger 
heated area, required more power to attain the required surface heat flux. However, due 
to the nature of Power Supply Configuration Three, the outer-to-inner heater resistance 
ratio produced a much higher outer surface heat flux than that for the inner heater. This 
problem was corrected by the addition of a series resistor to the outer heater leg of the 
power supply loop as shown in Figure 3.22. The following paragraphs describe the 
components of the inner and outer heater legs of the power supply loop independently. In 
addition, Figure 3.23, though not drawn to scale, provides a general, less detailed 
overview of the order in which the components of Power Supply Configuration Three 


99 






were connected. 


DC POWER 
SUPPLY 



Figure 3.22: Electrical Diagram of Power Supply Configuration Three 



Figure 3.23: Component Diagram of Power Supply Configuration Three 
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The inner heater leg of Power Supply Configuration Three began at the positive 
node of the DC power supply and ran through an Empro DC current shunt (Component 
18 of Figure 3.1) and on to a flexible, 1000 MCM, copper cable. Next, power was 
applied directly to the inner heater using an Inner Heater DC Electrode Plate (Figure 
3.24). The 1000 MCM copper cable bolted directly to this electrode plate using a single 
!4”-13 bolt; the hole for this bolt was drilled and tapped in the center of the plate (hole 
not shown in Figure 3.24). After passing through the inner heater, the inner heater leg of 
the power supply loop was completed by a second Inner Heater DC Electrode Plate that 
was connected to the negative node of the DC power supply via a second 1000 MCM 
copper cable. Note that the voltage drop across the inner heater was taken using taps that 
were connected to the electrode plates. The DAS system multiplied this voltage drop by 
the current reading taken by the Empro DC current shunt to determine the power supplied 
to the inner heater. In addition, using this power level and an operator defined inner 
heater heated surface area, the DAS system continuously calculated and displayed the 
surface heat flux on the inner heater, which allowed proper adjustment of the power level 
prior to the commencement of an OFI experiment. 

The outer heater leg of Power Supply Configuration Three began at the positive 
node of the DC power supply and ran through an Empro DC current shunt (Component 
17 of Figure 3.1) and on to a fifteen foot long, flexible, 1.5” (38.1 mm) diameter, 1500 
MCM, insulated, copper cable. Next, a part of the supplied electrical power was 
dissipated in a variable series resistor. This series resistor was constructed from a !4” 
(6.35 mm) diameter brass tube that was cooled by a water bath and a house water flow- 
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through line. Power from the 1500 MCM copper cable was applied to the brass heater 
using an Inner Heater DC Electrode Plate (Figure 3.24). Moving the electrical 
connection point in order to make the resistor length longer or shorter altered the 
resistance of the series resistor as required. Next, power was applied to the outer heater 
using a flat, tin coated, copper mesh that was connected to an Outer Heater DC Electrode 
Plate (Figure 3.21). The copper mesh connected to both the brass tube electrode plate 
and the outer heater electrode plate using eight 10-32 screws that were fitted with 
washers, forced through the mesh, and threaded into the holes shown in Figures 3.21 and 
3.24. After electricity passed through the outer heater, the outer heater leg of the power 
supply loop was completed by a second Outer Heater DC Electrode Plate that was 
connected to the negative node of the DC power supply via a second 1500 MCM copper 
cable. The 1500 MCM copper cable bolted directly to the electrode plate using a single 
'/2”-13 bolt; the hole for this bolt was drilled and tapped in the center of the plate (hole 
not shown in Figure 3.24). Note that the voltage drop across both the outer heater and the 
series resitor was taken using taps that were connected to the electrode plates. The DAS 
system multiplied these voltage drops by the current reading taken by the Empro DC 
current shunt to determine the power supplied to the outer heater and the power 
dissipated by the series resistor. In addition, using the outer heater power level and an 
operator defined outer heater heated surface area, the DAS system continuously 
calculated and displayed the surface heat flux on the outer heater, which allowed proper 
adjustment of the power level prior to the commencement of an OFI experiment. 

Finally, a measurement was made of the total power level supplied by the DC 
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power supply using an Ohio Semitronics DC wattmeter (Component 16 of Figure 3.1). 
In both OFI and CHF experimentation this redundant reading became invaluable since 
the two DC current shunts proved unreliable. Due to an unidentified ground in the DAS 
system, the readings from both DC current shunts became too erratic for experimental 
analysis. Therefore, the power applied to the test section was determined using an energy 
balance and the total power reading from the DC power supply. This procedure is 
described in Sections IV.7.2 and IV.8.2. 



NOTES: 

1. Extremely smooth (honed) finish required 
on bored hole. 

2. All dimensions in millimeters unless 
otherwise specified. 


Figure 3.24: Diagram of Inner Heater DC Electrode Plate 
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III.6.4 Power Supply Configuration Four 

Power Supply Configuration Four (shown in Figure 3.25) was used for all 
CHF experiments performed in this investigation. The power supply loop was identical 
to that of Power Supply Configuration Three, however, the power supplied by the RPT 
DC power supply was not controlled using the two rheostats described in Section III.6.2. 
Instead, power was controlled using a Hewlett Packard Data Acquisition / Switch Unit 
(Component 20 in Table 3.3) that was connected to a Compaq DeskPro 575 personal 
computer (HP Benchlink Data Logger Software in a Microsoft Windows 95 operating 
environment). An input signal was provided to the HP unit from the Compaq personal 
computer and the HP unit, in turn, provided a signal to the power supply, which 
controlled the power input to the test section. The HP unit also monitored the outer 
heater thermocouples mounted on Outer Heater K (Figure 3.17) and automatically 
terminated the experiment by cutting the power supply signal if the outer heater surface 
temperature exceeded 300°C. The HP unit could sample at rates up to 600 readings per 
second on a single channel and scan at rates up to 250 channels per second making it 
ideal for monitoring the outer heater temperature during a CHF experiment. All 
temperature readings from the outer heater thermocouples were recorded along with a 
real time stamp by the HP Benchlink software package on the Compaq personal 
computer for future analysis. Figure 3.26, though not drawn to scale, provides a general, 
less detailed overview of the order in which the components of Power Supply 
Configuration Four were connected. 

During the first CHF experiment, all four DC electrode plates became extremely 
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hot due to contact resistance between the plates and their respective heaters. For this 
reason, an external, copper 54” (6.35 mm) diameter coolant tube was silver soldered to 
each of the four electrode plates attached to the test section. Each coolant tube was 10 
inches (254 mm) long and bent into a U-shape such that it would fit around the central 
row of 10-32 tapped holes with each leg of the U-shape running perpendicular to the 14”- 
20 electrode connection holes. The four coolant tubes were then attached in series using 
flexible, !4” (6.35 mm) diameter nylon tubing. Next, the house water flow-through line 
that cooled the adjustable series resistor described in Section III.6.3 was connected to the 
inlet of the coolant tube series circuit and the exit of the circuit was fed into the house 
water drain. . Due to the copper construction of the electrode plates, the U-shaped 
coolant tubes provided adequate cooling and no further problems were encountered. 
Note that the coolant circuit described in this paragraph is not shown in Figure 3.26 in 
order to maintain the relative simplicity of the diagram 
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Figure 3.25: Electrical Diagram of Power Supply Configuration Four 
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Figure 3.26: Component Diagram of Power Supply Configuration Four 


III.7 Mounting Configuration 

The annular test sections used in this investigation involved long flow passages 
with very narrow gaps between the independent cylindrical heaters. This geometry 
introduced significant difficulty to the design and construction of an effective mounting 
configuration. Significant effort was expended in the production of two different 
mounting configurations that would allow precise alignment of the multiple test section 
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components used in this research. The base for both mounting configurations was a 
standard triangular optical rail that was bolted to a Unistrut supported particleboard table 
(Figure 3.27, bolts not shown). The custom-made supports used by both mounting 
configurations to secure the Swagelok T-fittings that connected the test section to the 
GTMTF flow loop were attached to the mounting rail using standard, screw-tightened 
optical rail mounts. The custom made supports (Figure 3.28) were constructed in house 
out of steel and aluminum components and incorporated an adjustable, telescoping 
support column. This telescoping support allowed either end of the test section to be 
precisely adjusted to guarantee a level configuration. 


* ALL DIMENSIONS IN MILLIMETERS 

Figure 3.27: Diagram of Triangular Optical Mounting Rail 

The Swagelok T-fittings used in both mounting configurations were standard, 
stainless steel, V” compression fittings and served as the inlet and exit plenums for the 
test section coolant flow. Since these fittings were stainless steel, their connection to the 
GTMTF flow loop effectively grounded the mounting configuration, thus necessitating 
the need for non-conducting inserts that would isolate the inner and outer heaters from 
ground as well as from each other. The inserts used by both mounting configurations 
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also had to maintain the test section in a level, horizontal configuration. To this end, 
Mounting Configurations One and Two used both flexible rubber stoppers and precisely 
machined Teflon inserts to align and isolate the test sections. The following subsections 
characterize the design details of these components. 



SEPARATED VIEW DF STEEL BRACE PLATE AND ALUMINUM BLOCK 



Figure 3.28: Diagram of Custom Made Swagelok T-Fitting Mount 
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III.7.1 Mounting Configuration One 

Mounting Configuration One accommodated the AC powered outer heaters 
(Outer Heaters E and H) and, therefore, was only utilized during OFI experimentation. 
The Teflon insert and rubber stopper shown in Figure 3.29 were inserted into the 
Swagelok T-fittings in order to isolate the inner heater from the outer heater. Note that, 
as described earlier, the copper structure of the outer heater was grounded to the GTMTF 
via the Swagelok T-fittings. The NiChrome wire wrapped around the outer heater 
provided the required outer heater power input. The Teflon sleeve (Figure 3.29, upper 
right comer) was placed on the outer end (away from the heated region) of the Swagelok 
T-fitting. This Teflon sleeve was used to align the inner heater concentrically within the 
Swagelok T-fitting and electrically isolate the inner heater from the outer heater. As 
shown in Figure 3.29, the Teflon sleeve was attached to the Swagelok T-fitting by a 
stainless steel extender piece. This extender was 1.70 inches (43.18 mm) long and was 
sealed by a rubber stopper (Figure 3.29, upper left comer). The primary function of the 
rubber stopper was to prevent leakage from the system, but it also helped support the 
inner heater and maintain its alignment. 

Figure 3.30 shows the entire diagram of Mounting Configuration One. This 
figure clearly shows that the large diameter part of the outer heater was completely 
wrapped with NiChrome heater wire. The steel Heater Wire Clamps held the wire on the 
outer heater and copper Electrode Rings were used to apply power to the outer heater. 
Note that the heated length of the outer heater was measured from the center of the first 
Electrode Ring to the center of the last electrode ring. The brass screws that carried 
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current to the outer heater and attached the Electrode Rings to the outer heater were 
located in the center of the Electrode Rings, which corresponded to the beginning of the 
heated section of the outer heater. 
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Figure 3.29: Diagram of Teflon and Rubber Inserts for Swagelok T-Fitting 
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Figure 3.30: Diagram of Mounting Configuration One 
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III.7.2 Mounting Configuration Two 


Mounting Configuration Two accommodated the DC powered, tube style outer 
heaters (Outer Heaters I - L) used in the highest heat flux OFI experiments and all CHF 
experiments. The Teflon components and the rubber stopper shown in Figure 3.31 were 
inserted into the Swagelok T-fittings in order to isolate the inner heater from the outer 
heater and the outer heater from ground. The angled Teflon piece (Figure 3.31, lower 
right comer) was placed on the outer end (away from the heated region) of the Swagelok 
T-fitting. This angled Teflon component was used to seat the innermost insert securely 
within the Swagelok T-fitting. In addition, the angled Teflon insert incorporated a rubber 
O-ring that helped prevent leakage from the system. The innermost Teflon insert (Figure 
3.31, upper right comer) provided the bulk of the support needed to keep the inner heater 
straight and isolated from the outer heater. As shown in Figures 3.31 and 3.32, this insert 
effectively created the annular flow passage required by this investigation by 
simultaneously supporting the inner heater while providing a concentric seat for the outer 
heater. The final insert placed in the Swagelok T-fitting was a simple rubber stopper 
(Figure 3.31, upper left comer) that fit tightly over the outer heater on the heated side of 
the Swagelok T-fitting. The primary function of this rubber stopper was to maintain the 
watertight integrity of the configuration, but it also helped support the outer heater. 
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.. * all dimensions in millimeters 

VIEW DF ASSEMBLY 

Figure 3.31: Diagram of Teflon and Rubber Inserts for Swagelok T-Fitting 
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Figure 3.32: Diagram of Mounting Configuration Two 
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III.8 Insulation Configuration 


Due to limitations imposed by Mounting Configuration Two, only experiments 
that utilized an AC powered outer heater (Outer Heaters E and H) were insulated. As 
shown in Figure 3.33, Insulation Configuration One consisted of a cylindrical sleeve of 
calcium silicate (maximum temperature of 2300°F / 1260°C) with calcium silicate end 
caps. Each end cap had a one inch (25.4 mm) diameter hole cut in its center to allow for 
connection of the test section to the flow loop and the entire insulation assembly was cut 
lengthwise to allow placement over the test section. Once in place, the insulation covered 
the entire heated length of the test section and was of adequate diameter to clear the outer 
heater electrodes. The insulation was held in place by mounts placed on the optical rail 
described in Section III.7. These mounts were constructed in house and are not shown in 
Figure 3.33. 

Note that the calcium silicate cylindrical sleeve was formed by the manufacturer 
(Thermal Ceramics) using a binder that was only rated to 500°F / 260°C. This binder 
was slowly baked off in a high temperature oven prior to using the insulation in 
experimentation. Previous experimental experience showed that this binder, if not 
properly removed, had a tendency to smoke and catch fire while the test section was 
energized. 
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* ALL DIMENSIONS IN MILLIMETERS 


Figure 3.33: Diagram of Insulation Configuration One 


III.9 Non-Condensable Gas Status 

The amount of non-condensable gas in the coolant was set using the components 
described in Section III.3.6.2 and the procedures outlined in Section IV.3. Each 
experimental setup specifically included the Non-Condensable Gas Status so that 
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comparisons could be made between experiments that utilized degassed water, air- 
saturated water, or nitrogen-saturated water as a working fluid. 


118 



CHAPTER IV 


EXPERIMENTAL PROCEDURES 


This chapter provides a complete description of system preparation and control 
procedures. Many of these procedures are similar to those used by Kennedy, et al. 
(1997), Roach (1998), and Blasick (1999). 


IV.l Initial System Preparation 

As described in Chapter III, the test section used in the GTMTF was comprised of 
two pieces: an inner heater held concentrically within an outer heater. The mounting 
configurations described in Chapter III required leveling and alignment prior to the filling 
of the system. Each of the two braces affixed to the optical mounting rail could be 
swiveled and vertically adjusted independently. In order to ensure that the Swagelok 
fittings, which supported the two pieces of the test section, were aligned and horizontal, a 
half inch diameter solid stainless steel rod was placed in the braces. This solid rod forced 
the four contact points on the two braces to align. Also, a level was placed on the bar in 
order to even out the heights of the two braces. Once in place, the supports on the braces 
were tightened so that the braces would remain stationary. At this point, the Swagelok 
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fittings were placed into the braces and the outer heater was installed. Next, the inner 
heater was carefully fed through the outer heater and centered horizontally using a ruler. 
The appropriate teflon seals were then placed over the ends of the inner heater and the 
Swagelok fittings were tightened into place making the test section leg of the GTMTF 
water tight. The teflon seals, coupled with the teflon supports located within the 
Swagelok fittings, effectively centered the inner heater concentrically within the outer 
heater. Next, electrical isolation between the inner heater, the outer heater, and the metal 
frame (ground) had to be checked. This was accomplished using a Fluke 25 hand held 
multi-meter. Once electrical isolation was ensured, the appropriate power cables were 
attached to the inner and outer heaters. 

The GTMTF was frequently modified in order to repair broken components or to 
change flow conditions. Most modifications required opening the watertight system to 
the atmosphere. After any such occasion, the integrity of the system had to be visually 
checked for leaks before resuming experiments. 


IV.2 System Filling Procedure 

The system filling procedure was nearly identical for all experimental setups. 
Also, the primary components used while filling the system did not change among the 
four flow loop configurations used. Any minor differences in the filling procedure 
associated with a specific flow loop configuration will be clarified as necessary. 
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First, all internal system valves within the flow loop were opened in preparation 
for filling the system. Also, the saturation loop was isolated by closing the saturation 
loop inlet valve (Component SS in Figure 4.1), the saturation loop exit valve (Component 
TT), and the saturation tank pressure regulation valve (Component S). The saturation 
loop had to be isolated since the filling procedure required slight pressurization of the 
system using a peristaltic pump. The saturation tank gas bleed valve (Component HH), 
located within the saturation loop, was always left open to the atmosphere to allow the 
water-saturation gas to bubble through the saturation tank (Component FF) during 
experimentation. This means that any attempt to slightly pressurize the system during 
filling resulted in system overflow through the saturation tank gas bleed valve (located 
ten feet above the ground) if the saturation loop was not isolated. After isolating the 
saturation loop, the system was filled with distilled water using a variable speed 
peristaltic pump (Component AA) attached to the system fill valve (Component BB). 
The pump was switched on and the system relief valve (Component CC) was opened to 
allow air to escape the peristaltic pump line. Once water began flowing out of the system 
relief valve, it was closed and the system fill valve (Component BB) was opened and 
water entered the degassing tank (Component X). At this point, the degassing tank gas 
bleed valve (Component Z) located atop the degassing tank gas separation column 
(Component Y) was opened to allow any air in the degassing tank to be displaced by the 
incoming water. The degassing tank gas bleed valve was left open as water entered the 
system until a steady stream of water could be seen exiting the valve. At this point, the 
degassing tank bleed valve was closed. 
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Figure 4.1: Schematic Diagram of the Georgia Tech MicroChannel Test Facility 













Note that as the system was being filled, all flow loop valves were open and water could 
flow through the entire flow loop during the filling procedure; the saturation loop was 
still isolated at this point. As water continued to flow into the system, the toggle valves 
(Components Q) located on the differential pressure transducer (Component 8) were 
periodically opened to release air pockets from the system. These toggle valves were 
always opened simultaneously to prevent possible damage to the differential pressure 
transducer. In addition, the degassing tank gas bleed valve (Component Z) and 
accumulator bleed valve (Component PP) were opened periodically to release air pockets 
from the system. These four valves were intentionally placed at the highest positions in 
the flow loop in order to capture and bleed off air bubbles. Note that throughout the 
filling process, care was taken to prevent the system pressure from exceeding 20 psig. 
The peristaltic pump used to fill the system was not rated for high pressure and thus could 
have been damaged if forced to operate above 20 psig. 

Once the system was filled and cleared of most air pockets, the four valves 
mentioned above were closed and the system was allowed to pressurize to 20 psig. At 
this point, the peristaltic pump was turned off and the system fill valve (Component BB) 
was closed. Next, either the centrifugal pump (Component A) for Flow Loop 
Configurations 1, 2, or 3, or the positive displacement pump (Component C) for Flow 
Loop Configuration 4 was used to circulate water through the flow loop and cause any 
remaining air pockets to move to the high points in the system. Before switching on the 
pump, the system throttling valve (Component F) was completely closed while the pump 
bypass valve (Component E) was left open. This was done to prevent a sudden surge of 
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water from damaging the turbine flowmeters (Components 3 and/or 4). The pump bypass 
valve was left open to prevent “dead heading” of the pump. Once the pump was 
switched on, the system throttling valve was slowly opened to allow water to flow 
through the loop. The saturation loop was still isolated at this point. 

The pump was turned off after running for approximately five minutes. The 
toggle valves (Components Q), the degassing tank gas bleed valve (Component Z), and 
the accumulator bleed valve (Component PP) were once again opened to allow the 
release of air pockets. Also, the peristaltic pump was turned back on and the system fill 
valve (Component BB) was reopened. After clearing the air pockets from the four bleed 
valves, the system was re-pressurized to 20 psig, the peristaltic pump was turned off, and 
the system fill valve was closed. Next, the system throttling valve (Component F) was 
closed and the appropriate pump (Component A or C) was turned on. The system 
throttling valve was then slowly opened to allow circulation of water through the flow 
loop. Again, this was done to free any trapped air pockets within the system. 

This process of filling the flow loop, running the pump, and bleeding off air 
pockets was repeated until no bubbles could be seen in the flow through the rotameters 
(Components 1 and/or 2). Once the flow loop was completely free of air pockets and 
pressurized once again to 20 psig, the saturation loop was filled. First, the saturation loop 
inlet valve (Component SS) was opened to allow water to flow from the degassing tank 
into the saturation loop. Since the flow loop was pressurized to 20 psig and the saturation 
tank was open to the atmosphere via the saturation tank gas bleed valve (Component 
HH), water flowed into the saturation loop and no additional air passed into the primary 
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flow loop. Next, the peristaltic pump was turned on and the system fill valve 
(Component BB) was opened. Also, the saturation loop exit valve (Component TT) and 
the saturation tank pressure regulation valve (Component S) were opened. Note that the 
saturation tank pressure regulation valve was simply used to set the pressure at the test 
section exit. No water actually flowed through this line during experimental operation, 
but it still had to be filled and cleared of air pockets to ensure accurate experimental 
pressure settings. As the saturation loop filled, air was displaced from the saturation tank 
(Component FF) and released through the saturation tank gas bleed valve. The saturation 
tank blow-off valve (Component LL) located at the top of the saturation tank was also 
opened to release air. Once water began to flow out of the saturation tank blow-off valve, 
the peristaltic pump was turned off and both the system fill valve and the saturation tank 
blow-off valve were closed. At this point, the filling procedure was complete. Note that 
the system could not be pressurized above ambient using the peristaltic pump due to the 
partially opened saturation tank gas bleed valve (Component HH) located atop the 
saturation tank. Continued filling of the system would have lead to overflow through the 
saturation tank gas bleed valve. 


IV.3 Working Fluid Gas Saturation and System Pressurization 

Prior to the start of any set of OFI or CHF experiments, the system was filled / 
cleared of air pockets using the procedure outlined in Section IV.2 above. After the 
system was filled it had to be pressurized to a predetermined level (50, 65, 100, or 150 
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psia exit pressure) and the water had to be saturated with either air or nitrogen. Both of 
these conditions were satisfied using the saturation tank (Component FF of Figure 4.1) 
and a compressed gas supply cylinder (Component JJ). Note that OFI experiments were 
performed using water saturated with either air or nitrogen while CHF experiments 
utilized only nitrogen-saturated water to minimize corrosion of the test section and other 
components. 

First, the saturation tank recirculation pump (Component RR) was started. This 
pump drew water from the degassing tank (Component X) into the saturation tank 
(Component FF). Next, the saturation tank compressed gas supply cylinder (Component 
JJ) was opened to allow compressed air (Setups 1 - 5) or compressed nitrogen (Setups 6, 
8 - 13) to bubble through the saturation tank and saturate the water. Note that the 
saturation tank gas bleed valve (Component HH) was left partially opened to the 
atmosphere at all times. This ensured that compressed gas could continuously bubble 
through the saturation tank and keep the working fluid saturated. After becoming 
saturated, the water passed back into the primary flow loop via the saturation loop exit 
valve (Component TT). Two identical rotameters (Components 12 and 13) allowed 
observation of the flow into and out of the saturation tank. The flow could be adjusted 
and equalized using the various valves available in order to prevent either dryout or 
overflow of the saturation tank. The saturation tank pressure regulation valve 
(Component S) allowed the test section exit pressure to be set via the saturation tank. 
During operation, this valve was left completely open which effectively equalized the 
pressure between the saturation tank and the test section exit. The compressed gas 
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cylinder was fitted with a pressure regulator, which allowed precise pressurization of the 
system to the desired level. At this point the saturation loop was set and the system was 
pressurized for operation. 


IV.4 OFI Primary Flow Loop Preparation 

After saturating the working fluid and pressurizing the system as outlined in 
Section IV.3, the primary flow loop could be started. In preparation for starting the 
primary flow loop for OFI experimentation, the system throttling valve (Component F of 
Figure 4.1) was completely closed and the pump bypass valve (Component E) was 
completely opened. As stated earlier, these actions were taken in order to prevent a 
sudden surge of water from damaging the turbine flowmeters (Components 3 and/or 4) as 
the primary loop circulation pump was started. At this point, the startup procedure varied 
slightly for each of the three flow loop configurations used for OFI experimentation. 


IV.4.1 Flow Loop Configuration One 

Flow Loop Configuration One (shown in Figure 4.2) was used for high flow / 
high heat flux OFI experiments. One multi-stage centrifugal pump (Component A in 
Figure 4.2) was started and the system throttling valve (Component F) was slowly opened 
while the pump bypass valve (Component E) was slowly closed. 
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Figure 4.2: Schematic Diagram of Flow Loop Configuration One 
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At this point, air-saturated water was being pumped from the degassing tank through the 
centrifugal pump, up through a high capacity turbine flowmeter (Component 4) and a 
rotameter (Component 2), through the test section heat exchanger (Component O), 
through the test section, through the system heat exchanger (Component R) and finally 
back to the degassing tank for re-saturation. The flow rate through the test section was 
controlled using the coarse high flow line metering valve (Component K) and the fine 
high flow line metering valve (Component L). Also, the volumetric flow rate through the 
test section was monitored by the high capacity turbine flowmeter (Component 4). At the 
beginning of each OFI experiment that utilized Flow Loop Configuration One, both of 
the coarser and the fine high flow line metering valves were completely open. This 
means that at this point in the setup procedure, air-saturated water was flowing through 
the entire system at full capacity and it was pressurized to the appropriate level for 
experimentation. All system parameters were set for Flow Loop Configuration One 
except for test section inlet temperature which was simply adjusted using three heat 
exchangers (Components O, R, and DD) and/or two degassing tank cartridge style heaters 
(Components EE) as described in Section IV.6. 


IV.4.2 Flow Loop Configuration Two 

Flow Loop Configuration Two (shown in Figure 4.3) was used for low flow / low 
heat flux OFI experiments. One multi-stage centrifugal pump (Component A in Figure 
4.3) was started and the system throttling valve (Component F) was slowly opened while 
the pump bypass valve (Component E) was slowly closed. 
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Figure 4.3: Schematic Diagram of Flow Loop Configuration Two 


COOLANT FLOW 


11 

HH 

GG 


OUT IN 



OUT IN 
COOLANT FLOW 


CC 












At this point, air-saturated water was being pumped from the degassing tank through the 
centrifugal pump, after which point the flow branched. The test section line continued up 
through a rotameter (Component 1) and a low capacity turbine flowmeter (Component 3), 
through the test section heat exchanger (Component O), through the test section, then 
converged with the test section bypass line. The test section bypass line left the 
centrifugal pump, entered a high capacity turbine flowmeter (Component 4), passed 
through a rotameter (Component 2), and went through the test section bypass line before 
converging with the test section flow. The combined flow then passed through the 
system heat exchanger (Component R) and returned to the degassing tank for re¬ 
saturation. The flow rate through the test section was controlled using the low flow line 
metering valve (Component H), while the flow rate through the test section bypass line 
was controlled using the coarse and the fine high flow line metering valves (Components 
K and L). Also, the volumetric flow rate through the test section was monitored by the 
low capacity turbine flowmeter (Component 3) and the volumetric flow rate through the 
test section bypass line was monitored by the high capacity turbine flowmeter 
(Component 4). At the beginning of each OFI experiment that utilized Flow Loop 
Configuration Two, all three flow metering valves (Components H, K, and L) were 
completely open. This means that at this point in the setup procedure, air-saturated water 
was flowing through the entire system at full capacity and it was pressurized to the 
appropriate level for' experimentation. If needed, the two test section bypass line 
metering valves (Components K and L) could be slightly closed in order to increase flow 
through the test section. Care had to be exercised while adjusting these valves since the 
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low capacity turbine flow meter could not exceed a volumetric flowrate of 90 L/hr. Now, 
all system parameters were set for Flow Loop Configuration Two except for test section 
inlet temperature which was simply adjusted using three heat exchangers (Components 
O, R, and DD) and/or two degassing tank cartridge style heaters (Components EE) as 
described in Section IV.6. 

IV.4.3 Flow Loop Configuration Three 

Flow Loop Configuration Three (shown in Figure 4.4) was used for high flow / 
high heat flux OFI experiments and was merely a slight modification of Flow Loop 
Configuration One. Two multi-stage centrifugal pumps (Components A in Figure 4.4) 
were started and the system throttling valve (Component F) was slowly opened while the 
pump bypass valve (Component E) was left open. At this point, nitrogen-saturated water 
was being pumped from the degassing tank through the centrifugal pumps, up through a 
high capacity turbine flowmeter (Component 4) and a rotameter (Component 2), through 
the test section heat exchanger (Component O), through the test section, through the 
system heat exchanger (Component R) and finally back to the degassing tank for re¬ 
saturation. Water was also being pumped directly from the degassing tank, through the 
centrifugal pumps, and back to the degassing tank via the pump bypass line. This was 
necessary due to the high flow rate created by the two centrifugal pumps in series. The 
primary flow loop line could not withstand such a high volume of water. The flow rate 
through the test section was controlled using the coarse and the fine high flow line 
metering valves (Components K and L). 
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Figure 4.4: Schematic Diagram of Flow Loop Configuration Three 
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Also, the volumetric flow rate through the test section was monitored by the high 
capacity turbine flowmeter (Component 4). At the beginning of each OFI experiment 
that utilized Flow Loop Configuration Three, both flow metering valves (Components K 
and L) were completely open. This means that at this point in the setup procedure, 
nitrogen-saturated water was flowing through the entire system at full capacity and it was 
pressurized to the appropriate level for experimentation. All system parameters were set 
for Flow Loop Configuration Three except for test section inlet temperature which was 
simply adjusted using three heat exchangers (Components O, R, and DD) and/or two 
degassing tank cartridge style heaters (Components EE) as described in Section IV.6. 


IV.5 CHF Primary Flow Loop Preparation 

After saturating the working fluid and pressurizing the system as outlined in 
Section IV.3, the primary flow loop could be started. Flow Loop Configuration Four 
(shown in Figure 4.5) was used for all CHF experiments. In preparation for starting the 
primary flow loop for CHF experimentation, the system throttling valve (Component F of 
Figure 4.5) was completely closed and the pump bypass valve (Component E) was 
completely opened. As stated earlier, these actions were taken in order to prevent a 
sudden surge of water from damaging the turbine flowmeter (Component 3) as the 
primary loop circulation pump was started. At this point, the variable speed positive 
displacement pump (Component C) was started. The system throttling valve was then 
slowly opened and the pump bypass valve was slowly closed. Now, nitrogen-saturated 
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water was being pumped from the degassing tank through the positive displacement 
pump, up through a rotameter (Component 1) and a low capacity turbine flowmeter 
(Component 3), through the test section heat exchanger (Component O), through the test 
section, through the system heat exchanger (Component R) and finally back to the 
degassing tank for re-saturation. During CHF experiments, precisely setting the positive 
displacement pump’s rotational speed controlled the coolant flowrate which was 
monitored by the low capacity turbine flowmeter (Component 3). The low flow line 
metering valve (Component H) was slightly closed to provide the positive displacement 
pump with a moderate back pressure (100 kPa). This back pressure helped reduce flow 
oscillations resulting from boiling and two-phase pressure fluctuations encountered 
during CHF experiments. At this point in the setup procedure, nitrogen-saturated water 
was flowing through the entire system at a prescribed flowrate and it was pressurized to 
the appropriate level for experimentation. All system parameters were set for Flow Loop 
Configuration Four except for test section inlet temperature which was simply adjusted 
using three heat exchangers (Components O, R, and DD) and/or two degassing tank 
cartridge style heaters (Components EE) as described in Section IV.6. 
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Figure 4.5: Schematic Diagram of Flow Loop Configuration Four 













IV.6 Test Section Inlet Temperature Control 


The test section heat exchanger (Component O of Figure 4.1) was connected to a 
constant temperature bath which allowed precise setting of the heat exchanger 
temperature; however, the flow rate of the cooling water could not be adjusted. On the 
other hand, both the system heat exchanger (Component R) and the degassing tank heat 
exchanger (Component DD) were attached to the house water line (17° - 20°C) which 
allowed no control of heat exchanger temperature, but permitted adjustment of the 
cooling jacket water flowrate. The system heat exchanger was typically left at its 
maximum capacity due to the large amount of heat it had to remove from the water 
immediately after the test section. The degassing tank heat exchanger was, by design, a 
relatively low capacity heat exchanger even when its cooling jacket water flowrate was 
maximized and thus was used to fine-tune the water temperature as needed prior to being 
pumped through the test section heat exchanger. These three heat exchangers, working in 
concert, were used to control and maintain the test section inlet temperature at the desired 
value for about ninety percent of all OFI and CHF experiments run without the aid of the 
two degassing tank cartridge style heaters. However, during low mass flux OFI and CHF 
experiments, the water was often cooled too much by the system heat exchanger and thus 
had to be reheated prior to being pumped through the test section heat exchanger. The 
two degassing tank cartridge style heaters (Components EE) could be adjusted 
independently in order to accomplish this task. Use of all of these components allowed 
the test section inlet temperature to be set and maintained at the desired level (30°, 50°, 
65°, or 70°C). 
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IV.7 OFI Experimental Procedure and Analysis Techniques 


IV.7.1 OFI Experimental Procedure 

Once the system was running and all experimental parameters were set, an OFI 
experiment could be performed. Each OFI experiment began with the test section 
flowrate maximized and all test section power supplies energized and set to zero. The 
inner heater was consistently more reliable than the outer heater during OFI 
experimentation since it was powered by Power Supply Configuration Two, which 
utilized a DC power supply. For this reason, the surface heat flux of the AC supplied 
outer heater (Power Supply Configuration One) was slowly raised to the prescribed level 
first. Once the outer heater surface heat flux was set, the apparatus was allowed to run 
for five minutes before continuing the startup procedure. The inner heater surface heat 
flux was then slowly raised to the prescribed level. Note that outer heater failure within 
the first five minutes of the OFI experimental startup procedure accounted for most 
experimental failures. Thus, allowing the system to run for five minutes prior to 
increasing inner heater power permitted observation of the outer heater performance. 
Any perceived overheating (which resulted in glowing) of the outer heater triggered 
immediate manual termination of the experiment by the operator since the glowing 
indicated that the heater wire had separated from the outer heater and was about to melt. 
Note that as explained in Chapter III, the power input to the test section during OFI 
experimentation was supplied by two independent sources (Power Supply Configurations 
One and Two) for all but the last five OFI experiments. These final five OFI experiments 
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utilized Power Supply Configuration Three, which used a single DC power supply to 
provide power to both heaters. Thus, for the last five OFI experiments, the surface heat 
flux on both heaters was increased simultaneously to the designated power level for each 
particular OFI experiment with little probability of heater failure. After both heaters were 
set to the appropriate surface heat flux, the DAS system was activated. As explained in 
Section III.3.1.7, the DAS system acquired data from the GTMTF once every five 
seconds and, in turn, recorded the data to a prescribed text file on the Pentium computer 
for future analysis. 

With all systems operating correctly and data being acquired every five seconds 
by the DAS system, the mass flux through the test section was gradually decreased by 
increments of approximately 1-5 L/hr every five to ten seconds until the test section 
differential pressure reached a minimum and began increasing again. This increase in 
differential pressure indicated that the test section had passed through OFI and was now 
operating in the negatively-sloped region of the demand curve. While running the 
experiment, the operator could reduce the flowrate by smaller increments as the 
experiment neared OFI, thus ensuring an accurate reading of the OFI point. 


IV.7.2 OFI Experimental Analysis Techniques 

Once OFI was surpassed, the experiment was terminated and the demand curve 
for the experiment was generated. Once generated, the demand curve directly provided 
the OFI conditions, with OFI defined as the flowrate at which the test section differential 
pressure reached a minimum. Figure 4.6 shows a typical experimental demand curve 
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Figure 4.6: Typical Experimental OFI Demand Curve 


where the OFI point is marked with a solid vertical line. Note that the nominal 
conditions for the experiment are indicated in the upper left comer of the graph and the 
calculated OFI-to-saturation ratios are included in the lower right comer. A plot of this 
type was generated for every OFI experiment conducted; the OFI conditions for all 
experiments conducted in this investigation are listed in Appendix B. 

For each OFI experiment, several experimental parameters were recorded by the 
DAS system described in Section III.3.1.7. These parameters included time, volumetric 
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flow rate through the test section, total power output from each power supply, test section 
inlet and exit temperatures, and test section inlet and exit pressures. Using these 
parameters, inlet and exit water conditions could be determined and used to perform an 
energy balance on the test section. This energy balance served two purposes. First, since 
the power input to the test section was set and maintained at a prescribed value, the 
energy balance served to quantify the percentage of heat lost to the environment. Since 
the OFI point was taken to be a single data point from the experiment instead of an 
average of several data points, and since the power level fluctuated slightly throughout 
each OFI experiment, the heat lost to the environment was needed to determine the exact 
surface heat flux at OFI. Second, the density of the water at the test section inlet was 
required to determine the exact OFI mass flux. The following paragraphs outline the 
energy balance procedure used in analyzing all OFI data. 

Using the inlet and exit temperatures and pressures, the inlet and exit enthalpy and 
fluid density were determined using the Engineering Equation Solver (EES) algorithm 
shown in Appendix A. The nature of OFI experimentation dictated that each experiment 
begin with single-phase liquid water at the test section inlet and exit and end with single¬ 
phase liquid water entering the test section and a two-phase liquid-vapor mixture exiting 
the test section. As such, the exit enthalpy and density found using EES was only 
accurate while single-phase liquid water was exiting the test section. Once the water 
began to boil, the exit temperature and pressure were no longer independent and the exit 
quality was needed to determine exit enthalpy and density. For this reason, the following 
techniques were used to determine the exact surface heat flux at OFI. 
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As stated earlier, the power input to each heater was set and maintained at a 
prescribed level for each OFI experiment. This power level was recorded every five 
seconds by the DAS system and corrected during analysis of the OFI experiment in the 
following manner for the first 61 experiments. First, as stated in Section III.6.1, the 
wattmeter that measured the power supplied to the outer heater from Power Supply 
Configuration One was located before the voltage transformer. From early comparisons 
of the outer heater power level reported by the DAS system and independent 
measurements of voltage and current after the transformer, the transformer was known to 
cause electrical losses of roughly three to five percent. A three percent loss was decided 
upon as a conservative estimate of the power lost through the transformer. This means 
that the outer heater power reported in the DAS data file was subject to electrical losses 
of three percent. Thus, in an effort to equalize the inner and outer surface heat fluxes, the 
outer heater power level was set for each OFI experiment using an independent voltage 
reading taken after the transformer; however, the power level reported in the DAS data 
file had to be corrected. Therefore, the first power correction performed in OFI analysis 
was the subtraction of three percent of the outer heater power level reported in the DAS 
data file. All other power corrections further manipulated this adjusted outer heater 
power level. Note that the DC supplied inner heater was not subject to these losses. 
Second, the test section heat loss to the environment had to be estimated. In order to 
calculate this loss, an energy balance was performed using recorded flow conditions and 
compared to the DAS reported value of total electrical power supplied to the test section. 
As described earlier, the energy balance was only valid while both the test section inlet 
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and exit flow conditions were single-phase liquid. For this reason, the calculated percent 
power loss during each of the first ten single-phase data readings were averaged together 
and applied to the total electrical power reading, i.e., both the inner and outer heater 
power readings, at the OFI point. This provided the actual power absorbed in the 
working fluid at OFI. The average surface heat flux was then calculated by simply 
dividing the total power absorbed in the working fluid by the total heated surface area of 
the test section. This provided the OFI heat flux for the experiment. Note that the last 
five OFI experiments, which utilized Power Supply Configuration Three, were not 
subject to the electrical losses from the transformer as described above. 

The energy balance technique just described enumerated the amount of heat lost 
to the environment from the test section and, thus, had a significant impact on the “exact” 
calculated OFI heat flux for the experiment. A brief error analysis on the calculated heat 
loss to the environment, therefore, was necessary. Table 4.1 shows the effect of changing 
the calculated heat loss by ±50 % to three OFI experiments (selected from the first 61 
experiments) with different nominal surface heat flux values. 


Table 4.1: Effect of Heat Loss Estimation on OFI Calculations 


Nominal 
Heat Flux 

Tinlct 

P exit 

Calculated 
Heat Loss 

Resulting 

Q’WQ”sat 

- 50 % 
Heat Loss 

Resulting 

Q’WQ”sat 

+ 50 % 
Heat Loss 

Resulting 

Q”0Fl/Q”sat 

(MW/m 2 ) 

(°C) 

(kPa) 

(%> 


(%> 


<%) 



65.9 

684.5 

15.7 

0.931 

7.85 

1.018 

23.55 

0.845 

1.5 

50.8 

707.4 

1.9 

0.937 


0.946 

2.85 

0.928 

2.0 

26.6 

685.9 

6.9 

0.924 

3.45 

0.958 

10.35 

0.889 


The OFI correlation developed by this investigation and reported in Section V.1.11 states 
that the average ratio of OFI heat flux-to-saturation heat flux over all experiments was 








0.942 ±9.80 %. This range corresponds to a confidence interval of (0.8497, 1.0343), to 
which all except one of the trial heat loss calculations adhered. This brief error analysis 
of the calculated test section heat loss to the environment demonstrates that even a radical 
miscalculation of heat loss (±50 %) results in a relatively minor error in the resulting OFI 
calculations. 

Power Supply Configuration Three was used to heat the test section during the 
final five OFI experiments and employed a DC wattmeter that measured the total power 
output of the power supply. Heat was lost in two places by Power Supply Configuration 
Three and, thus, a different analysis technique was employed to determine the OFI point 
in these final five experiments. As stated earlier, the DAS system recorded the total 
power supply output in addition to the water conditions at the inlet and exit of the test 
section. While the coolant was single-phase, an energy balance performed on the coolant 
using inlet and exit enthalpy and mass flow rate could be compared to the recorded total 
power supply output in order to determine the percentage of power lost to the 
environment. This loss was then averaged over the first ten single-phase data points and 
applied to the two-phase regime of the experiment. The losses generally accounted for 
approximately 30 percent of the total power supply output. Power was lost primarily in 
two places. First, by design, some power was dissipated in the adjustable series shunt 
resistance used to equalize the inner and outer heater surface heat fluxes. Second, the test 
section was not insulated and thus some heat was lost to the environment directly from 
the test section. From calculations performed using separate current readings and voltage 
taps placed directly on the test section, the percentage of power lost to the environment 


144 



accounted for approximately five percent of the total power supply output, while the 
power lost through the series resistor accounted for roughly 25 percent of the total power 
supply output. In actuality, however, the exact amount of power lost by each of these 
two separate sources was not important since the quantity of interest was the power 
actually absorbed by the working fluid. With the power input to the test section 
evaluated, the average surface heat flux over the test section was found by simply 
dividing the total power absorbed in the working fluid by the total test section heated 
surface area. 

The energy balance technique used to analyze the final five OFI experiments 
enumerated the total amount of heat lost to the environment from both the test section 
and the adjustable series resistor. This calculation, thus, had a significant impact on the 
“exact” calculated OFI heat flux for the experiment. A brief error analysis on the 
calculated heat loss to the environment, therefore, was necessary. Table 4.2 shows the 
effect of changing the calculated heat loss by ±50 % to two OFI experiments (selected 
from the final five experiments) with different nominal surface heat flux values. 


Table 4.2: Effect of Heat Loss Estimation on OFI Calculations 


Nominal 
Heat Flux 

Tinki 

Pexit 

Calculated 
Heat Loss 

Resulting 

Q”0Fl/Q”sat 

- 50 % 
Heat Loss 
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QWQ”« 

+ 50 % 
Heat Loss 

Resulting 

Q’WQ”*u 

(MW/m 2 ) 

(°C) 


(%> 


(%> 


(%> 


2.0 

34.3 


24.4 

0.830 

12.2 

0.964 

36.6 

0.696 

2.5 

47.8 

699.3 

30.7 

0.950 

15.4 

LI 60 

I 46.1 • 

0.739 


This table of heat loss calculations does not adhere to the confidence interval of (0.8497, 


1.0343) as well as the data in Table 4.1 due to the large percentage of heat loss to the 


environment inherently associated with Power Supply Configuration Three. However, 
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this brief error analysis of the calculated test section heat loss to the environment 
demonstrates the gross miscalculations of heat loss that are necessary to skew the results 
of the OFI calculations. 

The OFI mass flux, average exit pressure, and average inlet temperature for each 
OFI experiment were found in the following manner. First, the average inlet temperature 
was simply the mean inlet temperature over the five data points centered on the OFI 
point. Likewise, the average exit pressure was the mean exit pressure over the five data 
points centered on OFI. The standard deviation of each of these mean values was also 
calculated for each experiment and can be found in Appendix B. The OFI mass flux 
(Gofi) for each OFI experiment was found by first calculating the mass flux at each data 
reading taken during the OFI experiment and then picking the value that corresponded to 
the minimum differential pressure over the duration of the experiment. The mass flux at 
each data reading was found by multiplying the volumetric flow rate through the test 
section at that instant by the inlet fluid density and then dividing by the flow area of the 
test section. 

The OFI correlations developed in this investigation compare the OFI mass flux 
to the saturation mass flux and the OFI heat flux to the saturation heat flux. The 
procedures above outline the methods used in determining OFI mass flux and OFI heat 
flux. The following procedures summarize the techniques used in determining saturation 
mass flux and saturation heat flux. They also set forth the methods used in comparing 
OFI conditions to saturation conditions. 

Saturation mass flux simply refers to the coolant mass which corresponds to 
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saturated liquid at the test section exit for a specified set of flow parameters, i.e., inlet 
temperature and pressure, exit pressure, flow area, heated area, and average surface heat 
flux. In order to calculate the saturation mass flux, the inlet enthalpy was found using the 
average inlet temperature described above and the inlet pressure at the OFI point. Note 
that the inlet pressure was not averaged since the enthalpy of a subcooled liquid is a much 
stronger function of temperature than of pressure. Any small pressure fluctuations 
encountered within fifteen seconds of OFI were not significant to analysis calculations 
and any large pressure fluctuations would have triggered fatal experimental errors 
requiring termination of the OFI experiment. Next, the saturation enthalpy of the 
working fluid (hf) at the average exit pressure was found. The difference of the two 
enthalpy values dictated the required change in enthalpy from inlet to exit necessary to 
saturate the working fluid. Finally, the saturation mass flux was calculated using the 
following equation. 


a. 



(4.1) 


Note that Equation 4.1 provided the coolant mass flux required to saturate the working 
fluid at the test section exit given the experimentally fixed value of OFI heat flux. Since 
the OFI mass flux corresponding to the OFI heat flux used in Equation 4.1 was known, a 
direct ratio comparison between the OFI mass flux and the calculated saturation mass 
flux was possible. This provided the means to develop the OFI mass flux correlation 
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located in Chapter V. 

Similarly, saturation heat flux refers to the average surface heat flux required for 
the working fluid to become saturated liquid at the test section exit for given certain flow 
parameters, i.e., inlet temperature and pressure, exit pressure, flow area, heated area, and 
coolant mass flux. Once again, the first step in the saturation heat flux calculation was to 
determine the inlet enthalpy using the average inlet temperature and the inlet pressure at 
the OFI point. Next, the saturation enthalpy of the working fluid at the average exit 
pressure was found. As stated above, the difference of these two enthalpy values dictated 
the required change in enthalpy from inlet to exit necessary to saturate the working fluid. 
Finally, the saturation heat flux was calculated using the following equation. 

= .*!*,]) (4.2) 

^heated 


Equation 4.2 provided the average surface heat flux necessary to saturate the working 
fluid at the test section exit based on the assumption that the experimentally determined 
value of OFI mass flux was actually fixed. Once again, since the OFI heat flux 
corresponding to the OFI mass flux used in Equation 4.2 was known, a direct ratio 
comparison between the OFI heat flux and the saturation heat flux was possible. This 
provided the means to develop the OFI heat flux correlation located in Chapter V. 
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IV.8 CHF Experimental Procedure and Analysis Techniques 


IV.8.1 CHF Experimental Procedure 

Once the system was running and all experimental parameters were set, a CHF 
experiment could be performed. After the DAS system was initialized, each CHF 
experiment began by energizing the power supply and setting its input signal to zero. As 
explained in Chapter III, the power input to the test section during CHF experiments was 
controlled using a Hewlett Packard Data Acquisition / Switch Unit. An input signal was 
provided to the unit from a Compaq DeskPro 575 personal computer (Microsoft 
Windows 95 operating environment) and the HP unit, in turn, provided a signal to the 
power supply, which controlled the power input to the test section. The HP unit also 
monitored the outer heater thermocouples and automatically terminated the experiment 
by interrupting the power supply signal if the outer heater surface temperature exceeded 
300°C. All temperature readings from the outer heater thermocouples were recorded 
along with a real time stamp by the Compaq personal computer for future analysis. 

With the DC power supply energized, the surface heat flux on both the inner and 
outer heaters was slowly raised from zero by increments of roughly 0.01 MW/m 2 per 
minute until CHF was attained. Once a definite pattern for CHF was determined, later 
experiments incremented power more quickly early in the experiment until the heat flux 
exceeded a value equal to two times the predicted OFI heat flux. At that point, the power 
was increased slowly until CHF was attained, after which the power supply was de- 
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energized, the experiment was terminated, and the data were analyzed. For purposes of 
this research, CHF was said to have occurred if either the outer heater or the inner heater 
surface temperature exceeded 300°C. 


IV.8.2 CHF Experimental Analysis Techniques 

Once CF1F was attained, the experiment was terminated and four plots were 
generated for each experiment: inner heater surface temperature versus time, outer heater 
surface temperature versus time, average surface heat flux versus time, and coolant mass 
flux versus time. The two plots of surface temperature were compared with the average 
surface heat flux in order to show the rapid rise in wall temperature associated with CHF. 
Note that CHF experiments were performed using a positive displacement pump making 
the phenomenon of OFI impossible. The predicted OFI value mentioned in Section 
IV.8.1 above was calculated using the empirically developed OFI correlation reported in 
Section V.1.11 and the averaged experimentally measured values of inlet temperature, 
exit pressure, and mass flux. The predicted OFI value was reported on the plots of 
surface wall temperature versus time in order to compare the predicted OFI heat flux 
corresponding to the inlet test parameters with the actual CHF point and to emphasize the 
fact that CHF for a specific set of test parameters occurred at a much higher surface heat 
flux than the predicted OFI heat flux for the same set of parameters. 

For each CHF experiment, several experimental parameters were recorded by the 
DAS system described in III.3.1.7. These parameters included time, volumetric flow 
rate, inner heater wall temperature, total power output from the power supply, test section 
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coolant inlet and exit temperatures, and test section inlet and exit pressures. Using these 
parameters, inlet and exit water conditions could be determined and used to perform an 
energy balance on the test section. Note that the exact CHF heat flux for each experiment 
was not known until this analysis was performed. First, using the inlet and exit 
temperatures and pressures, the inlet and exit enthalpy and fluid density were found using 
the Engineering Equation Solver (EES) algorithm shown in Appendix A. The nature of 
CHF experimentation dictated that each experiment begin with single-phase liquid water 
at the test section inlet and exit and end with single-phase liquid water entering the test 
section and a two-phase liquid-gas mixture exiting the test section. As such, the exit 
enthalpy and density found using EES was only accurate while single-phase liquid water 
was exiting the test section. Once the water began to boil, the exit temperature and 
pressure were no longer independent and the exit quality was needed to determine exit 
enthalpy and density. For this reason, the following modified energy balance technique 
was used to determine the surface heat flux in the two-phase regime. 

As stated earlier, the DAS system recorded the total power supply output in 
addition to the water conditions at the inlet and exit of the test section. While the coolant 
was single-phase, an energy balance performed on the coolant using inlet and exit 
enthalpy and mass flow rate could be compared to the recorded total power supply output 
in order to determine the percentage of power lost to the environment. This loss was then 
averaged over the last ten single-phase data points and applied to the two-phase regime of 
the experiment. The losses generally accounted for approximately 30 percent of the total 
power supply output. Power was primarily lost in two places. First, by design, some 
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power was dissipated in the adjustable series shunt resistance used to equalize the inner 
and outer heater surface heat fluxes. Second, the test section was not insulated and thus 
some heat was lost to the environment directly from the test section. From calculations 
performed on the OFI experiments that utilized Setup 12, the percentage of power lost to 
the environment accounted for approximately five percent of the total power supply 
output, while the power lost through the series resistor accounted for roughly 25 percent 
of the total power supply output. In actuality, however, the exact amount of power lost 
by each of these two separate sources was not important since the quantity of interest was 
the power actually absorbed by the working fluid. 

The energy balance technique used to analyze all 35 CHF experiments 
enumerated the total amount of heat lost to the environment from both the test section 
and the adjustable series resistor. This calculation, thus, had a potentially significant 
impact on the “exact” calculated CHF heat flux for the experiment. An error analysis 
investigating the effect of altering the heat loss values was not performed, however, since 
this CHF investigation was performed only for comparison with OFI and no new CHF 
correlation was developed. Nevertheless, from the heat loss error analysis performed on 
the last five OFI experiments, the estimation of heat loss to the environment allows 
considerable room for error without significant alteration of the overall analysis results. 

With the total power input to the test section evaluated as described above, the 
average surface heat flux over the test section was found by simply dividing the total 
power absorbed in the working fluid by the total heated surface area. This average 
surface heat flux was then plotted against time and compared with the inner and outer 
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heater wall temperatures in order to show the wall temperature spike associated with 
CHF. Figure 4.7 shows a typical plot of heater wall temperature versus time compared 
with average surface heat flux. As stated in Section III.5.5, the outer heater wall 
temperature was measured one centimeter from the end of the test section heated length 
by three separate thermocouples placed 120° apart around the test section. The readings 
from these thermocouples were recorded on a Compaq personal computer via the Hewlett 
Packard Data Acquisition / Switch Unit, not on the DAS system. Since CHF was most 
likely to occur on the outer heater, the HP unit was set to read each thermocouple three 
times per second (the fastest available setting). This sampling rate was much faster than 
that of the DAS system. For analysis purposes, the outer heater wall temperature data 
had to be imported into the Quattro Pro analysis spreadsheet and matched with the DAS 
system real time data column. Once this was done, the data could be easily plotted and 
compared with the average surface heat flux data for the experiment. Note that a 
temperature spike on the upper most thermocouple on the outer heater terminated every 
CHF experiment, so typically this was the only thermocouple data plotted for an 
experiment. 

Figure 4.7 shows that the inner heater wall temperature remained above the 
saturation temperature (164.4°C at 100 psia) of the coolant water for the twenty minutes 
preceding CHF; however, the outer heater wall temperature did not. This fact can be 
directly attributed to the placement of the outer heater thermocouples. The outer heater 
produced heat volumetrically and lost heat to the environment since insulation was not 
used during CHF experiments. These two facts combined to make the outer surface of 
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the outer heater (the surface on which the temperature was measured) the coolest surface 
on the test section. Even though the inner surface temperature of the outer heater was 
inevitably well above the saturation temperature of the coolant water, the outer surface of 
the outer heater registered a temperature that was roughly 20°C less than the saturation 


temperature of the coolant water. Due to the nature of this investigation, precise 
knowledge of the inner and outer heater temperatures was unnecessary. Instead, the 


factor of interest in this investigation was the abrupt spike in the surface temperature of 
either heater associated with CHF; detection of this temperature spike was the sole 
purpose of the four thermocouples employed. 


Heater Wall Temperature vs. Time 




Figure 4.7: Typical Experimental Plot of Heater Wall Temperature Versus Time 
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The average mass flux, inlet temperature, and exit pressure used in the predicted 
OFI calculation outlined above were found in the following manner for each CHF 
experiment. First, the average inlet temperature was simply the mean inlet temperature 
over the duration of the experiment. Likewise, the average exit pressure was the mean 
exit pressure over the duration of the experiment. The standard deviation of each of these 
mean values was also calculated for each experiment and can be found in Appendix C. 
Both the inlet temperature and exit pressure were measured at least once every five 
seconds by the DAS system for each experiment. Averaging these values over the entire 
experiment provided adequate accuracy due to the consistency of the inlet temperature 
and exit pressure maintenance systems. The average mass flux (G) for each CHF 
experiment was found by first calculating the mass flux at each data reading taken during 
the CHF experiment and then averaging these values over the last five minutes of the 
experiment. The mass flux at each data reading was found by multiplying the volumetric 
flow rate through the test section at that instant by the inlet fluid density and then 
dividing by the flow area of the test section. Figure 4.8 shows a typical plot of mass flux 
versus time. This plot reports the nominal conditions of the experiment as well as the 
average and standard deviation of the mass flux. Note that mass flux was not averaged 
over the entire duration of the experiment since it did not remain as steady as inlet 
temperature and exit pressure. Figures 4.9 and 4.10 show typical plots of inlet 
temperature and exit pressure versus time, respectively. These two plots were not 
generated for every CHF experiment performed. They are merely included in this 
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Figure 4.8: Typical Experimental Plot of Coolant Mass Flux Versus Time 
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Figure 4.9: Typical Experimental Plot of Test Section Inlet Temperature Versus Time 
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Test Section Exit Pressure vs. Time 
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Figure 4.10: Typical Experimental Plot of Test Section Exit Pressure Versus Time 


discussion to justify the averaging technique used for test section inlet temperature and 
exit pressure. 

The value considered to be the CHF point for each experiment was actually an 
average of the last five calculated surface heat flux values up to and including the heat 
flux matching the wall temperature spike that terminated the experiment. A standard 
deviation associated with this mean value was also calculated and reported with the data 
in Appendix C. Due to the gradual rate of power increase used during CHF 
experimentation and the relatively fast sampling rate of the DAS system, the standard 
deviations were very low for every CHF value generated (< 1.43 %). These CHF values 
were compared with several CHF correlations in the literature and used to show that OFI, 
rather than CHF, is the limiting phenomenon in microchannel heat transfer. 
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CHAPTER V 


RESULTS AND DISCUSSION 


In this chapter, the experimental results obtained in this investigation are 
presented. Section V.l provides the results for the OFI experiments, while Section V.2 
contains the results of the CHF experiments. 


V.l Onset of Flow Instability 

This investigation, combined with the study performed by Blasick (1999), 
produced a total of 195 OFI experiments. Appendix B lists the OFI conditions for all 
experiments conducted; the data include the non-condensable gas content, the annular 
gap width, the inlet temperature, the exit pressure, the coolant mass flux, the average 
surface heat flux, and the inner-to-outer heat flux ratio. Among the 195 OFI experiments 
listed in Appendix B, 129 experiments were conducted by Blasick, while 66 were 
performed in this investigation. All 195 OFI data points were used to develop an OFI 
correlation that covers a heat flux range of 0.125 - 2.4 MW/m 2 . A brief summary of 
Blasick’s findings will first be discussed, followed by the results obtained here and the 
empirical correlations developed. 
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V.1.1 Brief Discussion of the Findings of Blasick ( 1999 } 

Blasick (1999) investigated the effect of non-condensables on OFI by conducting 
experiments with either a fully-degassed working fluid (zero non-condensable gas 
content) or a working fluid saturated with air at the test section inlet temperature and exit 
pressure. Figure 5.1 shows the effect of non-condensables on OFI for two different heat 
fluxes and three different exit pressures. For each pressure, it can be seen that OFI 
occurred slightly earlier (that is, at a slightly higher mass flux value for a given applied 
heat flux) when the working fluid was air-saturated rather than fully-degassed. However, 
this effect became less pronounced with increased exit pressure because the air in the 
saturated working fluid came out of solution more readily at lower pressures, thus 
increasing the void fraction and leading to premature OFI. Roach (1998) also observed 
this effect while investigating OFI in micro-rod bundle coolant channels. It was 
concluded that the use of air-saturated working fluid provided a more conservative 
estimate of OFI than a degassed working fluid, and would, therefore, be appropriate for 
use in developing design correlations since the amount of non-condensables present in 
practical systems during accident conditions may not be accurately known. Using this 
logic, the majority of experiments in both this investigation and the investigation of 
Blasick (1999) were run using air-saturated water. Figure 5.1 also shows the effect of 
test section exit pressure on OFI. For the same inlet temperature, a reduction in pressure 
resulted in decreased subcooling which caused OFI to occur earlier since OFI is related to 
bulk exit saturation. 
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Unlike the work of Blasick (1999), this investigation did not specifically examine 
the effect of non-condensables on OFI; however, two different gases were used to 
saturate the working fluid (air and nitrogen). As described in Section III.3.6.3, while 
running air-saturated experiments at higher heat fluxes, the inner heaters became pitted 
and showed signs of severe corrosion. In an effort to minimize destruction of the inner 
heaters, the saturation gas was changed from air to nitrogen. With the decrease in 
dissolved oxygen content, the damage to the inner heaters was significantly reduced. 
From Blasick’s results, air saturation was known to produce conservative OFI results as 
compared to fully-degassed conditions; nitrogen saturation was assumed to produce the 
same effect, because of its comparable solubility. 
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o Air Saturated (344 kPa) 

v Air Saturated (687 kPa) 

a Air Saturated (1.03 MPa) 

■ Degassed (344 kPa) 

▼ Degassed (687 kPa) 

A Degassed (1.03 MPa) 


Figure 5.1: Effect of Non-Condensable Gas Content on OFI (from Blasick, 1999) 

Blasick (1999) also investigated the effect of inner-to-outer power distribution on 
OFI by conducting experiments with inner-to-outer surface heat flux ratios ranging from 
0 to infinity. Figure 5.2 shows the effect of power distribution on OFI mass flux at two 
different values of average surface heat flux. Note that the average surface heat flux used 
in Figure 5.2 was not the simple numerical average of the separate heat fluxes applied to 
the inner and outer surfaces. Since the surfaces had different areas, the average heat flux 
was calculated by dividing the total power input by the total heated surface area. The 
resulting value represents the heat flux one would obtain had the same total power input 
been applied to the test section to produce identical heat fluxes on the inner and outer 


161 















surfaces. From the plot in Figure 5.2, it can be seen that decreasing exit pressure results 
in earlier occurrence of OFI, as previously noted. Furthermore, the power distribution 
between the inner and outer surfaces of the annulus evidently had an insignificant effect 
on OFI as long as the total power input was held constant. As a result of this finding, the 
majority of the OFI experiments in this investigation, as well as the experiments 
performed by Blasick (1999), were run with a heat flux ratio of unity (equal heat flux on 
the inner and outer surfaces of the annulus), which is the expected nominal condition for 
the APT target. 



- 

Ratio = 0 (344 kPa) 

- 

Ratio =0.36 (344 kPa) 

Ratio = 1 (344 kPa) 

□ 

Ratio = 8.33 (344 kPa) 

▼ 

Ratio = 0 (687 kPa) 

▼ 

Ratio = 0.36 (687 kPa) 

Ratio = 1 (687 kPa) 


Ratio = 8.33 (687 kPa) 

▲ 

Ratio = 0(1.03 MPa) 

▲ 

Ratio = 0.36 (1.03 MPa) 

« Ratio = 1 (1.03 MPa) 

A 

Ratio = 8.33 (1.03 MPa) 


Figure 5.2: Effect of Inner-to-Outer Surface Heat Flux Ratio on OFI 
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It should be noted that because of test section design constraints, the data of 
Blasick (1999) were limited to moderate heat flux levels up to 1.0 MW/m 2 . The 
remaining OFI subsections of this chapter present the results obtained in this 
investigation, which extended the heat flux range of the Blasick (1999) data to 2.4 
MW/m 2 . 


V.1.2 Effect of Annular Gap Width on OFI 

The effect of annular gap width on OFI was investigated by running experiments 
with annular gap widths ranging from 0.660 mm to 1.001 mm. These experiments 
provided the means to assess the possible effect of manufacturing tolerances on OFI. 
Figure 5.3 shows the effect of annular gap width on the OFI mass flux. It was found that, 
at equivalent heat fluxes, OFI occurred at a slightly lower mass flux in channels with 
larger gap widths (1.001, 0.914, and 0.940 mm). 
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□ 1.001mm -+• 0.914 mm A 0.760 mm x 0.869 mm v 0.940 mm X 0.762 mm 


Figure 5.3: Effect of Annular Gap Width on OFI 


V.1.3 Effect of Inner-to-Outer Heat Flux Ratio and Heated Length on OFI 

Unlike the work of Blasick (1999), this investigation did not specifically examine 
the effect of inner-to-outer surface heat flux ratio on OFI, however, two different surface 
heat flux ratios (zero and unity) were utilized to collect the OFI data located in Table B.2. 
From Blasick’s results, inner-to-outer surface heat flux ratio was known to have an 
insignificant effect on OFI as long as the total power input to the working fluid remained 
the same. For this reason, the bulk of OFI experiments were performed with a heat flux 
ratio of unity, however, while testing the AC powered outer heaters, several experiments 
were run that utilized an unheated inner heater. The experimental values of OFI collected 
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during these experiments were compared with the results reported by Blasick (1999). 
Figures 5.4 and 5.5 show the OFI experimental data collected at a heat flux ratio of zero 
(shown in Table B.2) compared with the OFI heat flux and OFI mass flux correlations 
developed for a heat flux ratio of unity. These two correlations are reported below in 
Section V. 1.11; the zero heat flux ratio data appear to match the correlation’s predictions, 
which supports the conclusion of Blasick (1999). 

Figures 5.4 and 5.5 also demonstrate a more significant and previously 
unobserved effect. All data points included in the two figures were generated using a 
heat flux ratio of zero; however, four different heated lengths were used on the outer 
heater (53, 55, 123, and 180 mm). All data obtained in this investigation for these 
different heated lengths fell along the line generated from correlating the data with a heat 
flux ratio of unity. In other words, both the heat flux ratio and the test section heated 
length had an insignificant effect on OFI, only the total power input to the test section 
was important. 
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V.1.4 Effect of Surface Heat Flux on OFI 

Figure 5.6 illustrates the effect of increasing surface heat flux on OFI. The test 
section inlet temperature and exit pressure were held nominally constant at 30 °C and 100 
psia, respectively, for these three experiments. As expected, the data showed that as the 
surface heat flux was increased, the mass flux at which OFI occurred increased. 



0 0.5 1 1.5 2 2.5 3 


Inlet Velocity (m/s) 



Figure 5.6: Effect of Surface Heat Flux on OFI 
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V.1.5 Effect of Test Section Inlet Temperature on OFI 


Figure 5.7 illustrates the effect of test section inlet temperature on OFI. For a 
given set of parameters, demand curve data were obtained at nominal inlet temperatures 
of 30, 50, and 65 °C. As observed in the figure, increasing the coolant inlet temperature 
resulted in decreased subcooling of the working fluid which caused OFI to occur at a 
higher mass flux. 



V T = 27.0 C □ T = 50.8 C A T = 63.2 C 


Figure 5.7: Effect of Test Section Inlet Temperature on OFI 


V.1.6 Effect of Test Section Exit Pressure on OFI 


Figure 5.8 illustrates the effect of test section exit pressure on OFI. For a given 
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set of parameters, demand curve data were obtained at nominal exit pressures of 50, 100, 
and 150 psia. For the same inlet temperature and surface heat flux values, OFI mass flux 
decreased as the test section exit pressure increased. This result is consistent with the 
findings of Blasick and others. The effect of pressure was quite significant; typically, the 
OFI mass flux for the 150 psia experiments differed from the 50 psia experiments by a 
factor of two, regardless of gap width. 



0 0.5 1 1.5 2 2.5 

Inlet Velocity (m/s) 



Figure 5.8: Effect of Test Section Exit Pressure on OFI 


V.1.7 Comparison with Kennedy’s OFI Model 

Kennedy, et al. (1997) conducted OFI experiments using four test sections with 
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either circular or microbundle interior subchannel cross-sections. Based on the results of 
those experiments, the empirical OFI correlations listed in Table 2.11 were developed. 
Kennedy, et al. (1997) determined that the heat flux necessary to cause OFI was 
approximately 90% of the heat flux required for the working fluid to reach bulk exit 
saturation. Likewise, Kennedy, et al. (1997) found that, for a given heat flux, the mass 
flux at OFI was approximately 1.11 times the mass flux required for the working fluid to 
reach bulk exit saturation. Figures 5.9 and 5.10 compare the OFI data obtained by 
Blasick (1999) and the OFI data from this investigation with the correlations of Kennedy, 
et al. (1997) These figures show that the Kennedy correlations provide an acceptable, 
though slightly conservative, fit to the data. 


Heat Flux Ratio = 1.0 



□ 

Gap = 1.001 mm 


Gap =0.914 mm 

OB 

Gap = 0.760 mm 

X Gap = 0.889 mm 

X 

Gap = 0.940 mm 

A 

Gap = 0.762 mm 

-fr 

Gap = 0.660 mm 



Figure 5.9: Comparison of OFI Heat Flux Data with Kennedy, et al. (1997) OFI Model 
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Heat Flux Ratio = t.O 



□ 

Gap a 1.001 mm 

V 

Gap = 0.914 mm 

OB Gap = 0.760 mm 

X Gap = 0.889 mm 

X 

Gap = 0.940 mm 

A 

Gap = 0.762 mm 

•+- Gap =0.660 mm 



Figure 5.10: Comparison of OFI Mass Flux Data with Kennedy, et al. (1997) OFI 

Model 


V.1.8 Comparison with Roach (19981 OFI Model 

Roach (1998) conducted experiments using the same four test sections as 
Kennedy, et al. (1997) in order to develop the empirical OFI correlations listed in Table 
2.12. Roach (1998) worked primarily with lower mass flux conditions and ultimately 
separated their correlations into subsets of fully-degassed and air-saturated. According to 
Roach (1998), the heat flux necessary to cause OFI under fully-degassed conditions was 
approximately 110% of the heat flux required for the working fluid to reach bulk exit 
saturation; 107% under conditions where the working fluid was air-saturated. Likewise, 
Roach (1998) found that, for a given heat flux, the mass flux at OFI was approximately 
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0.89 times the mass flux required for the working fluid to reach bulk exit saturation under 
fully-degassed conditions; 0.91 times under air-saturated conditions. Figures 5.9 and 
5.10 compare the OFI data obtained by Blasick (1999) and the OFI data from this 
investigation with the air-saturated correlations of Roach (1998). These figures show 
that, for the thin, annular channels studied by Blasick (1999) and this investigation, the 
Roach (1998) low-flow OFI correlations do not provide a good fit to the data; the 
correlations were generally non-conservative. 


Heat Flux Ratio = 1.0 



□ 

Gap « 1.001 mm 

V 

Gap =0.914 mm 

Bl Gap =» 0.760 mm 

X Gap =0.889 mm 

X 

Gap = 0.940 mm 

A 

Gap = 0.762 mm 

Gap = 0.660 mm 



Figure 5.11: Comparison of OFI Heat Flux Data with Roach (1998) Low-Flow, Air- 

Saturated OFI Model 
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Heat Flux Ratio = 1.0 
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Gap = 1.001 mm 


Gap =0.914 mm 

aa 

Gap *0.760 mm 

X Gap = 0.869 mm 

X 

Gap = 0.940 mm 


Gap = 0.762 mm 

-f- 

Gap - 0,660 mm 



Figure 5.12: Comparison of OFI Mass Flux Data with Roach (1998) Low-Flow, Air- 

Saturated OFI Model 


V.1.9 Comparison with the Saha and Zuber (1974) QSV Model 

Next, the OFI data were compared with predictions of the Saha-Zuber OSV 
model. Correlation against the Saha-Zuber OSV model gave the following results. 

q'' OFI =l.000*q" s _ z ±n.93% (5.1) 

G ofi = 0.998 * G s _ z ±15.65% (5.2) 

This suggests that the Saha-Zuber OSV model can be accurately used to predict OFI for 
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the thin, annular channels used in this investigation over the experimental range of 
parameters, with appropriate safety factors given the uncertainty of the correlations. 
Figures 5.13 and 5.14 show graphic representations of the OFI data against OSV 
predictions by the Saha-Zuber correlations. 


Heat Flux Ratio « 1.0 



I—I 1.001mm T7 0.914 mm K 0.760 mm X 0 889 mm X 0.940 mm ^ 0.762 mm 0.660 mm 


Figure 5.13: Correlation of OFI Heat Flux Data Against Saha & Zuber OSV Model 
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Heat Flux Ratio = 1.0 



r~1 1.001 mm V 0.914 mm m 0.760 mm X 0.889 mm X 0.940 mm ^ 0.762 mm -+- 0.660 mm 


Figure 5.14: Correlation of OFI Mass Flux Against Saha & Zuber OSV Model 


V.1.10 Comparison with the Unal (19751IPNVG Model 

Using a similar approach, correlation against the Unal (1975) IPNVG model 
yielded the following. 


Q"ofi = 1-082 * q" Unal ±\2>.\S% (5.3) 

G 0FI = 0-922* Gy,±17.50% (5.4) 

This relationship indicates that the Unal (1975) IPNVG model provides a conservative 
estimate of OFI. Again, graphic representations of the OFI data against these correlations 
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are shown in Figures 5.15 and 5.16. As with the Saha-Zuber correlations above, the 
uncertainty range (±) listed in the correlation equals twice the standard deviation of the 
percent error between the prediction and the raw OFI data. 


Heat Flux Ratio = t.O 



□ t oot mm V 0.914 mm 00 0.760 mm X 0.889 mm X 0.940 mm 0.762 mm -+- 0.660 mm 


Figure 5.15: Correlation of OFI Heat Flux Data Against Unal (1975) IPNVG Model 
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Figure 5.16: Correlation of OFI Mass Flux Data Against Unal (1975) IPNVG Model 


V.1.11 Development of OFI Correlation 

As described in Section IV. 7 . 2 , with OFI mass flux, OFI heat flux, inlet enthalpy, 
and exit pressure known for each experiment, an energy balance was used to calculate the 
corresponding heat flux and mass flux that would lead to saturation at the test section 
exit: 


_ it heated 

"' ° n 

= a m (h, [ 4 , ] - h[T M ,„p lnl , t ]) 

iheated 


(5.5) 


(5-6) 













Based on these energy balances, simple correlations were developed for the OFI heat flux 
(or mass flux) as functions of the corresponding saturation heat flux (or mass flux). 
Using this approach, Kennedy, et al. (1997) and Roach (1998) developed'correlations for 
OFI representing high and low mass fluxes through circular microchannels and micro-rod 
bundle coolant subchannels, respectively. The data obtained in this investigation were 
augmented by the data produced by Blasick (1999) and used to develop the following 
correlations. Note that only the experimental data with an inner-to-outer heat flux ratio 
of unity were used, since this was the nominal condition of the APT target. For a given 
mass flux or heat flux, the following correlations were developed to predict the 
corresponding heat flux or mass flux at OFI: 

<l" 0F i = 0-942 * tf" ja ,±9.80% (5.7) 

G ofi = 1.064 * G sat ± 11.55% (5.8) 

These correlations were obtained by minimizing the sum of the squares of the 
absolute error (i.e., the difference between the experimental and predicted values of q’oFi 
and Gofi- These correlations predict the data within ± 9.80% for OFI heat flux given a 
known mass flux and ± 11.55% for OFI mass flux given a known heat flux, with a 95% 
confidence level. The uncertainty ranges equal twice the standard deviation of the 
percent error between the prediction and the raw OFI data. Figures 5.17 and 5.18 
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compare predictions of the above correlations with the data. The correlations are valid 
for the following parameter ranges: 


Heat Flux: 

Mass Flux: 

Inlet Temperature: 

Exit Pressure: 

Non-Condensable Gas Content: 
Inner-to-Outer Heat Flux Ratio: 
Coolant Channel Geometry: 
Coolant Channel Gap Thickness: 


0.125 2.4 MW/m 2 

85 -» 2900 kg/m 2 sec 
30 -» 70°C 
50 -> 150 psia 
Saturated at Ti n i e t and Pexit 
1.0 

Annulus 

0.660 —> 1.001 mm 


Heat Flux Ratio = 1 



□ 

Gap = 1.001 mm 

V 

Gap =0.914 mm 

X 

Gap = 0.760 mm 

5£ Gap = 0.889 mm 

a 

Gap = 0.940 mm 

A 

Gap = 0.762 mm 

*+■ 

Gap = 0.660 mm 



Figure 5.17: OFI Heat Flux Correlation 
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Heat Flux Ratio = 1 



□ Gap * 1.001 mm 
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Gap = 0.914 mm 
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Gap =0.760 mm 

X Gap = 0.889 mm 

OB Gap = 0.940 mm 


Gap = 0.762 mm 

H- 

Gap =0.660 mm 



Figure 5.18: OFI Mass Flux Correlation 


V.2 Critical Heat Flux (CHF) 

In this section, the results of the CHF experiments conducted in this investigation 
are presented. A total of 35 low mass flux CHF experiments were conducted. The 
primary goal of these experiments was to show that OFI, rather than CHF, is the limiting 
phenomenon in microchannel heat transfer. Appendix C lists the CHF conditions for all 
experiments conducted; the data include the non-condensable gas content, the annular 
gap width, the inlet temperature, the exit pressure, the coolant mass flux, the average 
surface heat flux, and the calculated exit quality. 


180 




















V.2.1 Effect of Coolant Mass Flux on CHF 


Figure 5.19 illustrates the effect of increased coolant mass flux on CHF. In this 
figure, the test section inlet temperature and exit pressure were nominally maintained at 
30°C and 100 psia, respectively, while the coolant mass flux was varied. The data 
showed that, as coolant mass flux was increased, the surface heat flux at which CHF 
occurred increased. Roach (1998) also observed this intuitive result after performing a 
similar investigation of CHF in circular microchannels and in a micro-rod bundle 
subchannel. 



Figure 5.19: Effect of Coolant Mass Flux on CHF 
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V.2.2 Effect of Test Section Inlet Temperature on CHF 


Figure 5.20 illustrates the effect of test section inlet temperature on CHF. For a 
given set of parameters, CHF experiments were conducted at nominal inlet temperatures 
of 30, 50, and 65 °C. As observed in the figure, increased inlet temperature resulted in 
decreased subcooling of the working fluid which caused CHF to occur at a lower heat 
flux. The effect appears to become more significant at higher mass flux values. 



V Inlet T = 30 C O lnletT = 50C A lnletT=65C 


Figure 5.20: Effect of Test Section Inlet Temperature on CHF 


V.2.3 Effect of Test Section Exit Pressure on CHF 

Figure 5.21 illustrates the effect of test section exit pressure on CHF. For a given 
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set of parameters, CHF experiments were conducted at nominal exit pressures of 50, 100, 
and 150 psia. With inlet temperature and coolant mass flux set, the CHF heat flux 
increased as the test section exit pressure increased. This result was consistent with the 
findings of Roach (1998) and can be directly attributed to the higher saturation 
temperatures associated with higher system pressures. The effect appears to be more 
significant at elevated mass flux values. 



Figure 5.21: Effect of Test Section Exit Pressure on CHF 



Next, the CHF data were compared with predictions of the Caira, et al. (1993) 
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empirical CHF model (Table 2.15). Correlation against the Caira CHF model gave the 
following results. 


<ta* = 0.785 *q" Caira ±24.03% (5.9) 

This suggests that the Caira CHF model is not a conservative predictor of low mass flux 
CHF in horizontal, thin annuli. The large uncertainty value associated with the 
comparison (± 24.03%) detracts from the possible use of Equation 5.9 for CHF 
predictions (Figure 5.22). 



Figure 5.22: Correlation of Experimental CHF Data Against Caira, et al. (1993) CHF 

Model 
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V.2.5 Comparison with the Bowring (1972) Empirical CHF Correlation 


Next, the CHF data were compared with predictions of the Bowring (1972). 
empirical CHF model (Table 2.14). Correlation against the Bowring (1972) CHF model 
yielded the following results. 

q "chf = 0.586 * q" Bowring ±16.28% (5.10) 

As with the Caira, et al. (1993) comparison. Equation 5.10 suggests that the Bowring 
(1972) CHF model is not a conservative predictor of low mass flux CHF in horizontal, 
thin annuli. The lower uncertainty values associated with the comparison (± 16.28%) 
suggest that the model represented by Equation 5.10 may be used for design purposes 
(Figure 5.23). 
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Figure 5.23: Correlation of Experimental CHF Data Against Bowring (1972) CHF 

Model 

V.2.6 Comparison with Predicted OFI Correlation 

As stated earlier, the primary goal of this CHF study was to show that OFI, rather 
than CHF, was the limiting heat transfer phenomenon in microchannel heat transfer. 
Figure 5.24 shows a plot of experimentally determined CHF surface heat fluxes versus 
the corresponding OFI values predicted using Equation 5.7. The experimental conditions 
for the CHF experiments listed in Appendix B. A least square fit of the CHF data and 
the corresponding predicted OFI values produced the following correlation. 
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q" CHF = 3.047 * q" 0FI ± 12 . 63 % 


( 5 . 11 ) 


Figure 5.25 shows a plot of the experimental CHF data against this correlation. Equation 
5.11 and Figure 5.24 clearly show that for a given set of conditions, CHF occurs at a 
surface heat flux of more than three times the predicted OFI heat flux, thus proving that 
OFI, not CHF, is the limiting heat transfer phenomenon in microchannel heat transfer. 



Figure 5.24: Experimental CHF Data Versus Predicted OFI Heat Flux 
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CHAPTER VI 


CONCLUSIONS 


Ledinegg flow instability and critical heat flux experiments were conducted with 
water flowing through thin, concentric annuli in thirteen different experimental setups 
using the Georgia Tech MicroChannel Test Facility. Subcooled water at a controlled flow 
rate, pressure, temperature, and dissolved non-condensable gas concentration was 
delivered to the annular test section. The inner and outer surfaces of the annulus were 
heated electrically up to 2.5 MW/m 2 using four different variable power supply 
configurations. A total of 66 OFI and 35 CHF experiments were conducted in this 
investigation and combined with the 129 OFI experiments conducted in the investigation 
of Blasick (1999) for analysis and comparison with known two-phase models. The tables 
located in Appendices B and C provide all of the OFI and CHF data collected and used 
during this study. 


VI. 1 OFI Conclusions 

1. Inner-to-outer surface heat flux ratio was found to have an insignificant effect on 
OFI, only the total power input to the test section was important. This finding is 
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consistent with the conclusions of Blasick (1999). 

2. The test section heated length was found to have an insignificant effect on OFI. 
Once again, only the total power input to the test section was important, not the 
method of application. 

3. With all other operating parameters constant, OFI occurred at higher mass fluxes 
as the surface heat flux applied to the test section was increased. 

4. With all other operating parameters constant, OFI occurred at higher mass fluxes 
as the test section inlet temperature was increased. 

5. With all other operating parameters constant, OFI occurred at higher mass fluxes 
as the test section exit pressure was decreased. 

6. With all other operating parameters constant, OFI occurred at higher mass fluxes 
as the test section annular gap width was decreased. 

7. An integral OFI correlation was developed relating the heat flux required to cause 
bulk exit saturation of the working fluid with the heat flux required to initiate 
OFI. A second integral OFI correlation was developed for the coolant mass flux 
data. 


? V, = 0.942V',.,±9.80% (6.1) 

G on -1.064 •<?„ ±11.55% (6.2) 

The data fit each of these correlations within a 95% (2a) confidence level. The 
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parameter ranges used to develop these correlations were: 


Heat Flux: 

Mass Flux: 

Inlet Temperature: 

Exit Pressure: 

Non-Condensable Gas Content: 
Inner-to-Outer Heat Flux Ratio: 
Coolant Channel Geometry: 
Coolant Channel Gap Thickness: 


0.125 -> 2.4 MW/m 2 
85 —» 2900 kg/m 2 sec 
30 —> 70°C 
50 -> 150 psia 
Saturated at Ti n i et and P ex it 
1.0 

Annulus 

0.660 ->1.001 mm 


8. The Saha and Zuber (1974) OSV model gives an accurate fit to the data while the 
Unal (1975) IPNVG model provides a slightly conservative estimate of OFI. 

VI.2 CHF Conclusions 

1. With all other operating parameters constant, CHF occurred at lower heat fluxes 
as the coolant mass flux supplied to the test section was decreased. 

2. With all other operating parameters constant, CHF occurred at lower heat fluxes 
as the test section inlet temperature was increased. 

3. With all other operating parameters constant, CHF occurred at lower heat fluxes 
as the test section exit pressure was decreased. 

4. Neither the Caira, et al. (1993) CHF model nor the Bowring (1972) CHF model 
provides an accurate fit to the experimental CHF data. In addition, correlation of 
the experimental data against both models produces a relatively large amount of 
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error (24.03% and 16.28%). Of the two, however, the Bowring (1972) correlation 
produced the least amount of error, which lends to the conviction with which the 
Bowring (1972) correlation, with a multiplier of 0.586, can be used to predict 
CHF in thin, annular channels over the parameter range used in this investigation. 

5. Finally, OFI, rather than CHF, is the limiting heat transfer phenomenon in 
microchannel heat transfer. Comparison of the CHF experimental data to the OFI 
correlation shown in Equation 6.1 showed that, for a given set of operating 
conditions, CHF occurs at a surface heat flux of more than three times the 
predicted OFI heat flux. 

q" chf = 3.047 * q" 0FI ± 12.63% (6.3) 
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APPENDIX A 


EES PROPERTIES ALGORITHM 


This appendix lists the algorithm used in Engineering Equation Solver (EES) to 
determine the coolant properties at the test section inlet and exit. These properties were 
the viscosity, density, enthalpy, thermal conductivity, specific heat, and saturation 
temperature of the working fluid. For a particular set of OFI data, analysis began by 
copying each combination of inlet temperature and inlet pressure from the Quattro Pro 
analysis spreadsheet into an EES parametric table. The EES algorithm below was then 
run and the six fluid properties mentioned earlier were generated for each combination of 
temperature and pressure. These values were then copied from EES into the Quattro Pro 
analysis spreadsheet. An identical procedure was followed in order to determine the 
coolant properties at the test section exit. 

Two important facts about the algorithm must be noted. First, the pressure read 
by the DAS system and reported in the Quattro Pro analysis spreadsheet was gauge 
pressure. The EES algorithm below converted this gauge pressure into absolute pressure 
by adding 101.325 kPa to the value copied from the Quattro Pro analysis spreadsheet. 
Second, EES required input values of pressure and temperature for the properties 
calculated. Since both OFI and CHF experimentation required operation in a two-phase 
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flow regime, temperature and pressure were not independent during a portion of each 
experiment and, therefore, temperature had to be altered in order to determine the coolant 
properties. To this end, 0.5°C was subtracted from each data point during analysis. 


EES Algorithm 

viscosity=viscosity(WATER,T=T_corrected, P=P+101.325) 
rho=density(WATER,T=T_corrected, P=P+101.325) 
h=enthalpy(WATER,T=T_corrected, P=P+101.325) 
k=conductivity(WATER,T=T_corrected, P=P+101.325) 
cp=1000*specheat(WATER,T=T_corrected, P=P+101.325) 
temp_sat=temperature(WATER,P=P+101.325, x=1) 
T_corrected=temp-0.5 
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APPENDIX B 


SUMMARY OF OFI EXPERIMENTAL DATA 


The conditions corresponding to OFI for all experiments conducted by Blasick 
(1999) and in this investigation are tabulated in Tables B.l and B.2, respectively. The 
non-condensable gas content, annular gap width, inlet temperature, sample standard 
deviation of the inlet temperature, exit pressure, OFI mass flux, average surface heat flux, 
and inner-to-outer surface heat flux ratio are given in Table B.l. Table B.2 includes this 
same information along with the sample standard deviation of the exit pressure. 


Table B.l: OFI Experimental Data from Blasick (1999) 


Non-Condensable 
Gas Content 

Gap Width 

Tinlet 


P exit 

(absolute) 

OFI Mass Flux 

OFI Heat Flux 

Inner-to-Outer 
Heat Flux Ratio 


(mm) 

MifaMi 

mssm 

(kPa) 


(MW/m 2 ) 


Air Saturated 

1.001 


TiTiTif 

345 

356.54 

0.2501 

0 


1.001 


KEZEH 

348 



0.36 

HE3BE333S9I 

1.001 

50.5 

KEESI 

346 

■ 


1 

Air Saturated 

1.001 

50.0 

0.051 

346 

189.11 

0.1254 

1 

Air Saturated 

1.001 

49.5 

MT-TM 

346 

155.41 

0.1250 

8.33 

Air Saturated 

1.001 

49.6 

BilIM 

694 



0 


1.001 

51.1 

0.101 

686 

222.72 

0.2493 

0.36 


1.001 


0.143 

690 

248.12 

0.2505 

1 

Air Saturated 

1.001 

49.6 

0.094 

689 

108.42 

0.1240 

1 

Air Saturated 

1.001 

49.1 

TilHM 

693 

92.80 

0.1245 

8.33 

Air Saturated 

1.001 

50.2 


1041 

173.49 

0.2510 

0 

Air Saturated 

1.001 

50.1 

WEEHM 

1035 

166.25 

0.2506 

0.36 

Air Saturated 

1.001 

49.6 

KX&UH 

1026 

190.34 

0.2495 

1 

Air Saturated 

BHESOH 

50.5 

0.246 

1031 

219.16 

0.3062 

0 

Air Saturated 

HDES0HI 




245.75 


0.28 

Air Saturated 

1.001 

49.7 


1037 


0.3076 

0.71 

Air Saturated 

1.001 

49.1 

0.052 

1034 

100.04 

0.1254 

1 
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Table B.l (continued) 


Non-Condensable 
Gas Content 


Gap Width Timet a (Time) 


OFI Mass Flux 

OFI Heat Flux 

Inner-to-Outer 
Heat Flux Ratio 



Air Saturated 
Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


Fully Degassed 
Fully Degassed 


Fully Degassed 




Fully Degassed 


Fully Degassed 
Fully Degassed 
Fully Degassed 


Fully Degassed 
Fully Degassed 


Fully Degassed 
Fully Degassed 
Fully Degassed 


_F 

F 


■BB^I 

mmmi 

Mnll ' ll l Jii-UW-f I 


2 ! 

Air Saturated 
Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 
Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 



Air Saturated 
Air Saturated 


Air Saturated 
Air Saturated 


49.7 


49.2 

_ 


29.9 


68.8 


49.9 


50.1 


50.2 




0.168 

69 


0.052 

34 


303.58 


4 

696 

5 

730 


89.11 







49.8 


49.3 
















































































































































Table B.l (continued) 


Non-Condensable 
Gas Content 

Gap Width 

TinJct 

^(Tiijiet) 


OFI Heat Flux 


Inner-to-Outer 
Heat Flux Ratio 



Air Saturated 
Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 
Air Saturated 
Air Saturated 


1024.11 


923. 

1063 


1305.94 
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Table B.2: OFI Experimental Data (Stoddard) 


Non-Condensable 
Gas Content 

Gap 

Width 

T inlet 

O’(Tinlct) 

P exit 

(absolute) 

tf(Pexit) 

OFI Mass Flux 

OFI Heat Flux 

Inner-to-Outer 
Heat Flux Ratio 


(mm) 

(°C). 

>C) 

(kPa) 

(kPa) 

(kg/m 2 s) 

(MW/m 2 ) 



Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


0.940 49.3 


0.940 

0.940 


0.940 56.5 


0.940 


0.940 


0.940 


0.940 




1 696 

L_ 


Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


] Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 



690 

Air Saturated i 

(0.940 

iityp 


690 


j 

344 


n 

1034 



Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 
Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


Air Saturated 


Nitrogen Saturated 


Nitrogen Saturated I 0.940 






1350.83 


869.17 


1402.13 


779.70 


907.41 


1 

4.143 

9 

2.697 

3H 

1.695 


1377.32 

1 

1048.31 

1 

1085.26 

1 

1426.83 

1 

1633.30 

1 

1298.76 

1 

746.70 

1 

1663.18 

1 

1683.66 

1 


Nitrogen Saturated 


Nitrogen Saturated 





1 


1 

i 

1 


























































































































































Table B.2 (continued) 


Non-Condensable 
Gas Content 

Gap 

Width 

Tinlct 

O’CTjnict) 

Pexit 

(absolute) 

a(P«it) 

OFI Mass Flux 

OFI Heat Flux 

Inner-to-Outer 
Heat Flux Ratio 


(mm) 

■SSI 

WK5M 

(kPa) 

(kPa) 

(kg/m 2 s) 



Nitrogen Saturated 

0.762 

EQ 

0.050 

997 

2.892 

0.8491 

683.34 

1 

Nitrogen Saturated 

0.762 

wm 

0.083 

698 

4.309 

1.2181 


1 

Nitrogen Saturated 


mi 

0.212 

696 

if 

1.2185 

1251.79 

1 

Nitrogen Saturated 

WSEM 


0.189 



1.2088 

1499.75 

1 

Nitrogen Saturated 

■sKSS 

51.2 


1002 

MEM 

1.1909 

1030.09 

1 

Nitrogen Saturated 

mui.ym 

29.7 


697 

II.T.Tf 

1.4607 

1222.95 

1 


WEEM 

eei 

Bra 

698 


1.4620 

1407.01 

1 


WSEM 

KM 

HKirngjl 

697 

■*li» 

1.4527 

1682.18 

1 

mmmkmmm 

0.762 

48.9 

0.344 

990 


1.4691 

1259.97 

1 

uummmmm 

0.660 

34.3 

0.052 

694 

KSZQ1 

1.7302 

2086.20 

1 

BHI 

0.660 



691 

0.809 

1.7519 

2251.86 

1 

Nitrogen Saturated 1 

0.660 

32.6 


699 

■m 

2.3419 

2762.58 

1 


■iliMfl 



699 

0.632 

2.4177 

2897.07 

1 


jjl]TiW«i 

M&a 

0.000 I 

986 

0.786 

1.9925 

2008.69 

1 J ^ 
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APPENDIX C 


SUMMARY OF CHF EXPERIMENTAL DATA 


The conditions corresponding to CHF for all experiments conducted in this 
investigation are tabulated in Table C.l. The non-condensable gas content, annular gap 
width, inlet temperature, sample standard deviation of the inlet temperature, exit pressure, 
sample standard deviation of the exit pressure, coolant mass flux, sample standard 
deviation of the coolant mass flux, observed critical heat flux, sample standard deviation 
of the critical heat flux, and the calculated exit quality are given in Table C.l. 


Table C.l: CHF Experimental Data (Stoddard) 


Non-Condensable 
Gas Content 

Gap 

Width 

Tinlct 

O’(Tinlct) 

Pexit 

(absolute) 

CT(Pexii) 

Mass 

Flux 

ct(Gchf) 

Heat Flux 

<r(q”cHF) 

Xcxit 


(mm) 

■831 

WB3M 

(kPa) 

(kPa) 

(kg/m 2 s) 

(kg/m 2 s) 

(MW/m 2 ) 

(MW/m 2 ) 

(%> 

MJU.J.LL.I.'U.ILUgM 

0.660 

mm 

0.231 

690 


Wrlt'ihlM 

15.66 

0.4689 

0.0067 

KEBU 

Nitrogen Saturated 

0.660 

Egj 

0.206 

694 

msm 

■BEKMI1 

2.44 


0.0031 

EESBI 

Nitrogen Saturated 

0.660 

■am 

K TRIM 

688 

law 

■i:vjcgai 

1.39 


0.0002 

43.2 

Nitrogen Saturated 

0.660 

E3a 


344 

hues 


3.25 

Mittal 

0.0004 

35.0 

Nitrogen Saturated 

EI331 

50.2... 

■es 

1033 

wmm 

umm 

2.68 

0.4986 


W5EM 


BSEB1 

ESI 


689 

■ 

mu&m 

2.91 

0.3167 

■TTJZM 


Nitrogen Saturated 


Bill! 

0.112 

690 


100.68 

0.86 

■EEESH 



Nitrogen Saturated 

IlMf 

E9 

0.125 

692 

if 

99.19 



■ 

ESI 

Nitrogen Saturated 


49.6 

0.175 

354 



2.60 

0.2392 

0.0002 

44.4 1 


wmm 

■aw 


1037 



1.31 

MiJkWiJJB 

0.0014 

Eg 

Nitrogen Saturated 

■iKil 

»H:1 


695 



12.81 


0.0024 

1W 

Nitrogen Saturated 


49.7 


687 

3.976 

290.25 




KMW 

Nitrogen Saturated 


65.8 






Mil lif'iW 


WSSM 

Nitrogen Saturated 

0.660 

51.4 


345 



7.18 


0.0038 

38.6 

Nitrogen Saturated 

0.660 

49.2 


1039 



10.88 


0.0032 

47.9 



29.8 


692 

2.838 

151.16 

3.34 

0.3958 

0.0012 

44.3 

Nitrogen Saturated 

feMsI 

49.8 


689 

5.627 

146.59 

2.39 

0.3383 

0.0018 

39.9 

Nitrogen Saturated 

0.660 

65.3 

0.467 

691 

3.405 

150.97 

5.09 

0.3419 

0.0018 

41.8 
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Table C.l (continued) 


Non-Condensable 
Gas Content 

Gap 

Width 


IJjjggjj 



■E9SB 

ct(Gchf) 

Heat Flux 

Gr(q’W) 




(mm) 



BESBM 



(kg/m 2 s) 

(MW/m 2 ) 

(MW/m 2 ) 

(%) 

Nitrogen Saturated 

0.660 


0.639 

351 


■ttllEM 

2.67 

0.3031 


34.6 


0.660 

Hdl 

0.149 

1033 

■3HEM 


1.80 

0.3961 


gjllM 



KiEfl 

WWW 

694 



9.33 


MCliTilTS 

43.6 



49.5 

M»n$ai 

688 



10.64 



41.2 


| Nitrogen Saturated 

0.660 

65.6 j 

1.356 

689 

mm 

WrViWhM 

3.78 

HBM 

■ ■ 

HI! 



MMU 

1.128 

349 

mm 


7.21 



Bd'VH 









8.32 



EXB 

Nitrogen Saturated 

0.660 

EH 

SEEEBKi 

700 






g Eg 

Nitrogen Saturated 


EE1 


684 


331.58 

MW 

■ifcfcEkl 

0.0044 

45.9 

Nitrogen Saturated 



■SE9 

710 


328.99 

M5BM 

■iUrHftl 

0.0046 

40.6 

Nitrogen Saturated 

0.660 


Wi!M 

343 




0.7716 

0.0041 

43.1 

Nitrogen Saturated 

0.660 

EEI 


1029 



WBHM 


HEEES3H 


Nitrogen Saturated 


WthW 


705 

■HMW 

MBIM 

■119; 

■SE 

0.0022 

■aiKJ 



ea 

EESfl 


mm 

gViWUH 

MW 

Mil TTTliM 

0.0044 


Nitrogen Saturated 

0.660 


MiliWiM 




■ESEH 


BHB 

MVJM 

Nitrogen Saturated 

0.660 

BH 

H.W1 

343 

9.894 

ebb 

H9I 

mssmm 

hrb 

esb 


0.660 

Wb*m 

1.537 

1031 

3.737 

379.39 | 

13.12 

1.0686 | 

0.0034 

MIM 
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APPENDIX D 


OFI ERROR ANALYSIS 


This appendix presents an error analysis of the experimentally obtained OFI mass 
flux and OFI heat flux values for all experiments performed in this investigation. Blasick 
(1999) performed a similar analysis after which this examination is modeled. 


D.l Mass Flux Error Analysis 

Mass flux can be calculated as a function of the coolant mass flowrate and the test 
section cross-sectional area. 


G 



(D.l) 


Likewise, the variance of the coolant mass flux can be expressed by the following error 
propagation formula. 



(D.2) 


where 
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dG _ J_ _5G_ = __m_ 

am 4, 34 4 


(P-3) 


Since the variance of the coolant mass flowrate and the test section cross-sectional area 
were not directly known, these quantities were expanded into known components until all 
variances were in terms of experimentally measured variables (i.e., volumetric flowrate, 
power, test section diameter, and test section length). In this analysis, the calculated 
variance (a 2 ) and standard deviation (a) values were considered with units, while the 
uncertainty values (u) were considered in percentage terms. 

Mass flowrate (m ) was a function of the experimentally measured and controlled 

volumetric flowrate (V) and the coolant density (p). 

m = Vp (D.4) 

The variance of the mass flowrate can be expressed in the following manner. 


crl =1 


( dm > 

2 

crl + 

r dm' 

[dVj 

V 

Iw 


(D.5) 


The two partial derivatives in Equation D.5 have the following simplified form. 

dm 


dm _ 

~dt~ P 


and — = V 
dp 


(D.6) 


The first partial derivative in Equations D.6 had to be multiplied by the variance 
of the volumetric flowrate. The variance of the single experimental volumetric flowrate 
at OFI, measured using an FT-8 liquid flow transducer manufactured by EG&G/FTI, was 
calculated as the sum of the squares of the individual bias errors and precision errors 
associated with the measurement. The four individual bias errors were: 
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1 . The DAS voltage conversion error was calculated by determining the bit range 
over which the flow transducer was used. The calibration curve for the FT -8 flow 
transducer was: 

gph = 4.224 gph + 0.088135(^1^) ^ (D.7) 

m ^EXP16 

with EXP 16 gain of unity. Since the gain was unity, the millivolt reading on the 
DAS -8 board was equivalent to the millivolt reading on the EXP 16. The 
conversion from millivolts to bits was 2047 bits = 5000 millivolts, accurate to ± 1 
bit. The operational flowrate range associated with this investigation was 10.1 —> 
40.9 gph, which corresponded to a bit range of 28 -» 171 bits. Therefore, for the 
minimum flowrate of 10.1 gph, Ui-low = 1/28 = 0.0357 = 3.57%, and for the 
maximum flowrate of 40.9 gph, uj-High = 1/171 = 0.00585 = 0.585%. 

2. The 20.2 liter jar used to calibrate the FT -8 flowmeter was accurate to 10 mL and, 
therefore, U 2 = 0.010/20.2 = 0.000495 = 0.0495%. 

3. The stopwatch used during the FT -8 flowmeter calibration procedure was 
accurate to 0.001 seconds and the minimum time interval during calibration was 
150.753 (2 minutes, 30.753 seconds) and, therefore, u 3 = 0.001/150.753 = 
0.00000663 = 0.000663%. 

4. The standard error associated with the linear regression calibration curve fit was 
0.300158 gph. For the flowrate range mentioned above (10.1 -» 40.9 gph), U 4 -L 0 W 
= 0.0297 = 2.97% for the lowest flowrate and U 4 -High = 0.00734 = 0.734% for the 
highest flowrate. 
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The only precision error encountered in this investigation with respect to the 
volumetric flowrate came directly from the method of choosing the OFI point visually 
from the plot of AP versus flowrate. As described in Section IV.7.2, the OFI point was 
chosen to be the flowrate (L/hr) associated with the smallest pressure drop (kPa) through 
the test section. Since this value was chosen visually by the operator, some precision 
error was unavoidable. In an attempt to quantify this precision error, a representative set 
of OFI data was fit with two separate quadratic formulas: a fourth order fit (Equation 
D.8) and a fifth order fit (Equation D.9). 

A P = 255.5 - 7.304F + 8.071 x 10' 2 V 2 -3.646 x 10‘ 4 F 3 

(D.OJ 

+ 6.248 xlO'V 4 

A P = 766.2 - 28.36F + 0.4176F 2 - 2.979 x 10' 3 V 2 ^ ^ 

+1.049 x 10 -5 V 4 -1.449 x 10 -8 V 5 

A plot of the representative data set (Setup 5 with nominal conditions 1.5 MW/m heat 
flux, 30°C inlet temperature, 689 kPa exit pressure) is located in Figure D.l along with 
the two quadratic curve fits. 


205 





Coolant Volumetric Flowrate (L/hr) 

* Raw Data.Fourth Order Model Fifth Order Model 

Figure D.l: OFI Plot with Fourth and Fifth Order Quadratic Fits 


Obviously, both quadratic equations fit the data well. By taking the first derivative of 
each quadratic equation with respect to the volumetric flowrate and setting it equal to 
zero, the flowrate value associated with the minimum AP was easily calculated. For the 
fourth and fifth order fits, the respective coolant volumetric flowrates were 83.4 and 84.4 
L/hr for the minimum differential pressures. These values compare favorably to the 
visually chosen OFI volumetric flowrate of 80.0 L/hr, which corresponds to a Gofi value 
of 1022.2 kg/m 2 s (Table B.2). Since the fifth order fit shows the most discrepancy from 
the reported value, its value of coolant volumetric flowrate was used to calculate the 
precision error associated with determining the volumetric flowrate at OFI, therefore, 
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u 


5 


84.4-80.0 

80.0 


= 0.055 = 5.50%. 


By summing the squares of the four bias errors and the single precision error 
associated with the volumetric flowrate measurement at OFI, the error estimate for the 
volumetric flowrate at OFI was found. 


Table D.l: Summary of Volumetric Flowrate Uncertainties at OFI 




Uncertainty 

MORfllS SSfffflRfll 

Uncertainty 

mtimrnsm 




IkiMTSl 



Ui 

Bias 

3.57 % 

1.274 E-3 

0.585 % 

3.422 E-5 

U 2 

Bias 

0.0495 % 

2.450 E-7 

0.0495 % 

2.450 E -7 

mm 

Bias 

0.000663 % 

4.396 E-11 

0.000663 % 

4.396 E-11 

mm 

Bias 

2.97 % 

8.821 E-4 

0.734 % 

5.388 E-5 

mm 

Precision 

5.50 % 

3.025 E-3 

5.50% 

3.025 E-3 

Ellul 

— 

— 

5.181 E-3 

— 

3.113 E -3 


Using this analysis procedure, the error was 7.20 % or 0.727 gph at the lower end 
of the operational flowrate range and 5.58 % or 2.28 gph at the higher end of the 
operational flowrate range. For analysis purposes, the volumetric flowrate error estimate 
was 7.20 % or 2.28 gph, whichever was smaller for the volumetric flowrate value in 
question. Basically, this means that the 7.20 % error value was applied to the lower end 
of the operational flowrate range, while the 2.28 gph error value was applied to the higher 
end of the operational flowrate range. 

The second partial derivative in Equations D.6 had to be multiplied by the 
variance of the coolant density. The precision error associated with the coolant density 
came from the standard deviation of the inlet temperature to the test section during an 
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OFI experiment. The maximum inlet temperature standard deviation encountered over 
all OFI experiments was 1.094°C (from Table B.2), which occurred on an experiment 
that had an average inlet temperature of 52.0°C and an average absolute inlet pressure of 
386.8 kPa. This inlet temperature standard deviation thus created a range of possible 
inlet densities: 986.6 -» 987.7 kg/m 3 . The density value obtained when the average 
values of temperature and pressure were used was 987.2 kg/m . Thus, the representative 
deviation in density from the reported value for all OFI experiments was 0.6 kg/m 3 , 
which corresponded to a u p .i of 0.000608 or 0.0608 %. This value came from the 
standard deviation of the inlet temperature, which depended on the bias error of the inlet 
thermocouple (± 0.1 °C). The lowest inlet temperature encountered during all OFI 
experiments was 25.2°C, which produced a u p .2 of 0.00397 or 0.397 %. By combining 
these two terms as shown in Equation D.10, the final value of u p was 0.00402 or 0.402 %, 
which means that the standard deviation of the inlet temperature contributes very little to 
the error associated with the coolant density at the test section inlet. 

“, = lk-i 2 + v /= 0 - 00402 (D.10) 

Referring to Equation D.5, all values were known at this point such that the 
variance of the coolant mass flowrate could be determined. Since the value of 
importance to this investigation was coolant mass flux, the variance of the test section 
cross sectional area had to be determined for the error analysis to be complete, as shown 
in Equation D.2. The cross sectional area or flow area of the test section is a function of 
the inner diameter of the outer heater (outer wall of the annulus) and the outer diameter of 
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the inner heater (inner wall of the annulus). 


4,-f (pL-ol) CD-ii) 

The variance of the test section cross sectional area can be expressed in the following 
manner. 



/ 8A„ V 


dD. 




out / 






(D.12) 


where: 


dA 


dD, 


= *D. 


out 


and = (D.13) 

az>, 2 


The outer diameter (inner diameter of the outer heater) was measured using both 
precision ground pins and a Starrett dial caliper. The tolerance of the caliper was known 
to be <r n = 0.001 inch = 0.0254 mm, which created a maximum bias uncertainty 

estimate of u n = 0.001/0.304 = 0.00329 = 0.329 % for the smallest diameter outer 

heater used (Outer Heaters I, J, and K). This tolerance was used instead of the tolerance 
of the pins to provide a more conservative estimate of the error associated with the cross 
sectional area of the test section. The inner diameter (outer diameter of the inner heater) 
was also measured using a Starrett dial caliper (cr D . n — 0.001 inch = 0.0254 mm). All 

inner heaters were constructed of tubing that had been manufactured with a nominal outer 
diameter of 0.250 inch; thus, all inner heaters had roughly the same diameter. During 
construction, all inner heaters were found to have diameters slightly larger than 0.250 
inch, but 0.250 inch was utilized in this error analysis to maintain a conservative 
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estimate. The bias uncertainty estimate of the inner diameter was found to be u D = 
0.001/0.250 = 0.00400 = 0.400 %. 


Finally, combining Equations D.2 through D.13 yielded the following estimate of 
the coolant mass flux variance. 


&G = 


J-] \p 2 <r 2 y + vW p \+ 


f • V 

-m 


A 2 

\ A cs J 


* 

^ out 
V 2 J 




( 


—D- 
2 m 


'D:„ 


(D.14) 


By applying this formula along with the appropriate error estimates from Table D.2 to the 
set of experiments performed in this investigation, the average uncertainty in mass flux 
was calculated to be 7.46 %. Note that the a values for coolant volumetric flowrate in 
Table D.2 are given in terms of both gph and percentage since the error estimate of the 
flowrate in question depends on the magnitude of the flowrate. The mass flux error 
analysis results for the data with heat flux ratio of unity are listed in Table D.3 


Table D.2: Summary of Coolant Mass Flux Error Sources at OFI 


Gy 

Up 

cr D 


120% ox 2.2% gph 

0.402 % 

0.0254 mm 

0.0254 mm 
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Table D.3: OFI Coolant Mass Flux Error Analysis Results 


Sigma(Vol. Flow): 

2.40E-06 

(m A 3/sec) 

or 

7.20 

(%> 

Sigma(Density): 

0.402 

(%> 

Sigma(D_out): 

2.54E-05 

(m) 

Sigma(D in): 

2.54E-05 

(m) 


Setup 

Run 

Inlet Density 
@OFI 
(kg/m A 3) 

Inlet Density 
sigma 
(kg/m A 3) 

Vol Flowrate 
@OFI 
(m A 3/sec) 

Vol Flowrate 
sigma 
(m A 3/sec) 

Mass 

Flowrate 

(kg/s) 

Inner 

Diameter 

(mm) 

Outer 

Diameter 

(mm) 

A_cs 

(m A 2) 

G(ofi) 

(kg/m A 2 s) 

SigmaG 

(kg/m A 2 s) 

Uncertainty_G 

(%) 

5 

6 

991 

3.98 

1.18E-05 

8.50E-07 

0.0117 

6.401 

8.280 

2.167E-05 

539.66 

40.30 

7.47 

5 

7 

987 

3.97 

1.92E-05 

1.38E-06 

0.0190 

6.401 

8.280 

2.167E-05 

874.78 

65.32 

7.47 

5 

8 

988 

3.97 

1.07E-05 

7.68E-07 

0.0105 

6.401 

8.280 

2.167E-05 

486.07 

36.30 

7.47 

5 

9 

987 

3.97 

1.34E-05 

9.64E-07 

0.0132 

6.401 

8.280 

•2.167E-05 

609.72 

45.53 

7.47 

5 

10 

980 

3.94 

1.64E-05 

1.18E-06 

0.0161 

6.401 

8.280 

2.167E-05 

741.84 

55.40 

7.47 

5 

11 

998 

4.01 

1.14E-05 

8.18E-07 

0.0113 

6.401 

8.280 

2.167E-05 

522.93 

39.05 

7.47 

5 

12 

988 

3.97 

1.76E-05 

1.27E-06 

0.0174 

6.401 

8.280 

2.167E-05 

801.03 

59.82 

7.47 

5 

13 

989 

3.97 

- 

- 

- 

6.401 

8.280 

2.167E-05 


_ 


5 

14 

987 

3.97 

1.62E-05 

1.17E-06 

0.0160 

6.401 

8.280 

2.167E-05 

736.89 

55.03 

7.47 

5 

15 

980 

3.94 

1.82E-05 

1.31E-06 

0.0178 

6.401 

8.280 

2.167E-05 

821.50 

61.34 

7.47 

5 

16 

998 

4.01 

2.94E-05 

2.11E-06 

0.0293 

6.401 

8.280 

2.167E-05 

1350.83 

100.87 

7.47 

5 

20 

990 

3.98 

1.90E-05 

1.37E-06 

0.0188 

6.401 

8.280 

2.167E-05 

869.17 

64.90 

7.47 

5 

21 

987 

3.97 

3.08E-05 

2.22E-06 

0.0304 

6.401 

8.280 

2.167E-05 

1402.13 

104.70 

7.47 

5 

22 

988 

3.97 

1.71E-05 

1.23E-06 

0.0169 

6.401 

8.280 

2.167E-05 

779.70 

58.22 

7.47 

5 

23 

988 

3.97 

1.99E-05 

1.43E-06 

0.0197 

6.401 

8.280 

2.167E-05 

907.41 

67.76 

7.47 

5 

24 

980 

3.97 

- 

- 

... 

6.401 

8.280 

2.167E-05 

.. 

_ 


5 

25 

997 

4.01 

2.22E-05 

1.60E-06 

0.0222 

6.401 

8.280 

2.167E-05 

1022.25 

76.34 

7.47 

5 

26 

991 

3.98 

3.23E-05 

2.33E-06 

0.0320 

6.401 

8.280 

2.167E-05 

1475.87 

110.21 

7.47 

5 

27 

989 

3.98 

2.41E-05 

1.74E-06 

0.0239 

6.401 

8.280 

2.167E-05 

1101.23 

82.23 

7.47 

5 

28 

988 

3.97 

2.74E-05 

1.97E-06 

0.0271 

6.401 

8.280 

2.167E-05 

1249.29 

93.29 

7.47 

5 

29 

982 

3.95 

3.04E-05 

2.19E-06 

0.0299 

6.401 

8.280 

2.167E-05 

1377.32 

102.85 

7.47 

5 

30 

997 

4.01 

2.28E-05 

1.64E-06 

0.0227 

6.401 

8.280 

2.167E-05 

1048.31 

78.28 

7.47 

6 

32 

984 

3.96 

2.39E-05 

1.72E-06 

0.0235 

6.401 

8.280 

2.167E-05 

1085.26 

81.04 

7.47 

6 

33 

987 

3.97 

3.13E-05 

2.26E-06 

0.0309 

6.401 

8.280 

2.167E-05 

1426.83 

106.55 

7.47 













Table D.3 (continued) 


Setup 

Run 

Inlet Density 
@OFI 
(kg/m A 3) 

Inlet Density 
sigma 
(kg/m A 3) 

Vol Flowrate 
@OFI 
(m A 3/sec) 

Vol Flowrate 
sigma 
(m A 3/sec) 

Mass 

Flowrate 

(kg/s) 

Inner 

Diameter 

(mm) 

Outer 

Diameter 

(mm) 

A_cs 

(m A 2) 

G(ofi) 

(kg/m A 2 s) 

Sigma_G 

(kg/m A 2 s) 

Uncertainty_G 

<%) 

6 

34 

984 

3.95 

3.60E-05 

2.40E-06 

0.0354 

6.401 

8.280 

2.167E-05 

1633.30 

113.47 

6.95 

6 

35 

987 

3.97 

3.13E-05 

2.26E-06 

0.0309 

6.401 

8.280 

2.167E-05 

1298.76 

106.55 

8.20 

6 

36 

990 

3.98 

1.64E-05 

1.18E-06 

0.0162 

6.401 

8.280 

2.167E-05 

746.70 

55.76 

7.47 

6 

37 

990 

3.98 

3.64E-05 

2.40E-06 

0.0361 

6.401 

8.280 

2.167E-05 

1663.18 

114.33 

6.87 

6 

39 

991 

3.98 

3.68E-05 

2.40E-06 

0.0365 

6.401 

8.280 

2.167E-05 

1683.66 

114.54 

6.80 

8 

1 

996 

4.01 

1.31E-05 

9.42E-07 

0.0130 

6.401 

7.925 

1.715E-05 

759.85 

57.71 

7.59 

8 

2 

990 

3.98 

1.41E-05 

1.02E-06 

0.0140 

6.401 

7.925 

1.715E-05 

814.31 

61.85 

7.59 

9 

3 

982 

3.95 

1.54E-05 

1.11E-06 

0.0151 

6.401 

7.925 

1.715E-05 

879.55 

66.80 

7.59 

9 

4 

982 

3.95 

1.66E-05 

1.20E-06 

0.0163 

6.401 

7.925 

1.715E-05 

950.36 

72.18 

7.59 

9 

5 

988 

3.97 

1.60E-05 

1.15E-06 

0.0158 

6.401 

7.925 

1.715E-05 

919.51 

69.84 

7.59 

9 

6 

996 

4.00 

1.35E-05 

9.74E-07 

0.0135 

6.401 

7.925 

1.715E-05 

785.50 

59.66 

7.59 

9 

7 

990 

3.98 

1.41E-05 

1.02E-06 

0.0140 

6.401 

7.925 

1.715E-05 

815.52 

61.94 

7.59 

9 

8 

988 

3.97 

1.19E-05 

8.54E-07 

0.0117 

6.401 

7.925 

1.715E-05 

683.34 

51.90 

7.59 

9 

9 

997 

4.01 

1.76E-05 

1.26E-06 

0.0175 

6.401 

7.925 

1.715E-05 

1019.99 

77.47 

7.59 

9 

10 

989 

3.97 

2.17E-05 

1.56E-06 

0.0215 

6.401 

7.925 

1.715E-05 

1251.79 

95.07 

7.59 

9 

11 

981 

3.94 

2.62E-05 

1.89E-06 

0.0257 

6.401 

7.925 

1.715E-05 

1499.75 

113.90 

7.59 

9 

12 

988 

3.97 

1.79E-05 

1.29E-06 

0.0177 


7.925 

1.715E-05 

1030.09 

78.23 

7.59 

9 

13 

996 

4.01 

2.11E-05 

1.52E-06 

0.0210 


7.925 

1.715E-05 

1222.95 

92.88 

7.59 

9 

14 

989 

3.97 

2.44E-05 

1.76E-06 

0.0241 


7.925 

1.715E-05 

1407.01 

106.86 

7.59 

9 

15 

981 

3.94 

2.94E-05 

2.12E-06 

0.0289 


7.925 

1.715E-05 

1682.18 

127.76 

7.59 

9 

16 

989 

3.97 

2.19E-05 

1.57E-06 

0.0216 


7.925 

1.715E-05 

1259.97 

95.69 

7.59 

11 

1 

995 

4.00 

3.08E-05 

2.22E-06 

0.0307 


7.722 

1.465E-05 

2086.20 

161.45 

7.74 

11 

2 

988 

3.97 

3.35E-05 

2.40E-06 

0.0331 


7.722 

1.465E-05 

2251.86 

173.35 

7.70 

12 

4 

995 

4.00 

4.08E-05 

2.40E-06 

0.0406 


7.722 

1.465E-05 

2762.58 

179.98 

6.51 

12 

5 

989 

3.98 

4.31E-05 

2.40E-06 

0.0426 


7.722 

1.465E-05 

2897.07 

180.70 

6.24 

12 

6 

995 

4.00 

2.97E-05 

2.14E-06 

0.0295 


7.722 

1.465E-05 

2008.69 

155.45 

7.74 






D.2 Heat Flux Error Analysis 


Heat flux can be calculated as a function of the power applied to the test section 
and the test section heated area. 



Qjn 

~ 4 


(D.15) 


The variance of the surface heat flux can be expressed by the following error propagation 
formula. 



f 8 Q' a 

2 

<ri + 

(dQ" > 

2 

<J~ n + 

f 8Q") 


Qin 

^Qout) 

a out 

l dA s) 


'A, 


(D.16) 


where 

dQ^ = _d&_ = _l ag" ( O L +0j = Own 

dQ in dQ out 4 54 A) A] 


(D.17) 


The variance of the inner and outer heater power levels at OFI was calculated as the sum 
of the squares of the individual bias errors and precision errors associated with the 
measurement. The first 61 OFI experiments utilized an AC wattmeter manufactured by 
Ohio Semitronics with an external Omega AC current transducer (Components 15 of 
Table 3.3) to measure the power supplied to the outer heater, while the inner heater 
power was measured using an Ohio Semitronics DC wattmeter (Component 16 of Table 
3.3). The final five OFI experiments measured the inner and outer heater power levels 
using redundant sensors. Since only the total power level reading supplied by the Ohio 
Semitronics DC wattmeter (Component 16 of Table 3.3) was used for analysis of these 
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last five OFI experiments, only the tolerances of this wattmeter wsd considered in this 
error analysis. 

During the first 61 OFI experiments, there were three individual bias errors and 
no precision errors that contributed to the variance of the outer heater power: 

1. The manufacturer stated accuracy of the Ohio Semitronics AC wattmeter was 
± 0.2 % of the power reading, thus Ui = 0.002 = 0.2 %. 

2 . The manufacturer stated accuracy of the external Omega AC current transducer 
was ±1.5 %, thus U 2 = 0.015 = 1.5 %. 

3. The DAS voltage conversion error was calculated by determining the bit range 
over which the wattmeter was used. The calibration curve for the wattmeter was: 

Power (W) = -2.31W ± 6.4097(w F^J —(D. 18) 

mV EXP16 

with EXP 16 gain of unity. Since the gain was unity, the millivolt reading on the 
DAS -8 board was the same as on the EXP 16 board. The conversion from 
millivolts to bits was 2047 bits = 5000 mV, accurate to ±1 bit. The outer heater 
power over the first 61 OFI experiments ranged from a minimum of 1230 W to a 
maximum of 9311 W, which converted to a range of 80 —> 596 bits. Therefore, 
for the 1230 W power level, U 3 -L 0 W = 1/80 = 0.0125 = 1.25 % and for the 9311 W 
power level, U 3 -Hi g h = 1/596 = 0.00168 = 0.168 %. 

By summing the squares of the three individual bias errors associated with the 
outer heater power measurement at OFI, the error estimate for the outer heater power for 
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the first 61 OFI experiments was found. 


Table D.4: Summary of Outer Heater Power Uncertainties at OFI (first 61 experiments) 






Uncertainty 

BSSRI!S!&I 


| 



(@9311 W) 


Ui 

Bias 

1 

4.00 E -6 

0.2 % 

4.00 E -6 

U 2 

Bias 


2.25 E-4 

1.5 % 

2.25 E -4 

U3 

Bias 

1.25% 

1.56 E -4 

0.168 % 

2.82 E-6 

Sum 

— 

— 

3.85 E-4 

— 

2.32 E-4 


Using this analysis procedure, the error was 1.96 % or 24.1 W at the lower end of 
the operational power range and 1.52 % or 141.5 W at the higher end of the operational 
power range. For analysis purposes, the outer heater power error estimate was 1.96 % or 
141.5 W, whichever was smaller for the outer heater power value in question. Basically, 
this means that the 1.96 % error value was applied to the lower end of the operational 
power range, while the 141.5 W error value was applied to the higher end of the 
operational power range. 

During the first 61 OFI experiments, there were two individual bias errors that 
contributed to the variance of the inner heater power measurement: 

1. The manufacturer stated accuracy of the Ohio Semitronics DC wattmeter was 
±1.0 % of the power reading, thus ui = 0.01 = 1.0 %. 

2. The DAS voltage conversion error was calculated by determining the bit range 
over which the wattmeter was used. The calibration curve for this DC wattmeter 
was the same as that of the AC wattmeter used on the outer heater: 
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Power (W) = -2.31 W + 6.4097(iwr £XPI6 ) —(D.19) 

mV EXP!6 

with EXP 16 gain of unity. Since the gain was unity, the millivolt reading on the 
DAS-8 board was the same as on the EXP 16 board. The conversion from 
millivolts to bits was 2047 bits = 5000 mV, accurate to ±1 bit. The inner heater 
power over the first 61 OFI experiments ranged from a minimum of 3134 W to a 
maximum of 7637 W, which converted to a range of 201 —> 489 bits. Therefore, 
for the 3134 W power level, U 2 -L 0 W = 1/201 = 0.00498 = 0.498 % and for the 7637 
W power level, U 2 -High - 1/489 = 0.00204 = 0.204 %. 


By summing the squares of the two individual bias errors associated with the 
inner heater power measurement at OFI, the error estimate for the inner heater power for 
the first 61 OFI experiments was found. 


Table D.5: Summary of Inner Heater Power Uncertainties at OFI (first 61 experiments) 



Error Type 

Uncertainty 

Uncertainty 2 

Uncertainty 

Uncertainty 2 



IfellKliW 



m 7637 W) 

Ui 

Bias 

1.0% 

1.000 E-4 


1.000 E ^4 

U 2 

Bias 

0.498 % 

2.480 E -5 

0.204 % 

4.162 E-6 

Sum 

— 

— 

1.248 E -4 

— 

1.041 E -4 


Using this analysis procedure, the error was 1.12 % or 35.1 W at the lower end of 
the operational power range and 1.02 % or 77.9 W at the higher end of the operational 
power range. For analysis purposes, the outer heater power error estimate was 1.12 % or 
77.9 W, whichever was smaller for the inner heater power value in question. Basically, 
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this means that the 1.12 % error value was applied to the lower end of the operational 
power range, while the 77.9 W error value was applied to the higher end of the 
operational power range. 

For the final five OFI experiments, the power supplied to the test section was 
determined using an energy balance and the total amount of power supplied by the DC 
power supply as described in Section IV.7.2. Errors in the energy balance heat loss 
estimation were discussed independently in Section IV.7.2, however, the error associated 
with the measurement of the total power output from the power supply had to be 
quantified. This power level was measured using the Ohio Semitronics DC wattmeter 
(Component 16 of Table 3.3) mentioned earlier. During the final five OFI experiments, 
there were two individual bias errors that contributed to the uncertainty of the total power 
reading: 

1. The manufacturer stated accuracy of the Ohio Semitronics DC wattmeter was 
±1.0 % of the power reading, thus ui = 0.01 = 1.0 %. 

2. The DAS voltage conversion error was calculated by determining the bit range 
over which the wattmeter was used. The calibration curve for this DC wattmeter 
was the same as stated above: 

Power (W) = -2.31 W + 6.4097(/mF £ * p16 )——— (D.20) 

mV EXP16 

with EXP 16 gain of unity. Since the gain was unity, the millivolt reading on the 
DAS-8 board was the same as on the EXP 16 board. The conversion from 
millivolts to bits was 2047 bits = 5000 mV, accurate to ±1 bit. The total power 
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level over the final five OFI experiments ranged from a minimum of 14041 W to 
a maximum of 20160 W, which converted to a range of 897 -> 1288 bits. 
Therefore, for the 14041 W power level, U 2 -L 0 W = 1/897 = 0.00111 =0.111 % and 
for the 20160 W power level, U 2 -High = 1/1288 = 0.000776 = 0.0776 %. 


By summing the squares of the two individual bias errors associated with the total 
power measurement at OFI, the error estimate for the inner heater power for the last five 
OFI experiments was found. 


Table D.6: •Summary of Total Power Level Uncertainties at OFI (last five experiments) 



ES9 

BBSBI 

Uncertainty 2 

Uncertainty 

Uncertainty 2 



m 14041 W) 

((S) y 14041 W) 

m 20160 W) 

m 20160 W) 

Ui 

Bias 

1.0% 

1.000 E-4 

1.0% 

1.000 E-4 

u 2 

Bias 

0.111 % 

1.232 E-6 

0.0776 % 

6.022 E -7 

Sum 

— 

-- 

1.012 E-4 

— 

1.006 E -4 


Using this analysis procedure, the error was 1.01 % or 141.8 W at the lower end 
of the operational power range and 1.00 % or 201.6 W at the higher end of the 
operational power range. For analysis purposes, the total power error estimate was 1.01 
% or 201.6 W, whichever was smaller for the total power value in question. Basically, 
this means that the 1.01 % error value was applied to the lower end of the operational 
power range, while the 201.6 W error value was applied to the higher end of the 
operational power range. 

The test section heated surface area was a function of the inner diameter of the 
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annulus, outer diameter of the annulus, inner heater heated length, and outer heater heated 
length. 

A s = heated, oufiout ^heated .in ) (D-21) 

The variance of the test section heated surface area can be expressed in the following 
manner. 


< = 


< dA. > a 


\ dD ou,J 




f dA s V2 


a k + 


8A . 


V 


\8L h ' 0Ut j 


L h. 


f ^ 

SA s 

SL h,in J 


O* 

L h. 


(D.22) 


where 


dDou, 


— 7lL 


'h,out 


and 


dA 


dD, 


~ = &h,in 


(D.23) 


dA . 


dL 


= AD n , 


and 


dA. 


( h,out 


dL 


= AD ; . 


'h,in 


As stated earlier, the relevant annulus diameters were measured using a Starrett dial 
caliper, which had a known tolerance of cr = 0.001 inch. The heated lengths of all test 
sections were measured using Fowler vernier calipers, also with a known tolerance of a = 
0.001 inch. 

Finally, combining Equations D.16 through D.23 yielded the following estimate 
of the surface heat flux variance. 



(D.24) 
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By applying this formula along with the appropriate error estimates from Table 
D.7 to the set of experiments performed in this investigation, the average uncertainty in 
heat flux was calculated to be 1.25 %. Note that the a values for all power error 
estimates in Table D.7 are given in terms of both Watts and percentage since the error 
estimate of the power level in question depends on the magnitude of the power 
measurement. 


Table D.7: Summary of Heat Flux Uncertainties at OFI 


a Q 

Min 

(first 61) 

a Q 

Mout 

(first 61) 

a Q,^ 

(last 5) 

a D m , 

cr D 

'-'in 

Cr 

**h,out 

C>l 

L h,in 

1.96% or 
141.5 W 

1.12% or 
77.9 W 

1.01 % or 
201.6 W 

0.0254 mm 

0.0254 mm 

0.0254 mm 

0.0254 mm 


The outline below shows a sample error estimate for an OFI experiment that had 
the following experimentally measured parameters: D ou t = 7.925 mm, D in = 6.401 mm, 
Lh, out =180 mm, Lh, m = 180 mm, Qj n = 5424 W, and Q out = 6651 W, with an 
environmental power loss of 1.4 %. The heat flux error analysis results for the data with 
heat flux ratio of unity are listed in Table D.8 
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Measurement 

Uncertainty 

Influence Coefficient, f a 9 "'j 

(Uncertainty) 2 

Basis 

Inner Heater 
Power 

1.96%= 106.3 W 

■■■■■■ 

1.72 E -4 (MW/m 2 ) 2 

Resolution 

■ 

Outer Heater 
Power 

1.12 % = 74.5 W 

mmmmmm 

8.45 E-5 (MW/m 2 ) 2 

Resolution 


Outer Diameter 

0.0254 mm 

Q total ^ out 

_| 

Measurement 

x( D ou, L ou, + D in L in )' 

= -102.6 (MW/m 3 ) 

Inner Diameter 

0.0254 mm 

~~ Q total ^ in 

6.79 E -6 (MW/m 2 ) 2 

Measurement 

*( D ou, L ou, + D m L in f 

= -102.6 (MW/m 3 ) 

Outer L h 

0.0254 mm 

Q total ^ out 

1.32 E -8 (MW/m 2 ) 2 

Measurement 

n{ D ou,Lou, + D in Lin Y 

= -4.52 (MW/m 3 ) 

Inner L h 

0.0254 mm 

Q total ^ in 

8.60 E-9 (MW/m 2 ) 2 

Measurement 





2.70 E -4 (MW/m 2 ) 2 


BfSfjPHPPM 

■i 
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Table D.8: OFI Surface Heat Flux Error Analysis Results 


Sigma(L_h, out): 

2.54E-05 

(m) 

Sigma(L_h, in): 

2.54E-05 

(m) 

Sigma(D_out): 

2.54E-05 

(m) 

Sigma(D_in): 

2.54E-05 

("») 


Setup Run Inner Q_inner Outer Q_outer Power Loss Total Q_total Inner Inner Outer Outer A_s Q"(ofi) Sigma_Q" Uncertainty Q" 

Ppwer sigma Power sigma Environment Power sigma Diameter L_heated Diameter L_heated 




(W) 

(W) 

(W) 

(W) 

_ (%> 

_ 

(W) 

(mm) 

(mm) 

(mm) 

(mm) 

(m A 2) 

(MW/m A 2) 

(MW/m A 2) 

(%) 

5 

6 

3338 

65 

4010 

45 

14.4 

6290 

- 

6.401 

190 

8.280 

180 

8.503E-03 

0.7356 

9.51E-03 

1.29 

5 

7 

3345 

66 

3995 

45 

13.2 

6371 

- 

6.401 

190 

8.280 

180 

8.503E-03 

0.7445 

9.51E-03 

1.28 

5 

8 

3331 

65 

3961 

44 

14.0 

6271 

-- 

6.401 

190 

8.280 

180 

8.503E-03 

0.7328 

9.46E-03 

1.29 

5 

9 

3333 

65 

3980 

45 

12.8 

6377 

- 

6.401 

190 

8.280 

180 

8.503E-03 

0.7452 

9.48E-03 

1.27 

5 

10 

3335 

65 

4015 

45 

16.0 

6174 

- 

6.401 

190 

8.280 

180 

8.503E-03 

0.7220 

9.50E-03 

1.32 

5 

11 

3787 

74 

4507 

50 

16.1 

6958 

- 

6.401 

190 

8.280 

180 

8.503E-03 

0.8130 

1.07E-02 

1.32 

5 

12 

3801 

74 

4518 

51 

27.3 

6048 

- 

6.401 

190 

8.280 

180 

8.503E-03 

0.7074 

1.07E-02 

1.52 

5 

13 

- 

- 

— 

~ 

— 

-- 

- 

6.401 

190 

8.280 

180 

8.503E-03 

_ 

__ 


5 

14 

3781 

74 

4513 

51 

15.6 

7000 

- 

6.401 

190 

8.280 

180 

8.503E-03 

0.8182 

1.07E-02 

1.31 

5 

15 

3775 

74 

4471 

50 

15.7 

6951 

- 

6.401 

190 

8.280 

180 

8.503E-03 

0.8121 

1.07E-02 

1.32 

5 

16 

7637 

142 

9311 

78 

6.9 

15778 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.8441 

1.95E-02 

1.06 

5 

20 

4752 

93 

5850 

66 

3.3 

10252 

-- 

6.401 

190 

8.280 

180 

8.503E-03 

1.1980 

1.37E-02 

1.14 

5 

21 

4778 

94 

5821 

65 

3.4 

10239 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.1974 

1.37E-02 

1.15 

5 

22 

4775 

94 

5923 

66 

5.0 

10163 

-- 

6.401 

190 

8.280 

180 

8.503E-03 

1.1879 

1.38E-02 

1.16 

5 

23 

4763 

93 

5838 

65 

2.4 

10346 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.2089 

1.37E-02 

1.14 

5 

24 

- 

-- 

- 

-- 

— 

-- 

- 

6.401 

190 

8.280 

180 

8.503E-03 

_ 

__ 


5 

25 

5886 

115 

7070 

78 

0.6 

12878 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.4815 

1.68E-02 

1.13 

5 

26 

5697 

112 

7051 

78 

0.7 

12659 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.4793 

1.64E-02 

1.11 

5 

27 

5674 

111 

6977 

78 

1.7 

12436 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.4536 

1.64E-02 

1.13 

5 

28 

5687 

111 

6934 

78 

1.9 

12381 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.4478 

1.64E-02 

1.13 

5 

29 

5716 

112 

6971 

78 

1.4 

12509 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.4630 

1.65E-02 

1.13 

5 

30 

6677 

131 

8092 

78 

12.9 

12864 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.5039 

1.83E-02 

1.22 

6 

32 

4676 

92 

5943 

67 

2.3 

10375 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.2128 

1.37E-02 

1.13 

6 

33 

6635 

130 

8203 

78 

1.8 

14571 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.7037 

1.83E-02 

1.08 



Table D.8 (continued) 


Outer 

Power 

(W) 

Qouter 

sigma 

(W) 

Power Loss 
Environment 
(%) 

Total 

Power 

(W) 

Q_total 

sigma 

(W) 

Inner 

Diameter 

(mm) 

Inner 

L_heated 

(mm) 

Outer 

Diameter 

(mm) 

Outer 

L_heated 

(mm) 

A_s 

(m A 2) 

Q"(ofi) 

(MW/m A 2) 

SigmaQ" 

(MW/m A 2) 

UncertaintyQ 

(%) 

8131 

78 

3.1 

14320 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.6746 

1.83E-02 

1.09 

8175 

78 

3.5 

14284 

-- 

6.401 

190 

8.280 

180 

8.503E-03 

1.6689 

1.83E-02 

1.10 

4805 

54 

3.0 

8387 

-- 

6.401 

190 

8.280 

180 

8.503E-03 

0.9803 

1.12E-02 

1.14 

8180 

78 

1.9 

14535 

- 

6.401 

190 

8.280 

180 

8.503E-03 

1.6987 

1.83E-02 

1.08 

8554 

78 

1.0 

15534 

-- 

6.401 

190 

8.280 

180 

8.503E-03 

1.8165 

1.94E-02 

1.07 

3803 

43 

2.9 

6808 

- 

6.401 

180 

7.925 

180 

8.101E-03 

0.8400 

9.61E-03 

1.14 

3770 

42 

2.9 

6767 

- 

6.401 

180 

7.925 

180 

8.101E-03 

0.8352 

9.56E-03 

1.15 

3936 

44 

7.7 

6525 

- 

6.401 

180 

7.925 

180 

8.101E-03 

0.8057 

9.55E-03 

1.19 

4437 

50 

6.0 

7552 

-- 

6.401 

180 

7.925 

180 

8.101E-03 

0.9324 

1.09E-02 

1.17 

4448 

50 

4.7 

7662 

“ 

6.401 

180 

7.925 

180 

8.101E-03 

0.9462 

1.09E-02 

1.15 

4421 

50 

2.9 

7854 

” 

6.401 

180 

7.925 

180 

8.101E-03 

0.9692 

1.10E-02 

1.14 

4440 

50 

4.0 

7721 

“ 

6.401 

180 

7.925 

180 

8.101E-03 

0.9533 

1.09E-02 

1.15 

3918 

44 

2.9 

6878 

-- 

6.401 

180 

7.925 

180 

8.101E-03 

0.8491 

9.62E-03 

1.13 

5617 

63 

2.6 

9865 

-- 

6.401 

180 

7.925 

180 

8.101E-03 

1.2181 

1.37E-02 

1.13 

5653 

63 

2.7 

9874 

-- 

6.401 

180 

7.925 

180 

8.101E-03 

1.2185 

1.37E-02 

1.13 

5604 

63 

3.5 

9788 

~ 

6.401 

180 

7.925 

180 

8.101E-03 

1.2088 

1.38E-02 

1.14 

5542 

62 

3.6 

9652 

-- 

6.401 

180 

7.925 

180 

8.101E-03 

1.1909 

1.36E-02 

1.14 

6658 

75 

2.1 

11832 

- 

6.401 

180 

7.925 

180 

8.101E-03 

1.4607 

1.65E-02 

1.13 

6612 

74 

1.5 

11849 

-- 

6.401 

180 

7.925 

180 

8.101E-03 

1.4620 

1.64E-02 

1.12 

6582 

74 

1.8 

11768 

-- 

6.401 

180 

7.925 

180 

8.101E-03 

1.4527 

1.63E-02 

1.12 

6651 

74 

1.4 

11906 

- 

6.401 

180 

7.925 

180 

8.101E-03 

1.4691 

1.64E-02 

1.12 

- 

. -- 

- 

14078 

141 

6.401 

180 

7.722 

185 

8.107E-03 

1.7302 

1.79E-02 

1.04 

- 

-- 

~ 

14041 

140 

6.401 

180 

7.722 

185 

8.107E-03 

1.7519 

1.79E-02 

1.02 

-- 

-- 

- 

19528 

195 

6.401 

190 

7.722 

185 

8.308E-03 

2.3419 

2.43E-02 

1.04 

- 

-- 

- 

20160 

202 

6.401 

190 

7.722 

185 

8.308E-03 

2.4177 

2.50E-02 

1.04 

-- 

— 

— 

16614 

166 

6.401 

190 

7.722 

185 

8.308E-03 

1.9925 

2.06E-02 

1.04 



APPENDIX E 


CHF ERROR ANALYSIS 


This appendix presents an error analysis of the experimentally obtained CHF heat 
flux and CHF mass flux values for all experiments performed in this investigation. 


E.l Mass Flux Error Analysis 

Mass flux can be calculated as a function of the coolant mass flowrate and the test 
section cross-sectional area. 


G 



(E.l) 


Likewise, the variance of the coolant mass flux can be expressed by the following error 
propagation formula. 



(dG' 

2 

+ 

f 5G 1 

\drhj 



l dA cJ 



(E.2) 


where 


3G_ J_ 
dm A cs 


and 


dG _ m 

Hs 4 


(E.3) 
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Since the variance of the coolant mass flowrate and the test section cross-sectional area 
were not directly known, these quantities had to be expanded into known components 
until all variances were in terms of experimentally measured variables (i.e., volumetric 
flowrate, power, test section diameter, and test section length). In this analysis, the 
calculated variance (a 2 ) and standard deviation (a) values were considered with units, 
while uncertainty values (u) were considered in percentage terms. 

Mass flowrate ( m ) was a function of the experimentally measured and controlled 

volumetric flowrate ( V ) and the coolant density (p). 

m = Vp (E.4) 


The variance of the mass flowrate can be expressed in the following manner. 

(E.5) 




f dm V 2 

r dm'' 

—- \<?v + 

{dVJ v 

l d Pj 




The two partial derivatives in Equation E.5 have the following form. 

dm 


dm 

■¥ =p 


and 


dp 


= V 


(E.6) 


The first partial derivative in Equations E.6 had to be multiplied by the variance 
of the volumetric flowrate. The variance of the volumetric flowrate at CHF, measured 
using an FTO-3 liquid flow transducer manufactured by EG&G/FTI, was calculated as 
the sum of the squares of the individual bias errors and precision errors associated with 
the measurement. The four individual bias errors were: 

1. The DAS voltage conversion error was calculated by determining the bit range 
over which the flow transducer was used. The calibration curve for the FTO-3 
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flow transducer was: 


gph = 0.432929 gph + 0.00495(? W F £ffl6 ) 

EXP 16 

with EXP 16 gain of unity. Since the gain was unity, the millivolt reading on the 
DAS -8 board was equivalent to the millivolt reading on the EXP16. The 
conversion from millivolts to bits was 2047 bits = 5000 millivolts, accurate to + 1 
bit. The operational flowrate range associated with this investigation was 1.3 
5.5 gph, which corresponded to a bit range of 72 —> 420 bits. Therefore, for the 
minimum flowrate of 1.3 gph, Ui.low = 1/72 = 0.0139 = 1.39%, and for the 
maximum flowrate of 5.5 gph, ui-High- 1/420 = 0.00238 = 0.238%. 

2. The 1000 mL graduated cylinder used to calibrate the FTO-3 flowmeter was 
accurate to 5 mL and, therefore, U 2 = 5/1000 = 0.005 = 0.5%. 

3. The stopwatch used during the FTO-3 flowmeter calibration procedure was 
accurate to 0.001 seconds and the minimum time interval during calibration was 
35 seconds and, therefore, U 3 = 0.001/35 = 0.0000286 = 0.00286%. 

4. The standard error associated with the linear regression calibration curve fit was 
0.116035 gph. For the flowrate range mentioned above (1.3 —> 5.5 gph), U 4 .L 0 W = 
0.0893 = 8.93% for the lowest flowrate and u 4 -Hi g h = 0.0211 = 2.11% for the 
highest flowrate. 

The only precision error encountered in this investigation with respect to the 
volumetric flowrate came directly from the standard deviation of the flowrate during a 
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CHF experiment. The maximum standard deviation of the volumetric flowrate was 0.83 
L/hr (15.66 kg/m 2 s), which occurred during an experiment that maintained the flowrate 
constant at 11.04 L/hr (207.74 kg/m 2 s), therefore, U 5 = 0.83/11.04 = 0.0752 = 7.52 %. 

By summing the squares of the four bias errors and the single precision error 
associated with the volumetric flowrate measurement at CHF, the error estimate for the 
volumetric flowrate at CHF was found. 


Table E.l: Summary of Volumetric Flowrate Uncertainties at CHF 



BiBIlBiBi 


E3S5S351 

Uncertainty 




wmmm 



(@5.5 gph) 

Ui 

Bias 

1 . 39 % 


0.238 % 

5.664 E-6 

U 2 

Bias 

0.5 % 

EEHU 

0.5 % 

2.500 E -5 

U 3 

Bias 

0.00286 % 

8.180 E-10 

0.00286 % 

8.180 E-10 

U 4 

Bias 

| 

7.974 E -3 


4.452 E-4 

u 5 

Precision 

| | 

5.655 E-3 



Sum 


— 

1.385 E -2 

— 



Using this analysis procedure, the error was 11.8 % or 0.153 gph at the lower end 
of the operational flowrate range and 7.83 % or 0.431 gph at the higher end of the 
operational flowrate range. For analysis purposes, the volumetric flowrate error estimate 
was 11.8 % or 0.431 gph, whichever was smaller for the volumetric flowrate value in 
question. Basically, this means that the 11.8 % error value was applied to the lower end 
of the operational flowrate range, while the 0.431 gph error value was applied to the 
higher end of the operational flowrate range. 

The second partial derivative in Equations E.6 had to be multiplied by the 
variance of the coolant density. The precision error associated with the coolant density 
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came from the standard deviation of the inlet temperature to the test section during a CHF 
experiment. The maximum inlet temperature standard deviation encountered over all 
CHF experiments was 3.070°C (From Table C.l), which occurred on an experiment that 
had an average inlet temperature of 61.9°C and an average absolute inlet pressure of 
793.9 kPa. This inlet temperature standard deviation thus created a range of possible 
inlet densities: 980.7 -> 983.9 kg/m 3 . The density value obtained when the average 
values of temperature and pressure were used was 982.3 kg/m 3 . Thus, the representative 
deviation in density from the reported value for all CHF experiments was 1.6 kg/m 3 , 
which corresponds to a u p _i of 0.00163 or 0.163 %. This value came from the standard 
deviation of the inlet temperature, which depended on the bias error of the inlet 
thermocouple (+ 0.1 °C). The lowest inlet temperature encountered during all CHF 
experiments was 29.8°C, which produced a u p -2 of 0.00336 or 0.336 %. By combining 
these two terms as shown in Equation E.7, the final value of u p was 0.00373 or 0.373 %, 
which means that the standard deviation of the inlet temperature contributes very little to 
the error associated with the coolant density at the test section inlet. 

u p = Vvi' + V^ 2 = °- 00373 = 0.373% (E.7) 

Referring to Equation E.5, all values were known at this point such that the 
variance of the coolant mass flowrate could be determined. Since the value of 
importance to this investigation was coolant mass flux, the variance of the test section 
cross sectional area had to be determined for the error analysis to be complete, as shown 
in Equation E.2. The cross sectional area or flow area of the test section was a function 
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of the inner diameter of the outer heater (outer wall of the annulus) and the outer diameter 
of the inner heater (inner wall of the annulus). 


4,-j K-of.) 


(E.8) 


The variance of the test section cross sectional area can be expressed in the following 
manner. 



at 


V 


dD, 


Unut 


out J 


(aO 2 

\ dD inj 




(E.9) 


where: 


at 

dD m 


= £d 


and = -~D in (E.10) 

dD ; . 2 ' 


The outer diameter (inner diameter of the outer heater) was measured using both 
precision ground pins and a Starrett dial caliper. The tolerance of the caliper was known 
to be a n = 0.001 inch = 0.0254 mm, which created a maximum bias uncertainty of 

u 0 ut 

u n = 0.001/0.304 = 0.00329 = 0.329 % for the outer heater used (Outer Heater K). 

u 0 ut 

This tolerance was used instead of the tolerance of the pins to provide a more 
conservative estimate of the error associated with the cross sectional area of the test 
section. The inner diameter (outer diameter of the inner heater) was also measured using 
a Starrett dial caliper ( <r D = 0.001 inch = 0.0254 mm). The inner heater was constructed 

of tubing that had been manufactured with a nominal outer diameter of 0.250 inch. 
During construction, the inner heater was found to have a diameter slightly larger than 
0.250 inch, but 0.250 inch was utilized in this error analysis to maintain a conservative 
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estimate. The bias uncertainty estimate of the inner diameter was found to be u D = 


0.001/0.250 = 0.00400 0.400 %. 

Finally, combining Equations E.2 through E.10 yielded the following estimate of 
the coolant mass flux variance. 


“g = 


1 


L cs J 




Jt 


D. 


+ 


£ d) 

J 




(E. 11) 


By applying this formula along with the appropriate error estimates from Table E.2 to the 
set of experiments performed in this investigation, the average uncertainty in mass flux 
was calculated to be 11.08 %. Note that the <j values for coolant volumetric flowrate in 
Table E.2 are given in terms of both gph and percentage since the error estimate of the 
flowrate in question depends in part on the magnitude of the flowrate. The mass flux 
error analysis results for the data with heat flux ratio of unity are listed in Table E.3. 


Table E.2: Summary of Coolant Mass Flux Error Sources at CHF 



Up 

cr D 

u out 


11.8 % or 0.431 gph 

0.373 % 

0.0254 mm 

0.0254 mm 


230 






Table E.3: CHF Coolant Mass Flux Error Analysis Results 


Sigma(Vol. Flow): 
or 

Sigma(Density): 

Sigma(Dj)ut): 

Sigma(D in): 

4.53E-07 

11.8 

0.373 

2.54E-05 

2.54E-05 

(m A 3/sec) 

(%) 

(%) 

(m) 

(m) 


D out ~ 

D in = 

A cs = 

7.722 

6.401 

1.465E-05 

(m) 

(m) 

(m A 2) 

Uncertainty_G 

(%) 

Setup 

Run 

Inlet Density 
@CHF 
(kg/m A 3) 

Inlet Density 
sigma 
(kg/m A 3) 

Vol Flowrate 
@CHF 
(m A 3/sec) 

Vol Flowrate 
sigma 
(m A 3/sec) 

Mass 

Flowrate 

(kg/s) 

G(chf) Sigma_G 

(kg/m A 2 s) (kg/m A 2 s) 

13 

1 

996 

3.72 

3.07E-06 

3.62E-07 

3.06E-03 

207.74 

25.27 

12.16 

13 

2 

989 

3.69 

2.67E-06 

3.15E-07 

2.64E-03 

179.38 

21.81 

12.16 

13 

3 

981 

3.66 

2.81E-06 

3.31E-07 

2.75E-03 

187.34 

22.78 

12.16 

13 

4 

988 

3.69 

2.96E-06 

3.49E-07 

2.92E-03 

198.86 

24.19 

12.16 

13 

5 

988 

3.68 

2.79E-06 

3.29E-07 

2.76E-03 

187.63 

22.81 

12.16 

13 

6 

996 

3.72 

1.57E-06 

1.85E-07 

1.56E-03 

106.28 

12.93 

12.17 

13 

7 

988 

3.69 

1.50E-06 

1.77E-07 

1.48E-03 

100.68 

12.24 

12.16 

13 

8 

981 

3.66 

1.49E-06 

1.75E-07 

1.46E-03 

99.19 

12.06 

12.16 

13 

9 

988 

3.69 

1.52E-06 

1.79E-07 

1.50E-03 

102.01 

12.40 

12.16 

13 

10 

988 

3.69 

1.48E-06 

1.74E-07 

1.46E-03 

99.42 

12.08 

12.15 

13 

11 

996 

3.71 

4.30E-06 

4.53E-07 

4.28E-03 

291.15 

31.83 

10.93 

13 

12 

988 

3.68 

4.32E-06 

4.53E-07 

4.27E-03 

290.25 

31.60 

10.89 

13 

13 

980 

3.66 

4.40E-06 

4.53E-07 

4.31E-03 

293.57 

31.39 

10.69 

13 

14 

987 

3.68 

4.39E-06 

4.53E-07 

4.33E-03 

294.68 

31.59 

10.72 

13 

15 

988 

3.69 

4.43E-06 

4.53E-07 

4.38E-03 

297.63 

31.66 

10.64 

13 

16 

996 

3.72 

2.23E-06 

2.63E-07 

2.22E-03 

151.16 

18.38 

12.16 

13 

17 

989 

3.69 

2.18E-06 

2.57E-07 

2.16E-03 

146.59 

17.83 

12.16 

13 

18 

981 

3.66 

2.26E-06 

2.67E-07 

2.22E-03 

150.97 

18.36 

12.16 

13 

19 

989 

3.69 

2.28E-06 

2.69E-07 

2.25E-03 

153.14 

18.62 

12.16 

13 

20 

988 

3.69 

2.18E-06 

2.58E-07 

2.16E-03 

146.74 

17.85 

12.16 

13 

21 

997 

3.72 

3.66E-06 

4.32E-07 

3.65E-03 

248.22 

30.18 

12.16 

13 

22 

989 

3.69 

3.67E-06 

4.33E-07 

3.63E-03 

246.55 

29.99 

12.16 

13 

23 

980 

3.65 

3.63E-06 

4.28E-07 

3.55E-03 

241.67 

29.40 

12.17 

13 

24 

988 

3.69 

3.65E-06 

4.31E-07 

3.61E-03 

245.44 

29.84 

12.16 





Table E.3 (continued) 


Setup 

Run 

Inlet Density 
@OFI 
(kg/m A 3) 

Inlet Density 
sigma 
(kg/m A 3) 

Vol Flowrate 
@OFI 
(m A 3/sec) 

Vol Flowrate 
sigma 
(m A 3/sec) 

Mass 

Flowrate 

(kg/s) 

G(ofi) 

(kg/m A 2 s) 

SigmaG 

(kg/m A 2 s) 

UncertaintyG 

(%) 

13 

25 

989 

3.69 

3.68E-06 

4.34E-07 

3.64E-03 

247.55 

30.09 

12.16 

13 

26 

994 

3.71 

5.00E-06 

4.53E-07 

4.98E-03 

338.46 

32.16 

9.50 

13 

27 

986 

3.68 

4.94E-06 

4.53E-07 

4.88E-03 

331.58 

31.85 

9.61 

13 

28 

981 

3.66 

4.93E-06 

4.53E-07 

4.84E-03 

328.99 

31.70 

9.63 

13 

29 

986 

3.68 

5.06E-06 

4.53E-07 

4.99E-03 

339.69 

31.93 

9.40 

13 

30 

987 

3.68 

5.03E-06 

4.53E-07 

4.97E-03 

338.09 

31.94 

9.45 

13 

31 

994 

3.71 

5.55E-06 

4.53E-07 

5.52E-03 

375.14 

32.45 

8.65 

13 

32 

987 

3.68 

5.53E-06 

4.53E-07 

5.46E-03 

371.22 

32.22 

8.68 

13 

33 

981 

3.66 

5.50E-06 

4.53E-07 

5.39E-03 

366.80 

31.99 

8.72 

13 

34 

987 

3.68 

5.57E-06 

4.53E-07 

5.49E-03 

373.54 

32.23 

8.63 

13 

35 

988 

3.68 

5.65E-06 

4.53E-07 

5.58E-03 

379.39 

32.31 

8.52 


N) 

K> 



The heat flux can be calculated as a function of the power applied to the test 


section and the test section heated area. 

Q"=9l2!SL (E.12) 

A 

The variance of the surface heat flux can be expressed by the following error propagation 
formula. 



(E.13) 


where 


dor 

total 



dQ" 

dA s 



(El 4) 


For all 35 CHF experiments, the power supplied to the test section was 
determined using an energy balance and the total amount of power supplied by the DC 
power supply as described in Section IV.8.2. Errors in the energy balance heat loss 
estimation were also discussed independently in Section IV.8.2, however, the error 
associated with the measurement of the total power output from the power supply had to 
be quantified. This power level was measured using the Ohio Semitronics DC wattmeter 
(Component 16 of Table 3.3) mentioned earlier. During CHF experimentation, there 
were two individual bias errors that contributed to the uncertainty of the total power 


reading: 

1. The manufacturer stated accuracy of the Ohio Semitronics DC wattmeter was 
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±1.0% of the power reading, thus Ui = 0.01 = 1.0%. 

2. The DAS voltage conversion error was calculated by determining the bit range 
over which the wattmeter was used. The calibration curve for this DC wattmeter 
was the same as stated in Appendix D: 

Power { W) = -2.31 W + 6.4097(mF £W>16 )— (E ‘ 15) 

mV EXP16 

with EXP 16 gain of unity. Since the gain was unity, the mV reading on the DAS- 
8 board was the same as on the EXP 16 board. The conversion from millivolts to 
bits was 2047 bits = 5000 mV, accurate to ±1 bit. The total power level over all 
35 CHF experiments ranged from a minimum of 1918 W to a maximum of 8879 
W, which converted to a range of 123 -> 568 bits. Therefore, for the 1918 W 
power level, U 2 -low = 1/123 = 0.00813 = 0.813 % and for the 8879 W power level, 
u 2 -Hi g h = 1/568 = 0.00176 = 0.176 %. 


The only precision error encountered in this investigation with respect to the total 
applied power level came directly from the standard deviation of the power level around 
CHF. The maximum standard deviation of the total applied power was 55.7 W (0.0067 
MW/m ), which occurred during an experiment where CHF occurred at a power level of 
3895.6 W (0.4689 MW/m 2 ), therefore, u 5 = 55.7/3895.6 = 0.0143 = 1.43 %. 

By summing the squares of the two bias errors and the single precision error 
associated with the total power level at CHF, the error estimate for the total power level 
at CHF was found. 
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Table E.4: Summary of Total Power Level Uncertainties at CHF 



■SSS9 


Uncertainty 2 

Uncertainty 

Uncertainty 2 




CHEnisn 

| 

1 

m 8379 W) 

Ui 

Bias 

wmmm 

1.000 E-4 

1.0 % 

1.000 E-4 

U2 

Bias 

0.813 % 



3.098 E-6 

U3 

Precision 

1.43 % 



2.045 E -4 

Sum 

— 

— 

3.706 E -4 

— 

3.076 E -4 


Using this analysis procedure, the error was 1.93 % or 37.0 W at the lower end of 
the operational power range and 1.75 % or 155.4 W at the higher end of the operational 
power range. For analysis purposes, the total power error estimate was 1.93 % or 155.4 
W, whichever was smaller for the total power value in question. Basically, this means 
that the 1.93 % error value was applied to the lower end of the operational power range, 
while the 155.4 W error value was applied to the higher end of the operational power 
range. 


The test section heated surface area was a function of the inner diameter of the 
annulus, outer diameter of the annulus, inner heater heated length, and outer heater heated 
length. 


A- s Tr^Lfeated,out^'out ^heated,in^in) 


(E. 16) 


The variance of the test section heated surface area can be expressed in the following 


manner. 
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(E.17) 


where 
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As stated earlier, the relevant annulus diameters were measured using a Starrett dial 
caliper, which had a known tolerance of a = 0.001 inch. The heated lengths of all test 
sections were measured using Fowler vernier calipers, also with a known tolerance of a = 
0.001 inch. 

Finally, combining Equations E.13 through E.18 yielded the following estimate of 
the surface heat flux variance at CHF. 




r) k.+°eJ+ 


r ! 


-j) li^o J +(^) 2 <„, u , +(xDjcr 2 Lh J 


(E.19) 


By applying this formula along with the appropriate error estimates from Table E.5 to the 
set of experiments performed in this investigation, the average uncertainty in heat flux 
was calculated to be 1.92 %. Note that the a values for all power error estimates in Table 
E.5 are given in terms of both Watts and percentage since the error estimate of the power 
level in question depends on the magnitude of the power measurement. 


Table E.5: Summary of Heat Flux Uncertainties at CHF 


a Q<o,*> 

&D 

U OUl 

^D,„ 

CTr 

^ A.our 

cr, 

L h.ht 

1.93 % or 155.4 W 

0.0254 mm 

0.0254 mm 

0.0254 mm 

0.0254 mm 
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The outline below shows a sample error estimate for a CHF experiment that had 
the following experimentally measured parameters: D out = 7.722 mm, Dj n = 6.401 mm, 
Lh, out =185 mm, Lh, j n = 190 mm, and Q to tai = 8878 W. The heat flux error analysis 
results for the data are listed in Table E.6 


Measurement 

Uncertainty 

Influence Coefficient, ( Q q " y 

J 

(Uncertainty) 2 

Basis 

Total Power 

155.4 W 

■■DHI 

3.50 E -4 (MW/m 2 ) 2 

Resolution 

■■SSExtESI 

Outer Diameter 

0.0254 mm 

Q total ^ out 

3.60 E -6 (MW/m 2 ) 2 

Measurement 

out Lout + L> in L ln ') 

= -14.1 (MW/m 3 ) 

Inner Diameter 

0.0254 mm 

Q total ^ in 

3.80 E-6 (MW/m 2 ) 2 

Measurement 

*(■ D ou, L ou , + D m L in ) 2 
= -76.7 (MW/m 3 ) 

Outer L h 

0.0254 mm 

i 

Qtotai ^out 

6.28 E-9 (MW/m 2 ) 2 

Measurement 

r( D output +D in L tt,) 2 

= -3.12 (MW/m 3 ) 

Inner L h 

0.0254 mm 

Qtotai ^in 


Measurement 

x( D ou, L out + D in L in )" 

= -2.59 (MW/m 3 ) 



Sum of Variances = 

3.50 E -4 (MW/m 2 ) 2 




Uncertainty = 

0.0187 MW/m 2 = 1.75% 
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Table E.6: CHF Surface Heat Flux Error Analysis Results 


Sigma(Power): 

155.4 

(W) 

or 

1.93 

(%) 

Sigma(L_h, out): 

2.54E-05 

(m) 

Sigma(L_h, in): 

2.54E-05 

(m) 

Sigma(D_out): 

2.54E-05 

(m) 

Sigma(D in): 

2.54E-05 

(m) 


D_out = 

7.722 

(mm) 

L_h, out = 

185 

(mm) 

D_in - 

6.401 

(mm) 

L_h, in = 

190 

(mm) 

A s = 

8.308E-03 

(m A 2) 


Setup 

Run 

Total 
Power 
(W) ... 

Qtotal 

sigma 

(W) 

CT(chf) 

(MW/m A 2) 

Sigma_Q M 

(MW/m A 2) 

Sigma_Q 

(%> 

13 

1 

3896 

75 

0.4689 

9.05E-03 

1.93 

13 

2 

3663 

71 

0.4409 

8.51E-03 

1.93 

13 

3 

3608 

70 

0.4343 

8.38E-03 

1.93 

13 

4 

3274 

63 

0.3940 

7.60E-03 

1.93 

13 

5 

4143 

80 

0.4986 

9.62E-03 

1.93 

13 

6 

2631 

51 

0.3167 

6.11E-03 

1.93 

13 

7 

2181 

42 

0.2625 

5.07E-03 

1.93 

13 

8 

1918 

37 

0.2309 

4.46E-03 

1.93 

13 

9 

1987 

38 

0.2392 

4.62E-03 

1.93 

13 

10 

2996 

58 

0.3606 

6.96E-03 

1.93 

13 

11 

6373 

123 

0.7671 

1.48E-02 

1.93 

13 

12 

6371 

123 

0.7668 

1.48E-02 

1.93 

13 

13 

4437 

86 

0.5340 

1.03E-02 

1.93 

13 

14 

5159 

100 

0.6209 

1.20E-02 

1.93 

13 

15 

6650 

128 

0.8004 

1.54E-02 

1.93 

13 

16 

3288 

63 

0.3958 

7.64E-03 

1.93 

13 

17 

2811 

54 

0.3383 

6.53E-03 

1.93 

13 

18 

2841 

55 

0.3419 

6.60E-03 

1.93 

13 

19 

2518 

49 

0.3031 

5.85E-03 

1.93 

13 

20 

3291 

64 

0.3961 

7.64E-03 

1.93 

13 

21 

5336 

103 

0.6422 

1.24E-02 

1.93 

13 

22 

4834 

93 

0.5818 

1.12E-02 

1.93 

13 

23 

4497 

87 

0.5412 

1.04E-02 

1.93 

13 

24 

4405 

85 

0.5302 

1.02E-02 

1.93 
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Table E.6 (continued) 


Setup 

Run 

Total 

Power 

(W) 

Qtotal 

sigma 

(W) 

Q"(chf) 

(MW/m A 2) 

SigmaQ" 

(MW/m A 2) 

SigmaQ 

(%) 

13 

25 

5705 

110 

0.6867 

1.33E-02 

1.93 

13 

26 

7169 

138 

0.8629 

1.67E-02 

1.93 

13 

27 

6936 

134 

0.8348 

1.61E-02 

1.93 

13 

28 

6120 

118 

0.7366 

1.42E-02 

1.93 

13 

29 

6411 

124 

0.7716 

1.49E-02 

1.93 

13 

30 

8426 

155 

1.0142 

1.87E-02 

1.84 

13 

31 

8868 

155 

1.0674 

1.87E-02 

1.75 

13 

32 

7688 

148 

0.9253 

1.79E-02 

1.93 

13 

33 

6538 

126 

0.7869 

1.52E-02 

1.93 

13 

34 

6469 

125 

0.7786 

1.50E-02 

1.93 

13 

35 

8878 

155 

1.0686 

1.87E-02 

1.75 
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